24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 729  |  回复: 5
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

qcz2002

新虫 (小有名气)

[交流] 【求助】对于特定的函数用matlab做最小二乘法的曲线拟合

函数为n^2=x(1)+x(2)/[r^2-x(3)]-x(4)/[x(5)-r^2]
其中n 和r为已知的两组数据,我想得到其最小二乘法的拟合曲线以及5个参数x的值
我是这样弄得 结果老是无法运行
m文件
function f^2=f01(x,Data);
Data=[4358 2.42207;4861 2.40374;5461 2.40126;5893 2.39352;6563 2.37776;6678 2.37684];
r=Data(:,1);
n=Data(:,2);
z^2=x(1)+x(2)/[r^2-x(3)]-x(4)/[x(5)-r^2];
f^2=z^2-n^2;
输入命令
Data=[4358 2.42207;4861 2.40374;5461 2.40126;5893 2.39352;6563 2.37776;6678 2.37684];
r=Data(:,1);
n=Data(:,2);
x0=[1 0]';
x=lsqnonlin('f01',x0); 到这一步就出问题了,哪位好心的大侠帮帮我吧 应该怎么建m文件 程序错哪了 我刚刚接触matlab,我看了书上的例子,不过那个例子是一次的

[ Last edited by kuhailangyu on 2009-11-18 at 20:53 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qcz2002

新虫 (小有名气)

还是不行呢 可不可以帮我设计然后运行出来呀
3楼2009-11-19 15:55:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

zjliu1265

至尊木虫 (正式写手)

★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
nono2009(金币+3,VIP+0):谢谢应助!可能还需要进一步指导。 11-19 20:21
我来援助一下吧!
x0应该是含5个元素的向量,而你只写了2个。
我随意加了3个元素,不一定对,完全是为了
程序能运行。

另外:
1. 符号"^"是运算符,不要用于变量名字中;
2. 符号"^"和"/"在计算表达式中如果不是方矩阵,它们一般
需分别改写为".^"和"./"表示元素之间的乘除运算。

改后的程序是:
m文件
function f2=f01(x,Data);
Data=[4358 2.42207;4861 2.40374;5461 2.40126;5893 2.39352;6563 2.37776;6678 2.37684];
r=Data(:,1);
n=Data(:,2);
z2=x(1)+x(2)./[r.^2-x(3)]-x(4)./[x(5)-r.^2];
f2=z2-n.^2;
输入命令
Data=[4358 2.42207;4861 2.40374;5461 2.40126;5893 2.39352;6563 2.37776;6678 2.37684];
r=Data(:,1);
n=Data(:,2);
x0=[1 0 1 1 1]';
x=lsqnonlin('f01',x0);
2楼2009-11-19 09:23:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qcz2002

新虫 (小有名气)

Data=[435.8 1.5563;486.1 1.5504;546.1 1.5496;589.3 1.5471;656.3 1.542;667.8 1.5417]; 函数为n^2=x(1)+x(2)/[r^2-x(3)]-x(4)/[x(5)-r^2]
其中n 和r为已知的两组数据,x为5个参数值

谢谢帮帮我吧

[ Last edited by qcz2002 on 2009-11-19 at 16:01 ]
4楼2009-11-19 15:59:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qcz2002

新虫 (小有名气)

nono2009(金币+0,VIP+0):建议通过PM或“引用回复该帖”,以便别人及时收到你的message. 11-19 20:20
运行情况>> Data=[435.8 1.5563;486.1 1.5504;546.1 1.5496;589.3 1.5471;656.3 1.5420;667.8 1.5417];
>> r=Data(:,1);
>> n=Data(:,2);
>> x0=[1 0 1 1 1]';
>> x=lsqnonlin('f01',x0);
Optimization terminated: first-order optimality less than OPTIONS.TolFun,
and no negative/zero curvature detected in trust region model.
>> n2=x(1)+x(2)./[r.^2-x(3)]-x(4)./[x(5)-r.^2];
>> plot(r,n,'bol-');
??? Error using ==> plot
Error in color/linetype argument.
5楼2009-11-19 16:42:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600材料与化工调剂 +23 kikiki7 2026-03-30 23/1150 2026-04-03 23:43 by hemengdong
[考研] 化工求调剂 +11 荔香芝士椰奶 2026-04-03 11/550 2026-04-03 22:06 by 啵啵啵0119
[考研] 求调剂 +3 wos666 2026-04-03 3/150 2026-04-03 21:36 by lbsjt
[考研] 294求调剂 +6 Grey_Ey 2026-04-03 6/300 2026-04-03 20:46 by 欣喜777
[考研] 材料调剂 +15 一样YWY 2026-04-02 15/750 2026-04-03 15:39 by 晴空210210
[考研] 求调剂 +9 akdhjs 2026-03-31 11/550 2026-04-03 13:32 by akdhjs
[考研] 279求调剂 +6 qazplm0852 2026-04-02 6/300 2026-04-03 10:03 by 蓝云思雨
[考研] 调剂 +7 祉岷. 2026-04-02 7/350 2026-04-03 09:11 by 花呗还欠600
[考研] 一志愿大工学硕,求调剂 +4 yub0811 2026-04-02 4/200 2026-04-02 21:36 by 百灵童888
[考研] 求调剂 +7 Aniyaio 2026-04-02 7/350 2026-04-02 16:42 by zzsw+
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +10 哇呼哼呼哼 2026-04-01 11/550 2026-04-01 21:48 by chyhaha
[考研] 332求调剂 +8 Lyy930824@ 2026-03-29 8/400 2026-04-01 18:40 by 千钧澄玉宇啊
[考研] 材料专业调剂 +5 啦啦啦哭 2026-03-31 6/300 2026-04-01 16:48 by JourneyLucky
[考研] 一志愿武理材料工程302调剂环化或化工 +15 Doleres 2026-03-31 16/800 2026-04-01 09:49 by lfj11
[考研] 254材料与化工求调剂 +3 翰冬林楠 2026-03-30 4/200 2026-03-31 17:53 by yishunmin
[考研] 11408总分309,一志愿东南大学求调剂,不挑专业 +5 天赋带到THU 2026-03-29 6/300 2026-03-30 20:49 by dick_runner
[考研] 求佛 +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
信息提示
请填处理意见