[問題]求解10個係數的問題~求救...已回收
log(Pv)已知 log(Pv)=1.7320,2.4524,3.0197,3.5522,
4.0534,4.5257,4.9719,2.5268,3.2170,3.1449
t已知 t=298.15,310,320,330,340,350,360
w已知 w=0.1,0.2,0.3
求a0~a4,b0~b4的係數?
log(Pv)=((a0+a1*w+a2*w^2+a3*w^3+a4*w^4)+(b0+b1*w+b2*w^2+b3*w^3+b4*w^4)/t)
小弟因為要使用上式來代替三條線的數據
只是跑出來 不知是哪裡出錯了....懇請高手救命~謝謝!
明明應該有10個係數不知為何剩5個而已...
=============================================================================
clear
a=[1, 0.1, 0.1^2, 0.1^3, 0.1^4;
1, 0.1, 0.1^2, 0.1^3, 0.1^4;
1, 0.1, 0.1^2, 0.1^3, 0.1^4;
1, 0.1, 0.1^2, 0.1^3, 0.1^4;
1, 0.1, 0.1^2, 0.1^3, 0.1^4;
1, 0.1, 0.1^2, 0.1^3, 0.1^4;
1, 0.1, 0.1^2, 0.1^3, 0.1^4;
1, 0.2, 0.2^2, 0.2^3, 0.2^4;
1, 0.2, 0.2^2, 0.2^3, 0.2^4;
1, 0.3, 0.3^2, 0.3^3, 0.3^4];
B0=[1, 0.1, 0.1^2, 0.1^3, 0.1^4]/298.15;
B1=[1, 0.1, 0.1^2, 0.1^3, 0.1^4]/310;
B2=[1, 0.1, 0.1^2, 0.1^3, 0.1^4]/320;
B3=[1, 0.1, 0.1^2, 0.1^3, 0.1^4]/330;
B4=[1, 0.1, 0.1^2, 0.1^3, 0.1^4]/340;
B5=[1, 0.1, 0.1^2, 0.1^3, 0.1^4]/350;
B6=[1, 0.1, 0.1^2, 0.1^3, 0.1^4]/360;
B7=[1, 0.2, 0.2^2, 0.2^3, 0.2^4]/298.15;
B8=[1, 0.2, 0.2^2, 0.2^3, 0.2^4]/310;
B9=[1, 0.3, 0.3^2, 0.3^3, 0.3^4]/298.15;
B=[B0; B1; B2; B3; B4; B5; B6; B7; B8; B9];
C=[1.7320;
2.4524;
3.0197;
3.5522;
4.0534;
4.5257;
4.9719;
2.5268;
3.2170;
3.1449;
];
solution = (a+B)\C
Warning: Rank deficient, rank = 3, tol = 7.0436e-015.
solution =
4.9322
-19.0606
43.5598
0
0
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 125.230.41.109
討論串 (同標題文章)
以下文章回應了本文：
完整討論串 (本文為第 1 之 2 篇)：
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
