matlab動態模擬.....怎麼用MATLAB做一位置,每5筆矩陣資料又放到此位置 |
尚未結案
|
sharkneil
一般會員 發表:27 回覆:33 積分:12 註冊:2005-01-21 發送簡訊給我 |
我的問題是
我現在有62筆資料
如
x=[ 6.2016 0.4896 0 -0.2652 1.0000
5.9364 0.4896 -0.0204 0 1.0000
5.9364 0.4692 -0.0102 0 1.0000
5.9364 0.4590 -0.0306 0.0306 1.0000
...以下簡略
];
soc=[0
0.0152
0.0315
0.0477
0.0644
0.1805...也是62筆
]
y表示為x每五筆取一次的矩陣,如x1~x5為一矩陣,x2~x6為一矩陣,x3~x7為一矩陣,以此類推
c表示為soc每五筆取一次的矩陣,同上方式
製作一個動態的矩陣。 我現在要用到MATLAB做分析這些資料
假設我要每五筆與SOC做運算(a=inv(y)*c)
接下來f=y*a,並把f累積下來......並繪出曲線圖...
這怎麼用MATLAB寫
|
㊣
版主 發表:261 回覆:2302 積分:1667 註冊:2005-01-04 發送簡訊給我 |
ex:
x=[1 2 3 4 5 6 7 8 9 10]; y=[10 9 8 7 6 5 4 3 2 1]; for N=1:6, x1=[x(N) x(N 1) x(N 2) x(N 3) x(N 4)] y1=[y(N) y(N 1) y(N 2) y(N 3) y(N 4)] end 你試試看吧,這樣就把x,y的資料照你要的方式取出來了,再去做你要的運算即可.
------
------------------------------------------------------------------------- 走是為了到另一境界,停是為了欣賞人生;未走過千山萬水,怎知生命的虛實與輕重!? |
sharkneil
一般會員 發表:27 回覆:33 積分:12 註冊:2005-01-21 發送簡訊給我 |
引言: ex:x=[1 2 3 4 5 6 7 8 9 10]; y=[10 9 8 7 6 5 4 3 2 1]; for N=1:6, x1=[x(N) x(N 1) x(N 2) x(N 3) x(N 4)] y1=[y(N) y(N 1) y(N 2) y(N 3) y(N 4)] a=inv(x1)*y1 %得到a函數 f=a*x1 ^^^^^^^這個值f我要把它存起用,圖來顯示....我不知道怎麼寫! end 你試試看吧,這樣就把x,y的資料照你要的方式取出來了,再去做你要的運算即可. |
㊣
版主 發表:261 回覆:2302 積分:1667 註冊:2005-01-04 發送簡訊給我 |
|
sharkneil
一般會員 發表:27 回覆:33 積分:12 註冊:2005-01-21 發送簡訊給我 |
|
㊣
版主 發表:261 回覆:2302 積分:1667 註冊:2005-01-04 發送簡訊給我 |
引言: 現在因為我的資料是5筆為一組~~X軸要跟著變動~~譬如 FOR N=1:54, c=x*a % C為每次5筆的Y軸座標~~所以X軸也要做累計運算~怎麼寫好呢?? end把你的C與要畫的通通存到矩陣裡(c(N)),執行完迴圈再畫吧! 直接plot(c,f)就ok. c與f 緯度要一樣.
------
------------------------------------------------------------------------- 走是為了到另一境界,停是為了欣賞人生;未走過千山萬水,怎知生命的虛實與輕重!? |
sharkneil
一般會員 發表:27 回覆:33 積分:12 註冊:2005-01-21 發送簡訊給我 |
引言:我程式大概如下 for N=1:54, x1=[x(N,1) x(N 1,2) x(N 2,3) x(N 3,4) x(N 4,5)] x2=[x(N 1,1) x(N 2,2) x(N 3,3) x(N 4,4) x(N 5,5)] x3=[x(N 2,1) x(N 3,2) x(N 4,3) x(N 5,4) x(N 6,5)] x4=[x(N 3,1) x(N 4,2) x(N 5,3) x(N 6,4) x(N 7,5)] x5=[x(N 4,1) x(N 5,2) x(N 6,3) x(N 7,4) x(N 8,5)] xc=[x1; x2; x3; x4; x5] y1=[soc(N)] y2=[soc(N 1)] y3=[soc(N 2)] y4=[soc(N 3)] y5=[soc(N 4)] yc=[y1;y2;y3;y4;y5] a=inv(xc'*xc)*xc'*yc c(N)=xc*a loop 1:N plot(loop,C(N)) end 但是編譯時C(N)=xc*a就有問題了!!!好像不能這樣寫耶!!!引言: 現在因為我的資料是5筆為一組~~X軸要跟著變動~~譬如 FOR N=1:54, c=x*a % C為每次5筆的Y軸座標~~所以X軸也要做累計運算~怎麼寫好呢?? end把你的C與要畫的通通存到矩陣裡(c(N)),執行完迴圈再畫吧! 直接plot(c,f)就ok. c與f 緯度要一樣. |
本站聲明 |
1. 本論壇為無營利行為之開放平台,所有文章都是由網友自行張貼,如牽涉到法律糾紛一切與本站無關。 2. 假如網友發表之內容涉及侵權,而損及您的利益,請立即通知版主刪除。 3. 請勿批評中華民國元首及政府或批評各政黨,是藍是綠本站無權干涉,但這裡不是政治性論壇! |