Re: [問題]求解10個係數的問題~求救...已回收

看板MATLAB作者moussorgsky (法國號有氣質)時間16年前 (2009/09/24 15:40)推噓0(0推 0噓 0→)

留言0則, 0人參與討論串2/2 (看更多)

※ 引述《maxdeity (deity)》之銘言： : 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 ^^^^^^^ 這邊寫相反了 C\(a+B)才對 : Warning: Rank deficient, rank = 3, tol = 7.0436e-015. : solution = : 4.9322 : -19.0606 : 43.5598 : 0 : 0 -- 年輕是什麼？是風，是雲彩，也是天空，是一種心情，閃爍在生命的每一個轉折裡。 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.187.195