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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » matlab高手帮忙做个自定义复杂函数拟合，含有复数，
北京石油化工学院2026年研究生招生接收调剂公告
261/3返回列表123下一页
查看: 3185  |  回复: 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下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料307分求大佬组收留 +6 Hll胡 2026-04-07 6/300 2026-04-07 23:06 by jp9609
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 315求调剂 +17 小羊小羊_ 2026-04-02 18/900 2026-04-07 22:01 by lijunpoly
[考研] 环境专硕调剂 +12 会说话的肘子 2026-04-06 12/600 2026-04-07 17:23 by 可口可乐不加冰�
[考研] 081700化学工程与技术 一志愿中海洋 323 求调剂学校 +19 披星河 2026-04-03 19/950 2026-04-07 15:14 by 尽舜尧1
[考研] 285求调剂 +16 哦呦呼o 2026-04-04 16/800 2026-04-07 12:31 by 1018329917
[考研] 生物学调剂 可调剂到生物与医药 +3 李政莹 2026-04-06 3/150 2026-04-06 19:02 by macy2011
[考研] 269电子信息求调剂，可转专业 +5 独酌wl 2026-04-06 5/250 2026-04-06 17:23 by 土木硕士招生
[考研] 生物学求调剂 +5 15064154688 2026-04-03 5/250 2026-04-06 11:56 by lijunpoly
[考研] 372分，材料与化工，一志愿湖南大学，求调剂 +3 蓝笺片 2026-04-01 3/150 2026-04-06 09:04 by 无际的草原
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 280求调剂 +21 咕噜晓晓 2026-04-02 22/1100 2026-04-04 11:12 by 猪会飞
[考研] 336求调剂 +8 kiyy 2026-04-01 8/400 2026-04-03 19:41 by lijunpoly
[考研] 一志愿南昌大学324求调剂 +13 hanamiko 2026-04-01 13/650 2026-04-03 18:30 by ls刘帅
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +9 哇呼哼呼哼 2026-04-02 9/450 2026-04-03 12:05 by 1753564080
[考研] 建环，能源，土木老师路过看一看！！！ +5 嘿嘿uu 2026-04-01 5/250 2026-04-03 11:47 by znian
[考研] 求调剂 +3 心想事成可 2026-04-03 3/150 2026-04-03 11:22 by wangjy2002
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考博] 申博求助 +3 Reee1Llll 2026-04-01 3/150 2026-04-02 22:29 by 这是一个无聊的�
[考研] 求调剂 +7 Aniyaio 2026-04-02 7/350 2026-04-02 16:42 by zzsw+
信息提示
关闭
请填处理意见
关闭