版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

bwklkx

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 348
帖子: 71
在线: 32.3小时
虫号: 1138356
注册: 2010-11-04
性别: GG
专业: 凝聚态物性 II ：电子结构

[求助] 关于非线性方程已解决

用matlab解非线性方程的时候建立了两个m文件。
————————————————————————————————————————
以下是两个自定义的m文件。

%二分法求解方程
function [x,k]=demimethod(a,b,f,emg)
% a,b :求解区间的两个端点
% f :所求方程的函数名
% emg ：精度指标
% x：所求近似解
% k: 循环次数
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;

------------------------------------------------------------------------------------------------------
第二个文件，是建立函数的文件。
function fxz1=abcd(t)

xz=t

K11=14.8;
K33=15.8;
K=(K11-K33)/K33;
E0=8.85; %空气的介电常数
Ev=9.72; %垂直介电常数
Ep=4.11; %平行介电常数
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;%此处修改多少小区间数目
% n1,n分别为小区间、大区间的个数。
% h为其步长。
% s1，s2为函数在大区间、小区间的端点的函数值之和。
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
———————————————————————————————————————

下一步是执行的命令：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)所有循环的解，而不只是最后一个循环的解。

已解决，谢谢。

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

回复此楼

» 猜你喜欢

售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有4人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有4人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有5人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有5人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有5人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有4人回复
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
【求助】MATLAB 解方程组【已解决】已经有7人回复
【求助】用fortran求解大型线性方程组时出现的错误【已解决】已经有11人回复
【求助】使用MATLAB怎么实现拟合力场参数的程序？【已解决】已经有15人回复
【求助】线代的特征值与特征向量【已解决】已经有8人回复
【求助】求解非线性方程已经有8人回复
随机非线性微分方程如何离散化?? 已经有8人回复
用SAS进行非线性拟合已经有8人回复

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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

mengxc

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