 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 453 | 回復: 1 | ||
bwklkx銅蟲 (小有名氣)
|
[求助]
關于非線性方程 已解決
|
|
用matlab解非線性方程的時候建立了兩個m文件。 ———————————————————————————————————————— 以下是兩個自定義的m文件。 %二分法求解方程 function [x,k]=demimethod(a,b,f,emg) % a,b :求解區(qū)間的兩個端點 % f :所求方程的函數(shù)名 % emg ：精度指標 % x：所求近似解 % k: 循環(huán)次數(shù) fa=feval(f,a); fab=feval(f,(a+b)/2); k=0; while abs(b-a)>emg if fab==0 x=(a+b)/2; return; elseif fa*fab<0 b=(a+b)/2; else a=(a+b)/2; end fa=feval(f,a); fab=feval(f,(a+b)/2); k=k+1; end x=(a+b)/2; ------------------------------------------------------------------------------------------------------ 第二個文件，是建立函數(shù)的文件。 function fxz1=abcd(t) xz=t K11=14.8; K33=15.8; K=(K11-K33)/K33; E0=8.85; %空氣的介電常數(shù) Ev=9.72; %垂直介電常數(shù) Ep=4.11; %平行介電常數(shù) E=Ev-Ep; %介電各向異性 r=E/Ep; Vth=pi*sqrt(K33/(E*E0)); %閾值電壓 V=10; %施加電壓 d=10000; %液晶盒厚度單位nm no=1.484; ne=1.636; xo=0; l=532;%入射光波長，nm ZD=2;%衍射屏距觀察屏距離，m xm=1.5707; X=(sin(xm))^2; sxm1=12.665; a=xo;%此處修改預傾角 b=xz; n1=1000;%此處修改多少小區(qū)間數(shù)目 % n1,n分別為小區(qū)間、大區(qū)間的個數(shù)。 % h為其步長。 % s1，s2為函數(shù)在大區(qū)間、小區(qū)間的端點的函數(shù)值之和。 n=n1/2; h=(b-a)/n; for i=1:n Xk(i)=a+(i-1)*h; end s1=0; for i=2:n f1(i)=1/(Xk(i)*Xk(i)-1); s1=s1+f1(i); end s2=0; for i=1:n Xl(i)=Xk(i)+h/2; f2(i)=1/(Xl(i)*Xl(i)-1); s2=s2+f2(i); end sz=((b-a)/6)*(sqrt((1+r*X*sin(a)^2)*(1+K*X*sin(a)^2)/(1-X*sin(a)^2)... +4*s2+2*s1... +sqrt((1+r*X*sin(b)^2)*(1+K*X*sin(b)^2)/(1-X*sin(b)^2)))); oo=(0:0.1:0.5); for i=1:6 fxz1=0.5*sz/sxm1-oo(i); end ——————————————————————————————————————— 下一步是執(zhí)行的命令：t=demimethod(-pi/2,pi/2,'abcd',10^-7); 想要得到的是t的6個解，但是運算得出的只是一個解。運算出的是fxz1=0.5*sz/sxm1-oo(6)的解。請問如何修改文件，能夠解出fxz1=0.5*sz/sxm1-oo(i)所有循環(huán)的解，而不只是最后一個循環(huán)的解。 已解決，謝謝。 [ Last edited by bwklkx on 2011-6-30 at 20:53 ] |
回復此樓» 猜你喜歡
0856調(diào)劑，是學校就去
已經(jīng)有7人回復
-
0817 化學工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有14人回復
-
申博26年
已經(jīng)有3人回復
-
一志愿吉林大學材料學碩321求調(diào)劑
已經(jīng)有10人回復
-
0703化學調(diào)劑 ，六級已過，有科研經(jīng)歷
已經(jīng)有12人回復
-
294求調(diào)劑材料與化工專碩
已經(jīng)有13人回復
-
288求調(diào)劑
已經(jīng)有14人回復
-
307求調(diào)劑
已經(jīng)有8人回復
-
一志愿南京理工大學085701資源與環(huán)境302分求調(diào)劑
已經(jīng)有5人回復
-
288求調(diào)劑，一志愿華南理工大學071005
已經(jīng)有5人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- 求高人指點用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復
- 【求助】MATLAB 解方程組【已解決】
已經(jīng)有7人回復
- 【求助】用fortran求解大型線性方程組時出現(xiàn)的錯誤【已解決】
已經(jīng)有11人回復
- 【求助】使用MATLAB怎么實現(xiàn)擬合力場參數(shù)的程序？【已解決】
已經(jīng)有15人回復
- 【求助】線代的特征值與特征向量【已解決】
已經(jīng)有8人回復
- 【求助】求解非線性方程
已經(jīng)有8人回復
- 隨機非線性微分方程如何離散化??
已經(jīng)有8人回復
- 用SAS進行非線性擬合
已經(jīng)有8人回復
 |
2樓2011-07-06 14:32:15
相關版塊跳轉(zhuǎn)
我要訂閱樓主 bwklkx 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招收調(diào)劑碩士 +4 | lidianxing 2026-03-19 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 一志愿西安交通大學材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
|
[考研] 281求調(diào)劑（0805） +9 | 煙汐憶海 2026-03-16 | 19/950 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +7 | 小材化本科 2026-03-18 | 7/350 |
|
|
|
[考研] 一志愿985，本科211，0817化學工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 0703化學調(diào)劑 +8 | 妮妮ninicgb 2026-03-15 | 12/600 |
|
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
|
[考研] 化學工程321分求調(diào)劑 +15 | 大米飯！ 2026-03-15 | 18/900 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠意在機?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 268求調(diào)劑 +7 | 好運連綿不絕 2026-03-12 | 8/400 |
|
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿蘇州大學材料工程（085601）專碩有科研經(jīng)歷三項國獎兩個實用型專利一項省級立項 +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 材料與化工專碩調(diào)劑 +5 | heming3743 2026-03-16 | 5/250 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|