查看: 1129  |  回复: 3

perper

铜虫 (初入文坛)

[求助] matlab程序报错,请大神帮忙改正

M文件:
function f=myfunc(D,t)
syms n t;
x=1/n^2*exp((-1)*D*n^2*pi^2*t/0.0002^2);
s=symsum(x,n,1,100);
f=1-6/pi^2*s;                              
end
主程序:
  t=[180 420 780 1080 1380 1680 1920 2220 2520 2820 3120 3420 3720 4020 4380 4680 4980 5280];
c=[0.08762 0.20444 0.36545 0.48247 0.59 0.67602 0.75 0.80407 0.86 0.9 0.94 0.96831 0.98561 0.99 0.99968 0.999999 1 1];
D0=[0.001];
D=lsqcurvefit(@myfunc,D0,t,c);
f=myfunc(D,t);
错误:
Error using lsqcurvefit (line 247)
Function value and YDATA sizes are incommensurate.
是想通过非线性拟合求D。有t,c数据,方程是一个求和方程。
回复此楼

» 猜你喜欢

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

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

xiegangmai

版主 (职业作家)

我没头衔

优秀版主优秀版主优秀版主

perper: 回帖置顶 2014-05-26 18:49:49
两个问题:
(1)D0=[0.001];初始值应是二维的数组。
(2)myfunc不能用符号函数。
明德厚学、求是创新
2楼2014-05-26 17:05:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

perper

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by xiegangmai at 2014-05-26 17:05:13
两个问题:
(1)D0=;初始值应是二维的数组。
(2)myfunc不能用符号函数。

可是我只有一个未知数D啊,初值不是就一个?
BIGBOY
3楼2014-05-26 18:47:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

perper

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by xiegangmai at 2014-05-26 17:05:13
两个问题:
(1)D0=;初始值应是二维的数组。
(2)myfunc不能用符号函数。

哥们,很急啊 求帮忙
BIGBOY
4楼2014-05-26 18:50:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 perper 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 极限拉扯,如何抉择 +6 Isak0320 2026-10-12 8/400 2026-10-12 20:09 by caixiaoye1
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +4 wangyl_123 2026-10-10 8/400 2026-10-12 18:09 by wangyl_123
[教师之家] 课题组会帮你通过考核么? +3 好运来~9527 2026-10-10 6/300 2026-10-12 14:55 by 简单化xn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见