MATLAB 信号系统.docVIP

t1= -4*pi:0.1:4*pi; x1=(sinc(t1/pi)); plot(t1,x1); subplot(2,1,1); t2= -4*pi:0.1:4*pi; x2=(sinc(t2)); plot(t2,x2); subplot(2,1,2); %%%%%%%%%%%%%%%%%%%%%%%% t= -4*pi:0.1:4*pi; x1=(sinc(t/pi)); plot(x1); subplot(2,1,1); grid on t= -4*pi:0.1:4*pi; x2=(sinc(t)); plot(x2); subplot(2,1,2); grid on %显示网格 %%%%%%%%%%%%%%%%%%%%%%%%% syms w t; syms w t s; % sym w t syms w,t sym w sym t laplace(sin(w*t)) ans = w/(s^2+w^2) %%%%%%%%%%%%%%%%%%%%%%%%% ?syms t s w a; ?laplace(sin(w*t+a)) ans = cos(a)*w/(s^2+w^2)+sin(a)*s/(s^2+w^2) %%%%%%%%%%%%%%%%%%%%%% ?syms t s w a; ?laplace(sin(t)/t) ans = atan(1/s) %%%%%%%%%%%%%%%%%%%%%%%% ?syms t s w a; ?laplace(cos(t)/t) ans = laplace(cos(t)/t,t,s) %%%%%%%%%%%%%%%%%%% ?syms t s w a; %教材P250 4-2 ?laplace(sin(w*(t+a/w))) ans = cos(a)*w/(s^2+w^2)+sin(a)*s/(s^2+w^2) %%%%%%%%%%%%%%%%%%%%%%%%% ?syms t s w a; ?m=laplace(sin(w*t))*exp(s*a/w); ?m m = w/(s^2+w^2)*exp(s*a/w) ?simplify(m) ans = w/(s^2+w^2)*exp(s*a/w) ?ilaplace(m) ans = w*Heaviside(t+a/w)/(w^2)^(1/2)*sin((w^2)^(1/2)*(t+a/w)) %%%%%%%%%%%%%%%%%%%%%%%%% u（t）是单位阶跃信号t0时 u(t)=0 t0时 u(t)=1 heaviside.m %阶跃信号 function y = heaviside(t) y = t=0; %%%%%%%%%%%%%%%%%%%%%%%%% syms t f; f = sin(t); g = subs(f,t,-2*t) ezplot(g,[-3,0]) %%%%%%%%%%%%%%%%%%%%%% ?clear ?syms t f; f = sin(t); g = subs(f,t,-2*t) ezplot(g,[-3,0]) g = -sin(2*t) %%%%%%%%%%%%%%%%%%%%%%%% function y = heaviside(t) %单位阶跃信号 if t==0 y=0.5; elseif t0 y=-1; else y = 1; end %%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear a=[1,2,3]; b=[1,1]; sys=tf(b,a); t=[0:0.1:10]; e1=t.^2; r1=lsim(sys,e1,t); subplot(2,1,1); plot(r1); e2=exp(t); r2=lsim(sys,e2,t); subplot(2,1,2); plot(r2); %%%%%%%%见下图 %%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%% t=-4*pi:0.1/pi:4*pi; x1=(sinc(t/pi))/pi; x2=(sinc(3*t/pi)*3)/pi; x=x1.*sin(9*t)+x2.*cos(9*t); c1=2*cos(9*t).*sin(3*t); c2=pi*t; c=c1./c2; plot(t,c);grid on y=conv(c,x) %%%%%%%%%%%%%%%%%%% clear all; f=[0 0.19 0.2