24小时热门版块排行榜    

查看: 589  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 基金啊基金 +3 longfie172 2026-08-20 3/150 2026-08-21 00:04 by wuchongjun
[基金申请] 科研孤儿太难了 +14 我4大白菜 2026-08-20 15/750 2026-08-20 22:33 by 袁向阳007
[基金申请] 言之凿凿,明天官网会通知放榜? +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] filecode=后面第一个是大写字母 +10 wangze12014 2026-08-14 12/600 2026-08-19 16:56 by 苦难博士
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见