第一节课代码解释.docVIP

clear all; close all; dt=0.01; %定义变量dt的值为0.01 t=0:dt:5; %定义变量t的值为从0到5每间隔0.01取值 s1=exp(-5*t).*cos(20*pi*t); %定义函数s1 s2=cos(20*pi*t); %定义函数s2 E1=sum(s1.*s1)*dt; %定义变量E1,s1的能量 P2=sum(s2.*s2)*dt/(length(t)*dt); %定义变量P2，s2的功率 [f1 s1f]=T2F(t,s1); %以f1，s1f为参量调用T2F函数 [f2 s2f]=T2F(t,s2); %以f2，s2f为参量调用T2F函数 df=f1(2)-f1(1); %定义变量df E1_f=sum(abs(s1f).^2)*df; %定义变量E1_f，s1f的能量 df=f2(2)-f2(1); %定义变量df T=t(end); %定义变量T的值为t变量的最后一个数，即5 P2_f=sum(abs(s2f).^2)*df/T; %定义变量P2，s2f的功率 figure(1) %新建一个显示窗口 subplot(2，1，1) %在显示窗口中创建一个2行1列的图形区域且 在第一行显示下面的函数图形 plot(t,s1) %显示s1函数图形 xlabel(t); %将x轴标记为t ylabel(s1(t)); %将y轴标记为s1（t） subplot(2，1，2) %在显示窗口中创建一个2行1列的图形区域且 在第二行显示下面的函数图形 plot(t,s2) %显示s2函数图形 xlabel(t); %将x轴标记为t ylabel(s2(t)); %将y轴标记为s2（t） ===================================================================== function[f,sf]=T2F(t,st) %声明函数 %this is a function using the FFT function to calculate a signals Fourier Translation %input is the time and the signal vectors,the length of time must greater than 2 %output is the frequency and the signal spectrum dt=t(2)-t(1); %定义变量dt，其值等于变量t的间隔 T=t(end) %定义变量T，其值为变量t的最后一个值，即5 df=1/T; %定义变量df，其值为1/T N=length(st); %定义变量N，其值为变量st的长度 f=-N/2*df:df:N/2*df-df; %定义变量f的取值范围为-N/2*df到 N/2*df-df，取值间隔为df sf=fft(st); %定义变量sf为对st做快速傅里叶变换 sf=T/N*fftshift(sf); %定义变量sf，并赋值为T/N*fftshift(sf) 修改抽样间隔和程序最后两行 %可以观察量化输入输出关系 n=input(请输入量化级数,k=[8]); if isempty(n), ? ? ? ? n=8; end %产生题目要求的函数 endtime = 2;? ?? ?? ?? ?%结束时间 dt = 0.1;? ?? ?? ?? ?? ?%抽样间隔 t = 0:dt:endtime; a = zeros(size(t));? ???%产生一个长度为N的全0序列 for tt=1:1/dt? ? ? ? a(tt+[0:endtime/2-1]*(2/dt))=t(tt);? ?%t when 0 = t 1 end for tt=1/dt+1:2/dt ? ? a(tt+[0:endtime/2-1]*(2/dt))=2-t(tt);??% -t+2??when 1 = t 2 end; %A率非线性 A = 87.56; amax=