| 查看: 2205 | 回复: 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); 应该如何修改？ |
回复此楼» 猜你喜欢
光学工程学硕调剂信息
已经有32人回复
-
欢迎加入课题组
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有208人回复
-
散金币，求好运，祝面上顺利！
已经有79人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有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
40