Re: [請益]打上自己想要的函數

看板MATLAB作者 (成仔)時間13年前 (2012/10/21 23:14), 編輯推噓1(1018)
留言19則, 2人參與, 最新討論串2/2 (看更多)
: 題目: 畫出x from 0 to 30 的y(t)圖 : y(t)=sin((2/3)t))*exp(-0.1t) : : : 以下是我寫的: : : >>x=linspace(0,30); : >>y=sin((2/3)x))*exp(-0.1x); x是矩陣 所以要用 .* y= sin(2/3)*x).*exp(-0.1*x); : >>xlabel('Time(sec)'); : >>ylabel('Amplitude'); : >>print dmeta hw1 這邊用 plot (x,y) 即可 xlabel ylabel 在後 : 請問我的問題出在哪裡呢???QQ : : : 還有一個問題 就是如果我想要輸出y(x)的圖是: : g[x] 和 h[x] 的convolution sum的話 : 那我該怎麼打呢?? : 麻煩大家幫個忙了 謝謝!! g[x], h[x]是啥? rewirte x=linspace(0,30); y= sin(2/3)*x).*exp(-0.1*x); plot (x,y) xlabel('Time(sec)'); ylabel('Amplitude'); --- 我也是新手 @.@~ -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 111.255.179.239

10/21 23:29, , 1F
先謝謝您! ps g[x] h[x]是兩個不同的函數
10/21 23:29, 1F

10/22 10:32, , 2F
上面有人推 conv()
10/22 10:32, 2F

10/22 10:33, , 3F
做出 A= g(x) B = h(x) 下 conv(A,B) 應該可得
10/22 10:33, 3F

10/22 10:33, , 4F
A, B 是函數運作後產生數列
10/22 10:33, 4F

10/22 13:04, , 5F
好的!!!謝謝~~!
10/22 13:04, 5F

10/22 13:29, , 6F
我重打了一份如下:
10/22 13:29, 6F

10/22 13:30, , 7F
A=sin((1/2)*x);
10/22 13:30, 7F

10/22 13:31, , 8F
B=sin((1/5)*x);
10/22 13:31, 8F

10/22 13:31, , 9F
y=conv(A,B);
10/22 13:31, 9F

10/22 13:31, , 10F
plot (x,y)
10/22 13:31, 10F

10/22 13:32, , 11F
結果他說error using plot,vectors must be the same
10/22 13:32, 11F

10/22 13:32, , 12F
length. 怎麼會這樣子呢??QQ
10/22 13:32, 12F

10/22 15:02, , 13F
要先寫出你要x是多少 ex. x=linspace (0,30_);
10/22 15:02, 13F

10/22 15:03, , 14F
後面接下去 會產生 兩個一樣維度的矩陣 才可做conv
10/22 15:03, 14F

10/22 15:04, , 15F
最後才plot (x,y)
10/22 15:04, 15F

10/22 15:45, , 16F
恩恩!其實我有打在第一行...但是不行...QQ
10/22 15:45, 16F

10/22 21:53, , 17F
最後是plot (y) @@~ sorry 因為conv後會形成2x-1數據
10/22 21:53, 17F

10/22 21:54, , 18F
或是要將conv結果砍掉一半才能畫
10/22 21:54, 18F

10/22 22:20, , 19F
OK了!!!!太謝謝你了!!!^0^
10/22 22:20, 19F
文章代碼(AID): #1GX15Vzm (MATLAB)
文章代碼(AID): #1GX15Vzm (MATLAB)