[問題] 解微分方程式

看板MATLAB作者 (麒麟現在)時間14年前 (2011/11/01 20:46), 編輯推噓1(101)
留言2則, 2人參與, 最新討論串2/2 (看更多)
第一次用matlab就被擊倒了... 我想要用以下這個範例式 ------------------------------------------------- clear all; clc; n=[1:30]; y(1) = 6.16; y(2) = 8.77; for m=3:30; y(m)=0.1*y(m-1)+.72*y(m-2)+5; end x=fft(y); plot(n,x) ------------------------------------------------- 有以下幾個問題 1.y(m)=0.1*y(m-0.0002)+.72*y(m-0.0004)+5; 我需要在括號中放0.0002,但是他不讓我放非整數 是不是選錯範例的意思呢? 2.後來委曲求全把所有括號內的值都改成整數,但是他還是顯示 Attempted to access y(4); index out of bounds because numel(y)=2 是什麼意思呢? 是要我定義y(4)嗎?可是他不是應該幫我算嗎? 我目前就卡在這邊,希望版友好心幫我解惑 謝謝 -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.112.43.75

11/01 23:16, , 1F
矩陣的index當然要整數阿
11/01 23:16, 1F

11/03 04:59, , 2F
你應該要先開一個y = zeros(1,30)然後再填數字進去
11/03 04:59, 2F
文章代碼(AID): #1Eh-ejjI (MATLAB)
文章代碼(AID): #1Eh-ejjI (MATLAB)