应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab高手帮忙做个自定义复杂函数拟合，含有复数，
北京石油化工学院2026年研究生招生接收调剂公告
261/3返回列表123下一页
查看: 3189  |  回复: 25
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

delarzhang

木虫 (著名写手)

无

[求助] matlab高手帮忙做个自定义复杂函数拟合，含有复数，

需要对一个带复数的复杂函数进行拟合，并求出几个参量，本人是matlab菜鸟，完成后请一定给出源代码，具体函数请查看附件说明，多谢！
其它软件也可以，我只有matlab和origin可以用，我要处理的数据比较多，那就要麻烦您多算几次了，站内联系。filename]matlab拟合.doc[/filename]

[ Last edited by delarzhang on 2012-8-17 at 08:13 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
    本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com
  • 附件 1 : matlab拟合.doc
    • 2012-08-17 08:05:22, 99.5 K

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
天道酬勤
1楼 2012-08-17 08:06:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
delarzhang: 金币+50, ★有帮助 2012-09-04 07:21:58
delarzhang: 金币+50 2012-10-03 16:12:55
你自己试试吧，能拟合，但是结果不是很理想
CODE:
function delar
data=xlsread('delar.xlsx');
xdata=data(:,1);  %Im
ydata=data(:,2);  %
x0=[100,250,6,0.5];    %wp2,gama,T,S
options= optimset('MaxFunEvals',5000,'Maxiter',1000);
lb=[10,20,0,0];
%ub=[10000,10000,100,1000];
ub=[];
[x,resnorm]=lsqcurvefit(@myfun,x0,xdata,ydata,lb,ub,options)
y=myfun(x,xdata);
plot(xdata,ydata,'g',xdata,y,'r')
end

function F=myfun(x,xdata)  %x=[]
wl=734;
wt=533;
e8=5.5;
C=0.48;
eq1=(wl^2-wt^2)./(wt^2-xdata.^2-1i*x(3)*xdata);
eq2=x(1)./(xdata.^2+1i*xdata*x(2));
epslong=e8*(1+eq1-eq2);

delta=x(1)*x(2)*((wt^2-xdata.^2).^2+(xdata*x(3)).^2)  +  xdata.^2*x(3)*(wl^2-wt^2).*(xdata.^2+x(2)^2);
eq3=2*C*wt^2* (x(1)*x(2)*(wt^2-xdata.^2)   -  xdata.^2*x(3).*(xdata.^2+x(2)^2-x(1)))   ./delta;
eq4=C^2*(wt^4./delta)  .*(x(1)*(x(2)*(wl^2-wt^2)+x(3)*(x(1)-2*xdata.^2))  +xdata.^2*x(3).*(xdata.^2+x(2)^2))  ./(wl^2-wt^2);
A=1+eq3+eq4;
F=x(4)*A.*imag(-1./epslong);
end
赞一下(2人) 回复此楼
showmethemoney
6楼2012-08-20 10:59:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

love_ff

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by csgt0 at 2012-08-20 10:59:47
你自己试试吧，能拟合，但是结果不是很理想

function delar
data=xlsread('delar.xlsx');
xdata=data(:,1);  %Im
ydata=data(:,2);  %
x0=;    %wp2,gama,T,S
options= optimset('MaxFunEvals',5000,'Maxi ...

请问我用lsqcurvefit函数的时候我设定了上下线，但是好像没用。我把代码贴出来能帮我看看么？
而且，我拟合得到的系数都是复数，可否得到实数。麻烦了！
x0=[2,2,2,2,2];
xdata=w0';
ydata=xdl2;
lb=[1;1;1;1;1];
lu=[5;5;5;5;5];
[x,res]=lsqcurvefit(F,x0,xdata,ydata,lb,lu);
赞一下(1人) 回复此楼
24楼2014-01-02 10:26:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

感谢参与，应助指数 +1
复数拟合还没见过比1stOpt更强更好用的，建议自己试试。
赞一下 回复此楼
2楼2012-08-17 09:26:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

delarzhang

木虫 (著名写手)

无
引用回帖:
2楼: Originally posted by dingd at 2012-08-17 09:26:06
复数拟合还没见过比1stOpt更强更好用的，建议自己试试。

如果我自己有，就不用浪费时间在这里求助
赞一下 回复此楼
天道酬勤
3楼2012-08-17 10:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图
I是多少，m是多少
回复此楼
showmethemoney
4楼2012-08-17 17:49:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

delarzhang

木虫 (著名写手)

无
引用回帖:
4楼: Originally posted by csgt0 at 2012-08-17 17:49:20
I是多少，m是多少

Im是复数的虚部，Re复数的实部
赞一下 回复此楼
天道酬勤
5楼2012-08-17 20:13:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

delarzhang

木虫 (著名写手)

无
引用回帖:
6楼: Originally posted by csgt0 at 2012-08-20 10:59:47
你自己试试吧，能拟合，但是结果不是很理想

function delar
data=xlsread('delar.xlsx');
xdata=data(:,1);  %Im
ydata=data(:,2);  %
x0=;    %wp2,gama,T,S
options= optimset('MaxFunEvals',5000,'Maxi ...

谢谢，我安装的是matlab2012a，运行时显示resnorm是unused，不知道该怎么办？
赞一下 回复此楼
天道酬勤
7楼2012-08-21 07:43:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

delarzhang

木虫 (著名写手)

无
引用回帖:
6楼: Originally posted by csgt0 at 2012-08-20 10:59:47
你自己试试吧，能拟合，但是结果不是很理想

function delar
data=xlsread('delar.xlsx');
xdata=data(:,1);  %Im
ydata=data(:,2);  %
x0=;    %wp2,gama,T,S
options= optimset('MaxFunEvals',5000,'Maxi ...

运行过程中显示信息如下：
Error using xlsread (line 129)
XLSREAD unable to open file 'delar.xlsx'.
File 'D:\Program Files\MATLAB\R2012a\bin\raman\delar.xlsx' not found.

Error in delar (line 2)
data=xlsread('delar.xlsx');
赞一下 回复此楼
天道酬勤
8楼2012-08-21 07:49:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

引用回帖:
8楼: Originally posted by delarzhang at 2012-08-21 07:49:39
运行过程中显示信息如下：
Error using xlsread (line 129)
XLSREAD unable to open file 'delar.xlsx'.
File 'D:\Program Files\MATLAB\R2012a\bin\raman\delar.xlsx' not found.

Error in delar (line 2) ...

把你的数据在excel里保存为两列，命名为delar.xlsx，跟程序放到一起。
resnorm那个不用的话就删掉
赞一下 回复此楼
showmethemoney
9楼2012-08-21 09:02:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

delarzhang

木虫 (著名写手)

无
引用回帖:
9楼: Originally posted by csgt0 at 2012-08-21 09:02:43
把你的数据在excel里保存为两列，命名为delar.xlsx，跟程序放到一起。
resnorm那个不用的话就删掉...

请帮满看看，还是报错：
Error using snls (line 48)
Objective function is returning undefined values at initial point. lsqcurvefit
cannot continue.

Error in lsqncommon (line 175)
    [xC,FVAL,LAMBDA,JACOB,EXITFLAG,OUTPUT,msgData]=...

Error in lsqcurvefit (line 252)
[xCurrent,Resnorm,FVAL,EXITFLAG,OUTPUT,LAMBDA,JACOB] = ...

Error in delar (line 10)
[x,resnorm]=lsqcurvefit(@myfun,x0,xdata,ydata,lb,ub,options);
赞一下 回复此楼
天道酬勤
10楼2012-08-21 10:09:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 delarzhang 的主题更新
261/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323求调剂 +3 林zlu 2026-04-07 4/200 2026-04-07 23:21 by lbsjt
[考研] 计算机11408 287 求调剂 +3 LiLe5 2026-04-07 3/150 2026-04-07 23:15 by shanqishi
[考研] 材料工程322分 +14 哈哈哈吼吼吼哈 2026-04-01 14/700 2026-04-07 16:14 by 静花儿
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +12 哇呼哼呼哼 2026-04-01 13/650 2026-04-07 10:02 by zhen～
[考研] 求调剂 +4 wos666 2026-04-03 5/250 2026-04-06 15:22 by wos666
[考研] 一志愿安徽某211 0703化学总分339求调剂 +7 晚风不晚 2026-04-04 7/350 2026-04-06 14:06 by houyaoxu
[考研] 一志愿武汉理工大学080200机械工程308分，求调剂 +4 终不似从前 2026-04-05 4/200 2026-04-06 11:46 by 考研学校招点人
[考研] 301求调剂 +6 细胞相关蛋白 2026-04-02 10/500 2026-04-06 08:34 by jp9609
[考研] 一志愿C9的化学工程（085602） 340分，感觉校内调剂无望，求调剂 +12 万事宜臻 2026-04-04 12/600 2026-04-06 07:46 by 无际的草原
[考研] 求调剂求调剂 +8 121. 2026-04-02 8/400 2026-04-05 20:15 by lys0704
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 工科08-机械专硕-求调剂 +3 雷欧飞踢 2026-04-02 3/150 2026-04-05 18:49 by 蓝云思雨
[考研] 材料与化工306分找调剂 +23 沧海轻舟e 2026-04-02 27/1350 2026-04-04 21:52 by laoshidan
[考研] 085400电子信息319求调剂（接受跨专业调剂） +5 星星不眨眼喽 2026-04-03 6/300 2026-04-04 21:50 by hemengdong
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
[考研] 350一志愿北京航空航天大学08500材料科学与工程求调剂 +5 kjnasfss 2026-04-03 5/250 2026-04-03 22:29 by 无际的草原
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 工科 267求调剂 +5 wanwan00 2026-04-02 7/350 2026-04-03 14:14 by zhangdingwa
[考研] 复试调剂 +3 bvzz 2026-04-01 3/150 2026-04-03 09:47 by 蓝云思雨
[考研] 安全工程 285 求调剂 +3 Xinyu56 2026-04-01 4/200 2026-04-01 21:50 by 静静静静静静静�
信息提示
关闭
请填处理意见
关闭