电信0902石乐健.docVIP

姓名:石乐健 班级:电信0902 学号:200922240227 EX61 t=-1:0.001:1;%输出波形在-1到1之间，每隔0.01个单位取值 omega=2*pi;%角频率为2*pi y=square(2*pi*t,50);%产生矩形周期脉冲 plot(t,y),grid on%输出波形 xlabel(t),ylabel(周期方波信号)%x轴为t，y轴为‘周期方波信号’ axis([-1 1 -1.5 1.5])%输出波形x轴在-1到1之间，y轴在-1.5到1.5之间 n_max=[1 3 5 11 47];%分别去谐波数为1，3，4，11，47 N=length(n_max);%求n_max的长度 for k=1:N n=1:2:n_max(k); b=4./(pi*n); x=b*sin(omega*n*t); figure; plot(t,y); hold on; plot(t,x); hold off; xlabel(t),ylabel(部分和的波形)% x轴为t，y轴为‘部分和的波形’ axis([-1 1 -1.5 1.5]),grid on%输出波形x轴在-1到1之间，y轴在-1.5到1.5之间 title([最大谐波数=,num2str(n_max(k))])输出波形的标题为”最大谐波数=” end EX62 n=-30:30;tao=1;T=10;w1=2*pi/T;%输出波形在-30到30之间，周期为10，以及w1的表示形式 x=n*tao/T;fn=tao*sinc(x); subplot(311) stem(n*w1,fn),grid on%离散序列图 title(tao=1,T=10)%输出标题为tao=1，T=10 tao=1;T=5;w2=2*pi/T; x=n*tao/T;fn=tao*sinc(x); m=round(30*w1/w2); n1=-m:m; fn=fn(30-m+1:30+m+1); subplot(312) stem(n1*w2,fn),grid on title(tao=1,T=5) tao=2;T=10;w3=2*pi/T； x=n*tao/T;fn=tao*sinc(x); subplot(313) % 在标定位置上建立坐标系 stem(n*w3,fn),grid on% title(tao=2,T=10) EX631 t=-4:0.001:4;%输出波形在-4到4之间，每隔0.01个单位取值 omega=pi;角频率为pi y=sawtooth(0.5*pi*t,0.5);%输出周期三角波 y=abs(y);%对y去绝对值 plot(t,y),grid on%输出波形 xlabel(t),ylabel(三角波形)%输出波形的x轴为t，y轴为’三角波形’ axis([-2 2 -1.5 1.5])%输出波形的x取值范围在-2到2之间y取值范围在-1.5到1.5之间 n_max=[1 3 5 ]%谐波量分别去1，3，5 N=length(n_max);%求长度 for k=1:N n=1:2:n_max(k); b=4./(pi*n)./(pi*n); x=0.5+b*cos(omega*n*t); figure; plot(t,y); hold on; plot(t,x); hold off;%以上部分循环求y的值 label(t),ylabel(部分和的波形)%输出波形x轴为t，y轴为“部分和的波形” axis([-3 3 -1.5 1.5]),grid on%输出波形x轴在-3到3之间，y轴在-1.5到1.5之间 title([最大谐波数=,num2str(n_max(k))])%输出标题 end EX71 ft=sym(exp(-2*t)*Heaviside(t));%ft的表示形式 Fw=fourier(ft)%利用fourier函数求ft的傅立叶变换 Fw = 1/(2+i*w)%Fw的输出结果 EX72 syms t%定义符号表达式 Fw=sym(1/(1+w^2));%Fw的表示形式 ft=ifourier(Fw,t)%利用函数求Fw的傅立叶逆变换 ft = 1/2*exp(-t)*heaviside(t)+1/2*e