24小时热门版块排行榜    

查看: 590  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 00:55 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[教师之家] 跳槽后在研项目怎么办? +3 简单化xn 2026-08-22 7/350 2026-08-22 19:46 by 简单化xn
[基金申请] 什么时候开奖? +7 CrisMessi 2026-08-18 8/400 2026-08-22 18:29 by 淀粉搬运工
[基金申请] 只有每年这种时候来逛逛小木虫 +24 yaoyewhu2008 2026-08-20 26/1300 2026-08-22 17:43 by kammury
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 今天系统多次维护,明天很可能放榜! +10 zju2000 2026-08-16 11/550 2026-08-20 20:02 by cl479861084
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见