| 查看: 2281 | 回复: 4 | |||||
刘小嵩321银虫 (著名写手)
|
[求助]
matlab出现问题 已有2人参与
|
|
1 function deconvolution 2 data=load('0.10-645.dat'); 3 xdata=data(:,1); 4 ydata=data(:,2); 5 x0=[0.1]; 6 %x(1)为脉宽，x(2)为时间零点，x(3)为A1，x(4)为t1，x(5)为A2，x(6)为t2 7 %----------------------------------------------------------- 8 options=optimset('TolFun',1e-4,'TolX',1e-4,'MaxIter',1e9,'MaxFunEvals',1e9); 9 [x,resnorm,residual,exitflag,output]=lsqcurvefit(@fun,x0,xdata,ydata,options); 10 disp('the fitted parameter is  delta mu A1 tau1 A2 tau2)');disp(x);11 disp('the resnorm is:');disp(resnorm); 12 disp('the exit flag is:');disp(exitflag); 13 disp(output); 14 %--------------------------------------------------------------------------fit 15 Fp=fun(x,xdata); 16 irf=-1*2*sqrt(log(2))/0.06/sqrt(pi)*exp(-log(2)*(2*(xdata-x(2))/x(1)).^2)/20; 17 subplot(2,1,1); plot(data(:,1),Fp,'b-',xdata,ydata,'ro',xdata,irf,'g'); 18 subplot(2,1,2); plot(data(:,1),residual,'g*'); 19 %------------------------------------------------------------------------- 20 result=[xdata ydata Fp residual irf]; 21 save -ASCII -DOUBLE result.dat result; 22 y=[x resnorm]; 23 save -ASCII -DOUBLE parameters.dat y; 24 function F=fun(x,xdata) 25 delta1=x(1)/(sqrt(log(4))); 26 F1=x(3).*exp(delta1^2/x(4)^2/2).*exp(-(xdata-x(2))/x(4)).*(1+erf(((xdata-x(2))/delta1-delta1/x(4))/sqrt(2))); 27 F2=x(5).*(1+erf((xdata-x(7))/sqrt(2)/delta1)-exp(delta1^2/x(6)^2/2).*exp(-(xdata-x(7))/x(6)).*(1+erf(((xdata-x(7))/delta1-delta1/x(6))/sqrt(2)))); 28 F=F1+F2; 执行命令时出现： 错误使用 lsqcurvefit (line 181) LSQCURVEFIT requires the following inputs to be of data type double: 'LB'. 出错 Untitled2 (line 9) [x,resnorm,residual,exitflag,output]=lsqcurvefit(@fun,x0,xdata,ydata,options); 应该如何修改？ |
回复此楼» 猜你喜欢
如何从铁电相到顺电相。
已经有1人回复
-
各位大佬，求一份最新Wien2k版本作为学习用途
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有265人回复
-
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后
已经有42人回复
-
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向
已经有24人回复
-
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 运行matlab出现的问题
已经有0人回复
- matlab 版本问题
已经有2人回复
- 多脉波整流matlab仿真出现问题，求助
已经有2人回复
- matlab问题
已经有1人回复
- Matlab2010运行出现问题
已经有0人回复
- AMD64CPU系统装MATLAB2013出现问题，求助
已经有1人回复
- 你在matlab中出现这样的问题吗？
已经有12人回复
- 以下matlab代码出现问题
已经有10人回复
- matlab一运行就出现的问题
已经有7人回复
- 运行matlab出现错误
已经有4人回复
- matlab问题
已经有4人回复
- matlab中问题请教
已经有6人回复
- 用matlab绘曲线，但总是出现问题
已经有5人回复
- matlab 程序问题
已经有2人回复
- matlab安装问题
已经有2人回复
- matlab问题求助
已经有16人回复
- matlab仿真时出现的问题
已经有2人回复
- matlab问题
已经有1人回复
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 计算强帖: 12
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
X = lsqcurvefit(FUN,X0,XDATA,YDATA,LB,UB,OPTIONS) minimizes with the default parameters replaced by values in the structure OPTIONS, an argument created with the OPTIMSET function. See OPTIMSET for details. Use the Jacobian option to specify that FUN also returns a second output argument J that is the Jacobian matrix at the point X. If FUN returns a vector F of m components when X has length n, then J is an m-by-n matrix where J(i,j) is the partial derivative of F(i) with respect to x(j). (Note that the Jacobian J is the transpose of the gradient of F.) 改为[x,resnorm,residual,exitflag,output]=lsqcurvefit(@fun,x0,xdata,ydata,[],[],options); |
2楼2016-03-25 14:19:03
刘小嵩321
银虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 4937.6
- 散金: 419
- 红花: 3
- 帖子: 1130
- 在线: 59.7小时
- 虫号: 2468827
- 注册: 2013-05-17
- 性别: MM
- 专业: 光学
|
还是报错， 出错 Untitled2 (line 9) [x,resnorm,residual,exitflag,output]=lsqcurvefit(@fun,x0,xdata,ydata,[],[],options); 原因: Failure in initial user-supplied objective function evaluation. LSQCURVEFIT cannot continue. |
» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : QQ截图20160325165443.png
2016-03-25 16:55:21, 48.95 K
3楼2016-03-25 16:55:43
4楼2016-03-25 18:29:00
刘小嵩321
银虫 (著名写手)
- 应助: 13 (小学生)
- 金币: 4937.6
- 散金: 419
- 红花: 3
- 帖子: 1130
- 在线: 59.7小时
- 虫号: 2468827
- 注册: 2013-05-17
- 性别: MM
- 专业: 光学
5楼2016-03-28 07:08:06
100