ݷʱ : 2024/11/17 4:31:29һ µȫĶ
ΪǼչϵ: 1 ԪֱԪ 1,2,3...г cos չϵ d=B_sym;disp(d) % d ΪǼչϵ: 2,3,4,...Ԫ 1,2,3...г sin չϵ t=-8*a:0.01:T-a; f1=c(1)+c(2).*cos(2*pi*1*t/5)+0.*sin(2*pi*1*t/5); ; % f2=c(3).*cos(2*pi*2*t/5)+0.*sin(2*pi*2*t/5); ; % 2 г f3=c(4).*cos(2*pi*3*t/5)+0.*sin(2*pi*3*t/5); % 3 г f4=c(5).*cos(2*pi*4*t/5)+0.*sin(2*pi*4*t/5); ; % 4 г f5=c(5).*cos(2*pi*5*t/5)+0.*sin(2*pi*5*t/5); % 6 г f6=c(7).*cos(2*pi*6*t/5)+0.*sin(2*pi*6*t/5); % 6 г f7=c(7).*cos(2*pi*7*t/5)+0.*sin(2*pi*7*t/5); % 7 г f8=c(8).*cos(2*pi*8*t/5)+0.*sin(2*pi*8*t/5); % 8 г f9=c(9).*cos(2*pi*9*t/5)+0.*sin(2*pi*9*t/5); % 9 г f10=c(10).*cos(2*pi*10*t/5)+0.*sin(2*pi*10*t/5); % 10 г f11=f1+f2+f3; % +2 г+3г f22=f11+f4+f5+f6; % +2 г+3г+4г+5г+6г f33=f22+f7+f8+f9+f10; % +2 г+3г+4г+5г+6г +7г+8г+8г+10г subplot(2,2,1) plot(t,f1) hold on y=time_fun_e(t) %ʱ亯-ھ plot(t,y,'r:') title('+ֱ') axis([-4,4.5,-0.5,1.5]) subplot(2,2,2) plot(t,f11) hold on y=time_fun_e(t) plot(t,y,'r:') title('+2 г+3г') axis([-4,4.5,-0.5,1.5]) subplot(2,2,3) plot(t,f22), hold on y=time_fun_e(t)