ex1_6.m

来自「数字信号处理 实验 的matlab 源码」· M 代码 · 共 33 行

% 1-6
% Q1.34
clear,clf
m=input('调制度：');
fH=input('载波频率：');
fL=input('调制信号频率：');
n=0:100;%m=0.4;fH=0.1;fL=0.01;
xH=sin(2*pi*fH*n);
xL=sin(2*pi*fL*n);
y=(1+m*xL).*xH; % m 为调制度
figure(1)
stem(n,y),grid on,xlabel('时间序号 n '),ylabel('振幅')


% Q1.36
%p1.7
clear
n=0:100;a=pi/200;b=0;
arg=a*n.^2+b*n;x=cos(arg);
figure(2)
subplot(211),stem(n,x)
axis([0 100 -1.5 1.5]),xlabel('时间序号 n '),ylabel('振幅')
%w=2an+b n=0 w=b=0;n=100 w=pi
title('扫频正弦信号 w_m_i_n= 0  w_m_a_x= pi'),grid on

% Q1.38
%w=2an+b n=0 w_min=0.1;n=100 w_max=0.3 -->a=0.001 b=0.1
n=0:100;a=0.001;b=0.1;
arg=a*n.^2+b*n;x=cos(arg);
subplot(212),stem(n,x)
axis([0 100 -1.5 1.5]),xlabel('时间序号 n '),ylabel('振幅')
%w=2an+b n=0 w=b=0;n=100 w=pi
title('扫频正弦信号 w_m_i_n= 0.1  w_m_a_x= 0.3'),grid on