[自動轉寄] Re: [問題] 解聯立方程式

作者: qoodeva (發呆的花枝)
標題: Re: [問題] 解聯立方程式
時間: Sun Nov 13 17:55:51 2011

※ 引述《chungyuandye (養花種魚數月亮賞星星)》之銘言：
: ※ 引述《qoodeva (發呆的花枝)》之銘言：
: : Do[
: : zt=(aω^2 +c)+d e-d^2 q; (*這是方程式，a,b,c,e,q為常數,ω為未知數*)
: : sol=Solve[zt==0,d]; (求解d，但因為ω為未知數，會出現未精準值)
: : H=HH[[4;;6,4;;6]]; (*H是HH矩陣的子矩陣*)
: : HE=Eigenvalues[H];
: : HD=Det[H];
: : If[Sign[HD]=\=0, (*判定HD行列式不為零*)
: : If[
: : Sign[HD]=-Sign[zz],Print[ω] (*判定行列式大於零時，將對應的ω輸出*)
: : ];
: : ];
: : zz=HD;
: : ,{ω,1,12000,1}]
: zt=(aω^2 +c)+d e-d^2 q
: 你的aω^2中間少打一個等號或是空格
: 其他看起來沒有錯才對

喔喔，那只是打快了一點，實際上我是有空格的，不過謝謝你，
我今天重新試了好幾次，最後發現，我HH矩陣是寫在Do外面
，但是多打了個//MatrixForm，再變成矩陣顯示時似乎沒辦法執行Do裡面的code，
消去之後就有插值跑出來，但是我的研究還沒做完，
我也不確定這樣跑出來的數字正不正確，但是至少有東西在跑了，謝謝你黑。

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.116.78.76