24小时热门版块排行榜    

查看: 441  |  回复: 1

菜头_

金虫 (小有名气)

[求助] 关于非线性最小二乘的拟合 大神帮忙看看哪出错了

clc
clear
X=[100000.0, 100676.7, 101353.4, 102030.1, 102706.81, 103383.51, 104060.21, 104736.91, 105413.61, 106090.31, 106767.02, 107443.72, 108120.42, 108797.12, 109473.82, 110150.53, 110827.23, 111503.93, 112180.63, 112857.33, 113534.03, 114210.74, 114887.44, 115564.14, 116240.84, 116917.54, 117594.24, 118270.95, 118947.65, 119624.35, 120301.05, 120977.75, 121654.45, 122331.16, 123007.86, 123684.56, 124361.26, 125037.96, 125714.67, 126391.37, 127068.07, 127744.77, 128421.47, 129098.17, 129774.88, 130451.58, 131128.28, 131804.98, 132581.69, 133476.12, 134370.56, 135264.99, 136159.42, 137053.85, 137948.28, 138842.71, 139737.14, 140631.57, 141526.0, 142420.44, 143314.87, 144209.3, 145103.73, 145998.16, 146892.59, 147787.02, 148681.46, 149575.89, 150470.32, 151364.75, 152259.18, 153153.61, 154048.04, 154942.47, 155836.9, 156731.33, 157625.77, 158520.2, 159414.63, 160309.06, 161203.49, 162097.92, 162992.35, 163886.79, 164781.22, 165675.65, 166570.08, 167464.51, 168358.94, 169253.37, 170147.8, 171042.23, 171936.67, 172831.1, 173725.53, 174749.27, 175928.18, 177107.08, 178285.99, 179464.89, 180643.8, 181822.7, 183001.61, 184180.51, 185359.42, 186538.32, 187717.23, 188896.13, 190075.04, 191253.94, 192432.85, 193611.75, 194790.66, 195969.56, 197148.47, 198327.37, 199506.28, 200685.18, 201864.09, 203042.99, 204221.9, 205400.8, 206579.71, 207758.61, 208937.52, 210116.42, 211295.33, 212474.23, 213653.13, 214832.04, 216010.94, 217189.85, 218368.75, 219547.66, 220726.56, 221905.47, 223084.37, 224263.28, 225442.18, 226621.09, 227799.99, 228978.9, 230328.25, 231882.1, 233435.96, 234989.81, 236543.67, 238097.52, 239651.38, 241205.23, 242759.09, 244312.94, 245866.8, 247420.65, 248974.51, 250528.37, 252082.22, 253636.08, 255189.93, 256743.79, 258297.64, 259851.5, 261405.35, 262959.21, 264513.06, 266066.92, 267620.77, 269174.63, 270728.49, 272282.34, 273836.2, 275390.05, 276943.91, 278497.76, 280051.62, 281605.47, 283159.33, 284713.18, 286267.04, 287820.89, 289374.75, 290928.6, 292482.46, 294036.31, 295590.17, 297144.03, 298697.88, 300251.74, 301805.59, 303584.1, 305632.16, 307680.22, 309728.28, 311776.33, 313824.39, 315872.45, 317920.51, 319968.57, 322016.63, 324064.69, 326112.75, 328160.8, 330208.86, 332256.92, 334304.98, 336353.04, 338401.1, 340449.16, 342497.21, 344545.27, 346593.33, 348641.39, 350689.45, 352737.51, 354785.57, 356833.63, 358881.68,];
Y=[0.0003476636, 0.0003481631, 0.0003485057, 0.0003534742, 0.0003585374, 0.0003561982, 0.0003610355, 0.0003630067, 0.0003636565, 0.0003670694, 0.000369875, 0.0003720718, 0.0003770485, 0.0003788507, 0.0003803943, 0.0003837635, 0.0003846664, 0.0003882613, 0.00038865, 0.0003941239, 0.000393931, 0.0003957557, 0.0003978632, 0.0004007074, 0.0004029052, 0.0004039185, 0.0004067978, 0.0004077394, 0.0004082662, 0.0004120398, 0.0004136187, 0.0004163448, 0.0004180682, 0.0004181161, 0.0004197053, 0.0004237, 0.0004261216, 0.0004277627, 0.0004298843, 0.0004308783, 0.0004320549, 0.0004352948, 0.0004380516, 0.0004375443, 0.0004391876, 0.0004378034, 0.0004400687, 0.0004457715, 0.0004460298, 0.0004477484, 0.0004508044, 0.0004530239, 0.0004551622, 0.0004571168, 0.0004607805, 0.0004606528, 0.000465276, 0.0004653396, 0.0004643261, 0.0004681045, 0.0004686525, 0.0004774228, 0.0004766617, 0.000476559, 0.0004810295, 0.0004795355, 0.0004829681, 0.0004881095, 0.0004890235, 0.0004847109, 0.0004886249, 0.0004961468, 0.0004948185, 0.0004939894, 0.0005026379, 0.0004926372, 0.0005003095, 0.000503053, 0.0005055127, 0.0005037699, 0.0005059391, 0.0005073837, 0.0005076136, 0.0005113738, 0.0005093064, 0.0005098524, 0.0005157921, 0.0005134944, 0.0005148812, 0.0005107459, 0.0005145137, 0.0005116848, 0.0005183074, 0.0005178733, 0.0005152463, 0.0005141246, 0.0005195564, 0.0005187643, 0.0005274088, 0.0005248351, 0.0005243403, 0.0005268489, 0.0005258592, 0.0005303587, 0.0005280853, 0.0005331685, 0.0005380624, 0.0005373088, 0.0005490789, 0.0005431627, 0.0005543247, 0.0005442015, 0.0005585522, 0.0005626062, 0.0005659603, 0.0005707157, 0.0005813651, 0.000584626, 0.0005853017, 0.0005873509, 0.0005986107, 0.0005998682, 0.0006033591, 0.000602516, 0.0006131896, 0.0006222505, 0.000617344, 0.0006252857, 0.0006253377, 0.000635688, 0.0006300505, 0.0006321831, 0.0006388351, 0.0006479359, 0.0006447796, 0.000647399, 0.0006619523, 0.0006572286, 0.0006698229, 0.0006831476, 0.0006752985, 0.0006799055, 0.0006875366, 0.0007078557, 0.0007073916, 0.0007183332, 0.0007307543, 0.0007529389, 0.0007660364, 0.000763928, 0.0007806155, 0.0007935706, 0.0008028062, 0.000825684, 0.0008268275, 0.0008333289, 0.0008539869, 0.0008515532, 0.0008631907, 0.0008709452, 0.0008648916, 0.0008826287, 0.00089975, 0.000900333, 0.0009194366, 0.000922961, 0.000930595, 0.0009425183, 0.0009477795, 0.0009618909, 0.0009701521, 0.0009774935, 0.0009913497, 0.001001847, 0.001010384, 0.00102605, 0.001034779, 0.001045183, 0.001054803, 0.001064564, 0.001078737, 0.00109312, 0.001100417, 0.001116339, 0.001135811, 0.001153047, 0.001167981, 0.0011787, 0.001194414, 0.001216686, 0.001241228, 0.001257347, 0.001282052, 0.001308651, 0.001328122, 0.001347286, 0.001366008, 0.001380992, 0.001392525, 0.001403476, 0.001391937, 0.001386341, 0.001373025, 0.001358931, 0.001335692, 0.001315536, 0.001298277, 0.001277081, 0.001262265, 0.001250303, 0.001240685, 0.001232468, 0.001227612, 0.001225629, 0.001220269, 0.001220192, 0.00122092,];
f=inline('abs(1./(712.7584+j*(2.*pi.*X./(4*pi^2*3.241^2*1e10*c(1))-1./(2.*pi.*X*c(1))))','c','X');
opt=optimset('TolFun',1e-15);
c=lsqcurvefit(f,[6.657e-10], X, Y,opt);
A=c(1);
f是表示复数函数的取模 j是虚数单位
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

感谢参与,应助指数 +1
1.f少个)
2.查看lsqcurvefit的用法
clc
clear
X=[100000.0, 100676.7, 101353.4, 102030.1, 102706.81, 103383.51, 104060.21];
Y=[0.0003476636, 0.0003481631, 0.0003485057, 0.0003534742, 0.0003585374, 0.0003561982, 0.0003610355];
f=inline('abs(1./(712.7584+j*(2.*pi.*X./(4*pi^2*3.241^2*1e10*c(1))-1./(2.*pi.*X*c(1)))))','c','X');
opt=optimset('TolFun',1e-15);
c=lsqcurvefit(f,[6.657e-10], X, Y,[],[],opt);
A=c(1);
showmethemoney
2楼2013-07-05 17:36:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 菜头_ 的主题更新
信息提示
请填处理意见