Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2008-03-30, 21:42   #1
cckk0508
初级会员
 
注册日期: 2008-03-20
年龄: 38
帖子: 1
声望力: 0
cckk0508 正向着好的方向发展
默认 [求助]求助while问题

各位大侠帮我看看我的程序哪有问题!
clc;
clear;
syms Kd Ka n1 n2 y1 y2;
X=[Kd Ka];
N=[n1 n2];
y=[y1 y2];
Z=[10+20*X(1)/(X(2)-X(1))*(exp(-X(1)*8/4)-exp(-X(2)*8/4))-8.5]^2+[10+20*X(1)/(X(2)-X(1))*(exp(-X(1)*28/4)-exp(-X(2)*28/4))-7]^2+[10+20*X(1)/(X(2)-X(1))*(exp(-X(1)*36/4)-exp(-X(2)*36/4))-6.1]^2+[10+20*X(1)/(X(2)-X(1))*(exp(-X(1)*56/4)-exp(-X(2)*56/4))-7.2]^2;
x=[1 2];
old=subs(Z,[X(1) X(2)],[x(1) x(2)]);
Z1=[diff(Z,X(1));diff(Z,X(2))];
H=[diff(Z,X(1),2) diff(Z1(1),X(2));diff(Z1(2),X(1)) diff(Z,X(2),2)];
step=Z1.'*Z1\(Z1.'*H*Z1);
F=subs(step,[X(1) X(2)],[x(1) x(2)]);
N(1)=x(1)-F*subs(Z1(1),[X(1) X(2)],[x(1) x(2)]);
N(2)=x(2)-F*subs(Z1(2),[X(1) X(2)],[x(1) x(2)]);
new=subs(Z,[X(1) X(2)],[N(1) N(2)]);
while (abs(new-old)/old)>exp(-6),
x(1)=N(1);
x(2)=N(2);
old=new;
F=subs(step,[X(1) X(2)],[x(1) x(2)]);
N(1)=x(1)-F*subs(Z1(1),[X(1) X(2)],[x(1) x(2)]);
N(2)=x(2)-F*subs(Z1(2),[X(1) X(2)],[x(1) x(2)]);
new=subs(Z,[X(1) X(2)],[N(1) N(2)]);
end
X(1)=N(1);
X(2)=N(2);
为什么总说??? Error using ==> >=
Function '>=' is not defined for values of class 'sym'.
要怎么改改啊?谢谢各位!
cckk0508 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码


相似的主题
主题 主题作者 版面 回复 最后发表
BMP图像转换成时间序列 zlnhappy MATLAB论坛 0 2008-04-28 17:53
matlab编写板卡驱动 emirlee MATLAB论坛 2 2008-04-18 09:16
【求助】如何做FFT变换? ceeshy MATLAB论坛 1 2007-12-21 11:27
求救-----基于小波域的数字水印 tcgqb111 MATLAB论坛 0 2007-05-24 07:27


所有时间均为北京时间。现在的时间是 12:02


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.