24小时热门版块排行榜    

查看: 955  |  回复: 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的回帖
相关版块跳转 我要订阅楼主 不惑可乐 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博26年 +3 八6八68 2026-03-19 3/150 2026-03-19 19:43 by nxgogo
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +3 葵梓卫队 2026-03-18 5/250 2026-03-19 19:35 by 给你你注意休息
[考博] 东华理工大学化材专业26届硕士博士申请 +8 zlingli 2026-03-13 8/400 2026-03-19 16:32 by 轻松不少随
[考研] 266求调剂 +5 阳阳哇塞 2026-03-14 10/500 2026-03-19 15:08 by 阳阳哇塞
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +10 Liwangman 2026-03-15 10/500 2026-03-19 10:25 by 无际的草原
[考研] 330求调剂 +3 小材化本科 2026-03-18 3/150 2026-03-18 21:55 by 无懈可击111
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-17 7/350 2026-03-18 20:41 by Wangjingyue
[考研] 材料专业求调剂 +5 hanamiko 2026-03-18 5/250 2026-03-18 20:19 by 星空星月
[考研] 344求调剂 +6 knight344 2026-03-16 7/350 2026-03-18 20:13 by walc
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考研] 材料与化工专硕调剂 +5 heming3743 2026-03-16 5/250 2026-03-17 14:03 by 勇敢太监王公公
[考研] 考研调剂 +3 淇ya_~ 2026-03-17 5/250 2026-03-17 09:25 by Winj1e
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 289求调剂 +4 这么名字咋样 2026-03-14 6/300 2026-03-14 18:58 by userper
[考研] 295求调剂 +3 小匕仔汁 2026-03-12 3/150 2026-03-13 15:17 by vgtyfty
信息提示
请填处理意见