24小时热门版块排行榜    

查看: 490  |  回复: 1

樱木花道6

木虫 (知名作家)

[求助] 急:求助MATLAB拟合问题

function f=Dielectric1_1(c,fdata)
f=c(1)+c(2)./(1+c(3).*fdata.*fdata);

clc;clear;
fdata=[
8.23E8
8.57E8
8.91E8
9.24E8
9.58E8
9.92E8
1.03E9
1.07E9
1.11E9
1.15E9
1.2E9
1.24E9
1.29E9
1.33E9
1.38E9
1.44E9
1.5E9];
ydata=[
13.20549
12.96021
12.66066
12.37507
12.0785
11.75602
11.45932
11.12179
10.76602
10.33523
9.95384
9.46352
8.91796
8.369
7.63211
6.83466
6.01964];
c0=[10000 1000 10];
c=lsqcurvefit('Dielectric1_1',c0,fdata,ydata);
y=Dielectric1_1(c,fdata);
%x1=c(1)*c(2)*c(3)*c(3)-c(1)*c(2)*c(2)*c(3)*c(3)-2*c(1)*c(2)*c(3)*c(4);
%x2=(c(2)*c(3)+c(4)).^2;
%x3=3.088e-21*c(1)*c(1)*c(3)*c(3)*c(4)*c(4)*1e9*1e9;
plot(fdata,ydata,'o');hold on;plot(fdata,y,'r');grid on
回复此楼
NoPain,Nogain
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)


dbb627(金币+1): 欢迎交流 2012-03-05 15:39:56
用1stOpt吧,Matlab对初值依赖严重,猜初值是个问题:
CODE:
Variable fdata,ydata;
Function ydata=c1+c2/(1+c3*fdata^2);
Data;
8.23E8        13.20549
8.57E8        12.96021
8.91E8        12.66066
9.24E8        12.37507
9.58E8        12.0785
9.92E8        11.75602
1.03E9        11.45932
1.07E9        11.12179
1.11E9        10.76602
1.15E9        10.33523
1.2E9        9.95384
1.24E9        9.46352
1.29E9        8.91796
1.33E9        8.369
1.38E9        7.63211
1.44E9        6.83466
1.5E9        6.01964

均方差(RMSE): 0.0579937148077891
残差平方和(SSE): 0.0571756062725221
相关系数(R): 0.999627623709406
相关系数之平方(R^2): 0.999255386082914
决定系数(DC): 0.999255386082914
卡方系数(Chi-Square): 0.0028094926579597
F统计(F-Statistic): 9400.8449739933

参数        最佳估算
----------        -------------
c1        82.4460672941467
c2        -66.6397980875574
c3        -5.69597187742321E-20


2楼2012-03-05 12:05:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 樱木花道6 的主题更新
信息提示
请填处理意见