Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2009-06-04, 15:09   #1
alee817
初级会员
 
注册日期: 2009-05-25
帖子: 4
声望力: 0
alee817 正向着好的方向发展
默认 这个有点小错~可是却不知道怎么改,麻烦高手

clear;clc;
x=[-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180 210 240 270 300 330];
y=[0 0.02 0.03 0.05 0.07 0.1 0.15 0.20 0.26 0.30 0.28 0.25 0.20 0.16 0.1 0.05 0.03 0.01];
myfunc=inline('-2*pi*beta(1)*0.33/31684(x*exp(-pi*x*x/31684)-(x-204.08)*exp(-pi*(x-204.08)*(x-204.08)/31684))','beta','x');
beta=nlinfit(x,y,myfunc,[0.2]);
b=beta(1)
%test the model
xx=min(x):max(x);
yy=-2*pi*b*0.33/31684(x*exp(-pi*x*x/31684)-(x-204.08)*exp(-pi*(x-204.08)*(x-204.08)/31684));
plot(x,y,'o',x,yy,'r')
alee817 当前离线   回复时引用此帖
旧 2009-06-04, 16:13   #2
TTT_IOU
普通会员
 
注册日期: 2009-02-22
帖子: 91
声望力: 18
TTT_IOU 正向着好的方向发展
默认 回复: 这个有点小错~可是却不知道怎么改,麻烦高手

clear;clc;
x=[-180 -150 -120 -90 -60 -30 0 30 60 90 120 150 180 210 240 270 300 330];
y=[0 0.02 0.03 0.05 0.07 0.1 0.15 0.20 0.26 0.30 0.28 0.25 0.20 0.16 0.1 0.05 0.03 0.01];
myfunc=inline('-2*pi*beta*0.33/31684*(x.*exp(-pi*x.*x/31684)-(x-204.08).*exp(-pi*(x-204.08).*(x-204.08)/31684))','beta','x');
beta=nlinfit(x,y,myfunc,0.2);
b=beta
%test the model
xx=min(x):max(x);
yy=-2*pi*b*0.33/31684*(x.*exp(-pi*x.*x/31684)-(x-204.08).*exp(-pi*(x-204.08).*(x-204.08)/31684));
plot(x,y,'o',x,yy,'r')
inline的函数应和yy的一样,且 要注意矩阵相应项相乘时应加‘.’,望有帮助!!!
TTT_IOU 当前离线   回复时引用此帖
回复

主题工具
显示模式

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

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



所有时间均为北京时间。现在的时间是 00:06


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