我已衍生彈簧質量系統的組合的傳遞函數在s
域,它的形式爲計算傳遞函數係數
G = (as*s + bs +c)/(ps + q)
我有與位移涉及作爲輸入和力作爲測量數據輸出,我知道頻率和採樣率。
如何使用MATLAB計算常量a,b,c,p and q
的值?
我已衍生彈簧質量系統的組合的傳遞函數在s
域,它的形式爲計算傳遞函數係數
G = (as*s + bs +c)/(ps + q)
我有與位移涉及作爲輸入和力作爲測量數據輸出,我知道頻率和採樣率。
如何使用MATLAB計算常量a,b,c,p and q
的值?
您尚未提供足夠的信息以提供完整的答案(*),但聽起來您將要進行某種最小二乘曲線擬合。 fminsearch
將爲此工作,但有更好的選擇。我傾向於使用統計工具箱中的nlinfit
,但也有優化工具箱中的lsqcurvefit
和曲線擬合工具箱中的fit
。
(*)您的數據究竟是什麼?在已知頻率的正弦驅動力下,它是否作爲時間的函數而位移?但是,您必須擁有多個驅動頻率的數據,因爲您需要了解多個頻率的響應才能提取傳輸函數的所有5個參數,對不對?
感謝您的回覆。你是對的,我有超過1頻率的數據。這些數據是通過位移作爲輸入和力作爲時域輸出獲得的。所以我試圖對數據做一些簡單的處理,但是這個輸出是一個複數的列向量。我應該如何使用這個列向量來與理論導出的傳遞函數相關聯。也許我還沒有理解確切的概念。如果你能詳細闡述一下這件事,你真好。在此先感謝 – sat0408
我對tfestimate一無所知,所以如果您對此有特定疑問,我建議您打開一個新問題。但是這聽起來像是你真正的問題是關於如何解決這個特定實驗的數據分析問題。這實際上不是針對stackoverflow的主題,所以我建議在其他地方尋求幫助(教授/ ta是這是一個班級,同事如果這是工作等)。祝你好運! –
是的,這可以使用fminsearch
完成。欲瞭解更多信息,請閱讀function documentation
非常感謝您的即時回覆。爲了使用這個,我需要獲得上述形式的測量數據的傳遞函數,但是當我使用特徵函數時,我得到了一個複數的向量。你能請教我如何進一步進行。這將是非常有幫助的。 – sat0408
你確定這是你的系統模型嗎?你得出的模型是非因果的,這意味着它不能物理存在。假設你得出一個因果模型(例如1/G),你可以使用MATLAB系統識別UI。最簡單的方法是估算一個過程模型,您可以從那裏開始。 – HebeleHododo
順便說一下,非因果系統手段;系統的輸出取決於未來的輸入。因此在實踐中不可能存在。 – HebeleHododo
非常感謝回覆。想想我需要深入系統識別工具箱。關於因果關係的最後一個問題。你是否認爲這是非因果關係,因爲零的數量多於極點?從來沒有可能研究這樣的系統。按照我從彈簧阻尼器組合中得出的方式。再次感謝 – sat0408