untitled.m

来自「信号合成」· M 代码 · 共 21 行

/*fs=5120;
dt=1.0/fs;
n=1024;
T=dt*n;
t=linspace(0,5*T,n);
y1=sin(2*pi*10*t);
y2=square(t);
y3=square(3*pi*10*t);
y4=sawtooth(t);
y5=sawtooth(2*pi*10*t);
y7=sawtooth(2*pi*5*t,0.5);
y6=y1+y2+y3+y4+y5+y7;
fig1=subplot(3,1,1);
plot(t,y1,t,y2,t,y3,t,y4,t,y5,t,y7);grid;
fig2=subplot(3,1,2);
plot(t,y6);grid;
p=fft(y6,2048);
pyy=2*sqrt(p.*conj(p))/n;
f=linspace(0,fs/2,n/2);
fig3=subplot(3,1,3);
plot(f,pyy(1:n/2));grid;*/