[繪圖] 數位訊號繪圖

[zj44 (YY)]

這是使用內建spline
但是就是畫不出兩個圖
不知道出了甚麼bug....

clear; close all;

t = 0:0.001:1; xa = cos(20*pi*t);

Ts = 0.05; Fs = 1/Ts;

xa1=spline(Ts*n1,x1,t);

subplot(3,1,1);plot(t,xa1); axis([0,1,-1.1,1.1]);

ylabel('xa(t)'); title('Reconstruction of xa(t) when Ts=0.01');

Ts = 0.05; Fs = 1/Ts;

xa2 = spline(Ts*n2,x2,t);

subplot(3,1,2);plot(t,xa2); axis([0,1,-1.1,1.1]);

ylabel('xa(t)'); title('Reconstruction of xa(t) when Ts=0.05');

Ts = 0.1; Fs = 1/Ts;

xa3 = spline(Ts*n3,x3,t);

subplot(3,1,3);plot(t,xa3); axis([0,1,-1.1,1.1]);

ylabel('xa(t)'); title('Reconstruction of xa(t) when Ts=0.1');

xlabel('t,sec');

第二個使用

sinc函數

clear; close all;

t = 0:0.001:1; xa = cos(20*pi*t);

Ts = 0.01; Fs = 1/Ts;

xa1 = x1*sinc(Fs*(ones(length(n1),1)*t-(n1*Ts)'*ones(1,length(t))));

subplot(3,1,1);plot(t,xa1); axis([0,1,-1.1,1.1]);

ylabel('xa(t)'); title('Reconstruction of xa(t) when Ts=0.01');

Ts = 0.05; Fs = 1/Ts;

xa2 = x2*sinc(Fs*(ones(length(n2),1)*t-(n2*Ts)'*ones(1,length(t))));

subplot(3,1,2);plot(t,xa2); axis([0,1,-1.1,1.1]);

ylabel('xa(t)'); title('Reconstruction of xa(t) when Ts=0.05');

Ts = 0.1; Fs = 1/Ts;

xa3 = x3*sinc(Fs*(ones(length(n3),1)*t-(n3*Ts)'*ones(1,length(t))));

subplot(3,1,3);plot(t,xa3); axis([0,1,-1.1,1.1]);

ylabel('xa(t)'); title('Reconstruction of xa(t) when Ts=0.1');

xlabel('t,sec');

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.137.2.6
推 diva:是不是你第二張圖被close all關掉了
改掉還是一樣....
Error in ==> Untitled2222 at 4
xa1 = x1*sinc(Fs*(ones(length(n1),1)*t-(n1*Ts)'*ones(1,length(t))));

※ 編輯: zj44            來自: 140.137.2.6          (12/09 22:08)
推 flylee:x1 宣告在哪 ?
→ Rasin:n1, x1是什麼