24小时热门版块排行榜    

查看: 959  |  回复: 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的回帖

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的回帖
查看全部 3 个回答

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的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿吉林大学材料学硕321求调剂 +9 Ymlll 2026-03-18 12/600 2026-03-20 00:01 by 23Postgrad
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 生物学调剂招人!!! +3 山海天岚 2026-03-17 4/200 2026-03-19 21:34 by 怎么释怀
[考研] 能源材料化学课题组招收硕士研究生8-10名 +4 脱颖而出 2026-03-16 12/600 2026-03-19 16:17 by 脱颖而出
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 332求调剂 +3 ydfyh 2026-03-17 3/150 2026-03-19 10:14 by 功夫疯狂
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +7 步川酷紫123 2026-03-13 7/350 2026-03-18 17:12 by 尽舜尧1
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 一志愿西南交大,求调剂 +4 材化逐梦人 2026-03-18 4/200 2026-03-18 14:22 by 007_lilei
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 277调剂 +5 自由煎饼果子 2026-03-16 6/300 2026-03-17 19:26 by 李leezz
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 一志愿南京大学,080500材料科学与工程,调剂 +4 Jy? 2026-03-16 4/200 2026-03-17 11:02 by gaoqiong
[考研] 321求调剂 +5 大米饭! 2026-03-15 5/250 2026-03-16 16:33 by houyaoxu
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
信息提示
请填处理意见