N=37-13; %定义信号频率N和M
M=N3;
t=(0:0.0001:0.1); %定义信号的时间范围
x=Nsin(2piNt)+Msin(2piM*t); %生成信号
fs=800; %采样频率为800Hz
A=80; %定义采样点数
dt=1/fs; %采样间隔,采样间隔其实就可以理解为是采样信号的周期,周期=1/频率
T=(0:A-1)dt; %定义采样的每个时间点
x1=Nsin(2piNt)+Msin(2piM*t); %对信号进行采样
subplot(311);
plot(t,x);
ylim([-1 1])
title('原始信号')
subplot(312)
plot(t,x,T,x,'rp');
ylim([-1 1]);
title('采样过程')
subplot(313)
plot(T,x1);
ylim([-1 1])
xlabel('时间/s')
title('采样后的信号')
work1
错误使用 plot
向量长度必须相同。
出错 work1 (line 22)
plot(T,x1);
1、对给定周期信号进行不同频率(至少两个不同频率)的采样,绘制采
样后的时域波形。
(周期信号为 f(t)=Nsin(2πNt)+Msin(2πMt),其中 N=37-自己序
号,M=N×3,M 和 N 分别为两个正弦信号的频率,单位为 Hz) 2、完成采样后所得到信号的频谱, 绘制频谱图。(注意纵坐标幅值的计
算标注和频率的计算和标注)
y轴不要去限制,然后给你加了fft变换:
N=37-13; %定义信号频率N和M
M=N*3;
t=(0:0.00001:0.1); %定义信号的时间范围
x=N*sin(2*pi*N*t)+M*sin(2*pi*M*t); %生成信号
fs=800; %采样频率为800Hz、300
dt=1/fs; %采样间隔,采样间隔其实就可以理解为是采样信号的周期,周期=1/频率
T=0:dt:0.1; %定义采样的每个时间点
x1=N*sin(2*pi*N*T)+M*sin(2*pi*M*T); %对信号进行采样
subplot(411);
plot(t,x, 'b-');
title('原始信号')
subplot(412)
plot(t,x,'b-',T,x1,'rp');
title('采样过程')
subplot(413)
plot(T,x1, 'r');
xlabel('时间/s')
title('采样后的信号')
subplot(414);
Y = fft(x1);
L = length(T);
P2 = abs(Y/L);
P1 = P2(1:floor(L/2+1));
P1(2:end-1) = 2*P1(2:end-1);
f = fs*(0:floor(L/2))/L;
plot(f,P1, 'b-')
title('Single-Sided Amplitude Spectrum of X(t)')
xlabel('f (Hz)')
ylabel('|P1(f)|')
800Hz采样频率
300Hz采样频率