版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

bwklkx

銅蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 348
帖子: 71
在線: 32.3小時(shí)
蟲號(hào): 1138356
注冊(cè): 2010-11-04
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 關(guān)于非線性方程已解決

用matlab解非線性方程的時(shí)候建立了兩個(gè)m文件。
————————————————————————————————————————
以下是兩個(gè)自定義的m文件。

%二分法求解方程
function [x,k]=demimethod(a,b,f,emg)
% a,b :求解區(qū)間的兩個(gè)端點(diǎn)
% f :所求方程的函數(shù)名
% emg ：精度指標(biāo)
% 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;

------------------------------------------------------------------------------------------------------
第二個(gè)文件，是建立函數(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;%入射光波長(zhǎng)，nm
ZD=2;%衍射屏距觀察屏距離，m
xm=1.5707;
X=(sin(xm))^2;
sxm1=12.665;

a=xo;%此處修改預(yù)傾角
b=xz;
n1=1000;%此處修改多少小區(qū)間數(shù)目
% n1,n分別為小區(qū)間、大區(qū)間的個(gè)數(shù)。
% h為其步長(zhǎng)。
% s1，s2為函數(shù)在大區(qū)間、小區(qū)間的端點(diǎn)的函數(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個(gè)解，但是運(yùn)算得出的只是一個(gè)解。運(yùn)算出的是fxz1=0.5*sz/sxm1-oo(6)的解。請(qǐng)問如何修改文件，能夠解出fxz1=0.5*sz/sxm1-oo(i)所有循環(huán)的解，而不只是最后一個(gè)循環(huán)的解。

已解決，謝謝。

[ Last edited by bwklkx on 2011-6-30 at 20:53 ]

回復(fù)此樓

» 猜你喜歡

0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有17人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有5人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有7人回復(fù)
申博26年已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有10人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過，有科研經(jīng)歷已經(jīng)有12人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有13人回復(fù)
288求調(diào)劑已經(jīng)有14人回復(fù)
307求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑已經(jīng)有5人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；已經(jīng)有11人回復(fù)
【求助】MATLAB 解方程組【已解決】已經(jīng)有7人回復(fù)
【求助】用fortran求解大型線性方程組時(shí)出現(xiàn)的錯(cuò)誤【已解決】已經(jīng)有11人回復(fù)
【求助】使用MATLAB怎么實(shí)現(xiàn)擬合力場(chǎng)參數(shù)的程序？【已解決】已經(jīng)有15人回復(fù)
【求助】線代的特征值與特征向量【已解決】已經(jīng)有8人回復(fù)
【求助】求解非線性方程已經(jīng)有8人回復(fù)
隨機(jī)非線性微分方程如何離散化?? 已經(jīng)有8人回復(fù)
用SAS進(jìn)行非線性擬合已經(jīng)有8人回復(fù)

1樓 2011-06-29 22:08:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

mengxc

2樓2011-07-06 14:32:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 bwklkx 的主題更新