课程设计 ：Matlab来对通信系统和通信信号进行仿真 对相关算法进行验证

%本程序使用两种方式计算信号的功率和能量 clear; close all; dt=0.01; t=0:dt:5; s1=exp(-5*t).*cos(20*pi*t); s2=cos(20*pi*t); Es1=sum(s1.*s1.*dt); %时域直接计算，幅度平方的积分 Ps2=sum(s2.*s2.*dt)/(length(t)*dt); %时域直接计算，先求能量，再除以长度 [f1,sf1]=T2F(t,s1); [f2,sf2]=T2F(t,s2); df=f1(2)-f1(1); F_Es1=sum(abs(sf1).^2)*df; %频域模值平方对频率求和 df=f2(2)-f2(1); F_Ps2=sum(abs(sf2).^2.*df)/(length(t)*dt); %频域模值平方对频率求和除以总频宽 fprintf('s1的能量为：%f %f\n',Es1,F_Es1); fprintf('s2的功率为：%f %f\n',Ps2,F_Ps2); figure(1); subplot(221); plot(t,s1,'k-','linewidth',1.5); axis([0 2 -1 1]);grid on; xlabel('t');ylabel('s1(t)');title('s1的时域波形'); subplot(223); plot(t,s2,'k-','linewidth',1.5); axis([0 2 -1 1]);grid on; xlabel('t');ylabel('s2(t)');title('s2的时域波形'); subplot(222); plot(f1,sf1,'k-','linewidth',1.5); %axis([0 2 -1 1]); grid on; xlabel('f');ylabel('sf1(f)');title('s1的频域波形'); subplot(224); plot(f2,sf2,'k-','linewidth',1.5); %axis([0 2 -1 1]); grid on; xlabel('f');ylabel('sf2(f)');title('s2的频域波形');