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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab高手帮忙做个自定义复杂函数拟合，含有复数，
261/3返回列表123下一页
查看: 3070  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 8rmuugja8q 2026-02-22 8/400 2026-02-23 12:22 by alian_214
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低 +10 瞬息宇宙 2026-02-20 13/650 2026-02-23 11:23 by holypower
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 8/400 2026-02-23 09:35 by w4l55oybr1
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 8/400 2026-02-23 09:29 by w4l55oybr1
[教师之家] 为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？ +5 欢乐颂叶蓁 2026-02-21 8/400 2026-02-23 09:19 by 欢乐颂叶蓁
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 5/250 2026-02-23 08:04 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 6/300 2026-02-23 07:59 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 9/450 2026-02-23 07:49 by w4l55oybr1
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[基金申请] 面上可以超过30页吧？ +4 阿拉贡aragon 2026-02-22 4/200 2026-02-22 21:22 by 山西悬空寺空悬�
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-21 4/200 2026-02-22 16:27 by khieu8v8m0
[基金申请] “人文社科而论，许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的？ +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错，点赞！ +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
信息提示
关闭
请填处理意见
关闭