24小时热门版块排行榜    

查看: 956  |  回复: 2

yiran909

至尊木虫 (文坛精英)

孤独的跑者

[交流] matlab最小二乘曲线约束拟合编程问题 已有1人参与

各位老师,给指导下,我编写的matlab最小二乘曲线有约束拟合
其中一个约束条件a(1)+a(2)+a(3)=1 这个约束条件不知道该怎么设置
帮看看程序中的问题,编译通不过,谢谢
function f=polyfun(a,x)%x是独立变量应变,a是待定系数
f=a(1)*u(1)+a(2)*v(1)+a(3)*w(1)+a(4)*x^a(5)+a(6);%这个函数定义的是否正确呢?
end

clear;
clc;
%多晶拟合单晶体积分数
x=[0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 0.12 0.13 0.14];% x data
y=[58.63 68.47 79.80 92.49 102.40 113.67 124.95 134.85 146.86 156.03 163.77 174.41 182.15 ];% y data
u=[22.59 31.76 40.24 48.72 56.46 63.53 71.32 78.38 84.02 89.65 95.34 99.55 104.51];%[411]
v=[51.56 69.22 84.02 99.55 110.83 122.84 131.32 142.59 147.54 153.91 160.24 164.50 168.04];%[100]
w=[80.48 115.09 141.22 168.22 195.58 218.86 242.83 264.75 281.71 297.19 312.05 322.64 331.80];%[111]
a0=[0.2,0.5,0.3,50,0.3,80]; % initial values
lb=[0 0 0 0 0 0 ];
ub=[1 1 1 inf inf inf];
%a(1)+a(2)+a(3)=1 这个约束条件不知道该怎么设置
a=lsqcurvefit(@polyfun,a0,x,y,lb,ub)%help lsqcurvefit查看最小二乘拟合函数的用法
f=polyfun(a,x);
plot(x,y,'*',x,f,'r-');
subs(a)%输出函数系数值
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dingd

铁杆木虫 (职业作家)


小木虫: 金币+0.5, 给个红包,谢谢回帖
约束拟合问题,1stOpt倒是很简单方便:
CODE:
Parameter a(6)=[0,];
Variable x,y,u,v,w;
Function y=a1*u+a2*v+a3*w+a4*x^a5+a6;
         a1+a2+a3=1;
Data;
x=[0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 0.12 0.13 0.14];
y=[58.63 68.47 79.80 92.49 102.40 113.67 124.95 134.85 146.86 156.03 163.77 174.41 182.15 ];
u=[22.59 31.76 40.24 48.72 56.46 63.53 71.32 78.38 84.02 89.65 95.34 99.55 104.51];
v=[51.56 69.22 84.02 99.55 110.83 122.84 131.32 142.59 147.54 153.91 160.24 164.50 168.04];
w=[80.48 115.09 141.22 168.22 195.58 218.86 242.83 264.75 281.71 297.19 312.05 322.64 331.80];

均方差(RMSE):0.860758921413859
残差平方和(SSE):9.63177697031615
相关系数(R): 0.999759586120821
相关系数之平方(R^2): 0.999519230040476
确定系数(DC): 0.999519226423277
卡方系数(Chi-Square): 0.0417908673148741
F统计(F-Statistic): 2911.97367696933
约束条件: a1+a2+a3-1 = -1.21055721E-010

参数        最佳估算
----------        -------------
a1        0.999999999878933
a2        1.1216895450583E-14
a3        8.18442542490143E-26
a4        662.707342823122
a5        1.34795089754783
a6        31.6077209708818
matlab最小二乘曲线约束拟合编程问题
c108.jpg

» 本帖已获得的红花(最新10朵)

2楼2014-07-09 23:12:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yiran909

至尊木虫 (文坛精英)

孤独的跑者

送红花一朵
引用回帖:
2楼: Originally posted by dingd at 2014-07-09 23:12:11
约束拟合问题,1stOpt倒是很简单方便:

Parameter a(6)=;
Variable x,y,u,v,w;
Function y=a1*u+a2*v+a3*w+a4*x^a5+a6;
         a1+a2+a3=1;
Data;
x=;
y=;
u=;
v=;
w=;

均方差(RMSE):0.8607589214 ...

谢谢,我学习下。
3楼2014-07-10 09:09:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 不惑可乐 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 296求调剂 +3 www_q 2026-03-18 6/300 2026-03-19 22:28 by zhq0425
[考研] 0856调剂,是学校就去 +6 sllhht 2026-03-19 7/350 2026-03-19 19:50 by 制度的
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +5 枫桥ZL 2026-03-18 7/350 2026-03-19 14:52 by 功夫疯狂
[考研] 304求调剂 +3 曼殊2266 2026-03-18 3/150 2026-03-19 14:42 by peike
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-19 3/150 2026-03-19 13:36 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 一志愿天大材料与化工(085600)总分338 +5 蔡大美女 2026-03-13 5/250 2026-03-19 10:44 by 是小刘呀~
[考研] 材料专硕英一数二306 +5 z1z2z3879 2026-03-18 5/250 2026-03-19 07:43 by BruceLiu320
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 考研求调剂 +3 橘颂. 2026-03-17 4/200 2026-03-17 21:43 by 有只狸奴
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 0856专硕279求调剂 +5 加油加油!? 2026-03-15 5/250 2026-03-15 11:58 by 2020015
信息提示
请填处理意见