| 查看: 1920 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
chaiqing123金虫 (小有名气)
|
[求助]
MATLAB用于三元回归参数
|
|
|
function文件为下 function rcal = fun552(x,x1,x2) for n=1:11%m为所取的数据点个数 R=8.315; global x1 x2 x3 T; tao12=x(1)./R/T(n);%x(1)=g12-g22 tao21=x(2)./R/T(n); %x(2)=g21-g11 tao13=x(3)./R/T(n); %x(3)=g13-g33 tao31=x(4)./R/T(n); %x(4)=g31-g11 tao23=x(5)./R/T(n); %x(5)=g23-g33 tao32=x(6)./R/T(n); %x(6)=g32-g22 tao=[tao12,tao13,tao21,tao23,tao31,tao32]; G=exp(-0.3*tao); A(n)=(tao(3)*G(3)*x2(n)+tao(5)*G(5)*x3(n))/(x1(n)+G(3)*x2(n)+G(5)*x3(n)); B(n)=-x1(n)*(x2(n)*tao(3)*G(3)+x3(n)*tao(5)*G(5))/((x1(n)+G(3)*x2(n)+G(5)*x3(n))^2); C(n)=x2(n)*G(1)*(tao(1)-(x1(n)*tao(1)*G(1)+x3(n)*tao(6)*G(6))/(G(1)*x1(n)+x2(n)+G(6)*x3(n)))/(G(1)*x1(n)+x2(n)+G(6)*x3(n)); D(n)=x3(n)*G(2)*(tao(2)-(x1(n)*tao(2)*G(2)+x2(n)*tao(4)*G(4))/(G(2)*x1(n)+G(4)*x2(n)+x3(n)))/(G(2)*x3(n)+G(4)*x2(n)+x3(n)); r1cal(n)=exp(A(n)+B(n)+C(n)+D(n)); A1(n)=(tao(1)*G(1)*x1(n)+tao(6)*G(6)*x3(n))/(G(1)*x1(n)+G(6)*x3(n)+x2(n)); B1(n)=-x2(n)*(x1(n)*tao(1)*G(1)+x3(n)*tao(6)*G(6))/((G(1)*x1(n)+x2(n)+G(6)*x3(n))^2); C1(n)=x1(n)*G(3)/(x1(n)+G(3)*x2(n)+G(5)*x3(n))*(tao(3)-(x2(n)*tao(3)*G(3)+x3(n)*tao(5)*G(5))/(x1(n)+G(3)*x2(n)+G(5)*x3(n))); D1(n)=x3(n)*G(4)/(G(2)*x1+G(4)*x2(n)+x3(n))*(tao(4)-(x1(n)*tao(2)*G(2)+x2(n)*tao(4)*G(4))/(G(2)*x1(n)+G(4)*x2(n)+x3(n))); r2cal(n)=exp(A1(n)+B1(n)+C1(n)+D1(n)); rcal=[r1cal r2cal]; end script文件如下 xx1=[0.05 0.105 0.153 0.248 0.347 0.446 0.474 0.535 0.624 0.740 0.870];%不含离子液体的x1 xx2=1-xx1; x1=xx1/1.05;%离子液体浓度为5% x2=xx2/1.05; x3=1-x1-x2; global x1 x2 x3 T; y1=[0.397 0.483 0.521 0.601 0.635 0.676 0.685 0.713 0.745 0.825 0.897]; y2=1-y1; T=[364 360.2 358.5 356.3 355.3 354.5 354.3 354.1 353.9 353.7 353.6]; t=T-273.15; a=[7.23255 1750.286 235;7.16879 1552.601 222.419]; p1s=10.^(a(1,1)-a(1,2)./(t+a(1,3))); p2s=10.^(a(2,1)-a(2,2)./(t+a(2,3))); %A=7.23255; B=1750.286;C=235; %A1=7.16879;B1=1552.601;C1=222.419; %p1s=10.^(A-B./(t+C)) %p2s=10.^(A1-B1./(t+C1)) p3s=0; %a为aotonie方程的三个系数 P=101.3%总压 r1exp=y1.*P./(x1.*p1s); r2exp=y2.*P./(x2.*p2s); rexp=[r1exp r2exp] x0=[1,1,1,1,1,1]; xx=[xx1 xx1]; x=nlinfit(xx,rexp,'fun552',x0)  最小二乘法 这个有问题 但不知道怎么改.... |
回复此楼» 猜你喜欢
青基已中
已经有1人回复
-
散金币祈福
已经有186人回复
-
物理学I论文润色/翻译怎么收费?
已经有75人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有109人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有64人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- matlab反求未知输入参数
已经有4人回复
- 求回归方程参数估计的最优化算法matlab代码
已经有9人回复
- matlab非线性参数拟合问题
已经有7人回复
- 如何在matlab中绘制含有三个参数方程的三维图形
已经有6人回复
- 【求助】matlab 求解微分方程中的未知参数
已经有20人回复
- 有人学过<神经网络>和《模糊数学》还有《matlab》这三门课吗？名字都好虚幻~
已经有20人回复
- 【求助】三元二次方程的回归拟合
已经有5人回复
- 三元方程用matlab画图，求指令
已经有1人回复
chaiqing123
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 679.7
- 红花: 1
- 帖子: 71
- 在线: 65.2小时
- 虫号: 1704607
- 注册: 2012-03-20
- 性别: GG
- 专业: 化工热力学和基础数据
4楼2012-04-01 14:48:51
2楼2012-04-01 10:52:29
chaiqing123
金虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 679.7
- 红花: 1
- 帖子: 71
- 在线: 65.2小时
- 虫号: 1704607
- 注册: 2012-03-20
- 性别: GG
- 专业: 化工热力学和基础数据
3楼2012-04-01 14:46:38
10