%对变量定义为符号函数 syms t x=exp(-1*abs(t)); %定义双边指数函数 F=fourier(x); subplot(2,1,1) ezplot(x) subplot(2,1,2) ezplot(F) %对函数进行傅里叶变换 %定义二行一列第一个函数 %定义二行一列第二个函数 %对函数进行自动取值 %对函数进行自动取值 %定义二行一列第一个函数 %把h的模赋值给h1 %复数h的相角 %定义分子向量系数 b=[0 0 1]; %定义分母系数 a=[0.08 0.4 1]; [h,w]=freqs(b,a,100);%求频率函数 h1=abs(h); h2=angle(h); subplot(2,1,1) plot(w,h1); grid; xlabel('角频率'); ylabel('幅度'); title('H（jw）的幅频特性'); % 图形标题 subplot(2,1,2); plot(w,h2*180/pi); grid; xlabel('角频率'); ylabel('相位'); title('H（jw）的相频特性'); %图形标题 %输出幅频函数 % 显示格线 % x轴注解 % y轴注解 %输出相频函数 %显示格线 %x轴注解 %y轴注解 %定义二行一列第二个函数 指数函数的傅里叶变换 一个 RLC 二阶低通滤波器的幅频和 相频特性曲线 

%绘制原始信号 t=k/fs; N=1000;k=0:N-1; fm=10;fc=100;fs=1000;%定 义载波信号的频率 %定义周 期个数 %定义时 间间隔 x=sin(2.0*pi*fm*t); %定义 正弦函数 %定义四 行一列第一个函数 plot(t(1:200),x(1:200)) % 以二维形式输出 subplot(4,1,1) %产生’pm'调制信号 %调制信号 y=modulate(x,fc,fs,'pm'); subplot(4,1,2) %定义四行一列第二个函数 plot(t(1:200),y(1:200)) %以二维形式输出 xlabel('times(s)'); axis([0,0.2,-1,1]); 围 title('Modulated signal (pm)'); %图形标题 % x轴注解 %限定图形显示的范 正弦信号的调制 %调制信号 %产生’fm'调制信号 y=modulate(x,fc,fs,'fm'); subplot(4,1,3) %定义四行一列第三个函数 plot(t(1:200),y(1:200)) %以二维形式输出 xlabel('times(s)'); % x轴注解 axis([0,0.2,-1,1]); %限定图形显示的范围 title('Modulated signal (fm)'); %图形标题 %调制信号 %产生’am'调制信号 y=modulate(x,fc,fs,'am'); subplot(4,1,4) %定义四行一列第四个函数 plot(t(1:200),y(1:200)) %以二维形式输出 xlabel('times(s)'); % x轴注解 axis([0,0.2,-1,1]); %限定图形显示的范围 title('Modulated signal (am)'); %图形标题 

抽样与抽样定理的 matlab 实现 %定义模拟信号 %定义最大频率 %连续时间傅立叶变换 %频率从-Wmax到Wmax %Xa介于-Wmax至Wmax之间 %定义四行二列第一个函数 %以二维形式输出 %定义四行二列第二个函数 %以二维形式输出 %图形标题 %定义函数的周期和时间区间及增长速率 %定义离散信号 Dt=0.00005;t=-0.005:Dt:0.005;%定义时间区间及增长速率 xa=exp(-1000*abs(t)); Wmax=2*pi*2000; K=500;k=0:K;W=k*Wmax/K; Xa=xa*exp(-1j*t'*W)*Dt;Xa=real(Xa); % 连续时间傅立叶变换 W=[-fliplr(W),W(2:501)]; Xa=[fliplr(Xa),Xa(2:501)]; subplot(4,2,1); plot(t*1000,xa) subplot(4,2,2); plot(W/(2*pi*1000),Xa*1000) title('模拟信号') Ts=0.0002;n=-25:1:25; x=exp(-1000*abs(n*Ts)); w=pi*k/K;X=x*exp(-1j*n'*w); X=real(X); %连续时间傅立叶变换 w=[-fliplr(w),w(2:K+1)]; X=[fliplr(X),X(2:K+1)]; subplot(4,2,3); stem(n*Ts*1000,x) subplot(4,2,4); plot(w/pi,X) title('离散信号Ts=0.2毫秒') Ts=0.0005;n=-10:1:10; x=exp(-1000*abs(n*Ts)); w=pi*k/K;X=x*exp(-1j*n'*w);X=real(X);% 连续时间傅立叶变换 w=[-fliplr(w),w(2:K+1)]; X=[fliplr(X),X(2:K+1)]; subplot(4,2,5); stem(n*Ts*1000,x) subplot(4,2,6); %频率从-Wmax到Wmax %X介于-Wmax至Wmax之间 %定义四行二列第五个函数 %重构信号和原始信号之间的最大误差 %定义四行二列第六个函数 %重构信号和原始信号之间的最大误差 %定义四行二列第四个函数 %以二维形式输出 %图形标题 %定义函数的周期和时间区间及增长速率 %定义离散信号 %频率从-Wmax到Wmax %X介于-Wmax至Wmax之间 %定义四行二列第三个函数 

%以二维形式输出 %图形标题 %定义函数的周期和时间区间及增长速率 %定义离散信号 plot(w/pi,X) title('离散信号Ts=0.5毫秒') Ts=0.001;n=-5:1:5; x=exp(-1000*abs(n*Ts)); w=pi*k/K;X=x*exp(-1j*n'*w);X=real(X);% 连续时间傅立叶变换 w=[-fliplr(w),w(2:K+1)]; X=[fliplr(X),X(2:K+1)]; subplot(4,2,7); stem(n*Ts*1000,x) subplot(4,2,8); plot(w/pi,X) title('离散信号Ts=1毫秒') %频率从-Wmax到Wmax %X介于-Wmax至Wmax之间 %定义四行二列第七个函数 %重构信号和原始信号之间的最大误差 %定义四行二列第八个函数 %以二维形式输出 %图形标题