信号与系统课设

一、

1.正弦信号

A = input('input A=') ;% 给正弦信号的幅度A赋值

w = input('input w=') ; % 给正弦信号的频率w赋值

theta =input('input theta='); % 给正弦信号的初始相位theta 赋值disp(['这个信号是周期信号'])

T=2*pi/w

t = 0 : 0.01 : 3*T ; % 定义时间点

ft = A * sin( w * t + theta ) ; % th计算函数值

plot( t ,ft ) ; % 画图

title( '正弦信号' ) ; % 为图像加标题注释

grid on ; % 在图上画方格

2.复指数信号

j00 = sqrt( - 1 ) ; % 定义复数j

a = input('input a='); % 复指数信号赋值w = input('input w=');

K = input('input K=');

if a==0

disp('这是一个周期信号')

T=2*pi/w

else if a>0

disp('这不是一个周期信号')

else

disp('这不是一个周期信号')

end

end

t = -1.5*abs(a) : 0.01 : 1.5*abs(a) % 定义时间点

ft = K*exp( ( a + j00 * w ) * t ) ;

subplot( 2 , 2 , 1 ) ; plot( t , real( ft ) ) ; title( '实部' ) ; %画图subplot( 2 , 2 , 2 ) ; plot( t , imag( ft ) ) ; title( '虚部' ) ;

subplot( 2 , 2 , 3 ) ; plot( t , abs( ft ) ) ; title( '模' ) ;

subplot( 2 , 2 , 4 ) ; plot( t , angle( ft ) ) ; title( '相角' ) ;

3.实指数信号

A = input('input A=') % 给指数信号中的A赋值

a = input('input a=') % 给指数信号中的a赋值

t = 0 : 0.01 : 10 ; % 定义时间点

ft = A * exp( a * t ) ; % 计算函数值

plot( t , ft ) ; % 以时间为横轴，函数值为纵轴画图title( '实指数信号' ) ; % 为图像加标题注释

grid on ; % 在图上画方格

三、信号的频谱图

w = 0 : 0.025 : 5 ;

b = input('请依次输入bn,bn-1...b0:');

a = input('请依次输入am,am-1...a0:');

H = freqs( b , a , w ) ;

subplot( 2 , 1 , 1 ) ;

plot( w , abs( H ) ) ;

grid on ;

xlabel( '\omega (rad/s)' ) ; ylabel( '|H(j\omega)|' ) ;

title( 'H(j\omega)的幅频特性' ) ; subplot( 2 , 1 , 2 ) ;

plot( w , angle( H ) ) ;

grid on ;

xlabel( '\omega (rad/s)' ) ; ylabel( '\phi(\omega)' ) ;

title( 'H(j\omega)的相频特性' ) ;

四、连续时间信号卷积

dt=0.01;t = - 1 : 0.01 : 4 ;

f1 =input('input f1=');

f2 =input('input f2=');

f = conv( f1 , f2 ) ; % 卷积n = length( f ) ;

tt = ( 0 : n - 1 ) * 0.01-2 ; subplot( 3 , 1 , 1 ) ;

plot( t , f1 ) ;

title( 'f1(t)' ) ;

subplot( 3 , 1 , 2 ) ;

plot( t , f2 ) ;

title( 'f2(t)' ) ;

subplot( 3 , 1 , 3 ) ;

plot( tt , f ) ;

title( 'f1(t)*f2(t)' ) ;

离散时间信号卷积

nx=-1:5;

nh=-2:10;

x=input('input f1=');

h=input('input f2=');

y=conv(x,h);

ny1=nx(1)+nh(1);

ny=ny1+(0:(length(nx)+length(nh)-2)); subplot(311)

stem(nx,x,'fill'),grid on

axis([-4 16 0 3])

xlabel('n'),title('x(n)')

subplot(312)

stem(nh,h,'fill'),grid on

xlabel('n'),title('h(n)')

axis([-4 16 0 3])

subplot(313)

stem(ny,y,'fill'),grid on

xlabel('n'),title('y（n）=x(n)*h(n)') axis([-4 16 0 3])

六、采样信号

Ts=1;

dt=0.1;

t1=-4:dt:4;

ft=input('input ft=');

subplot(221)

plot(t1,ft),grid on

axis([-4 4 -0.1 1.1])

xlabel('Time(sec)'),ylabel('f(t)') title('抽样前信号：')

N=500;

k=-N:N;

w=pi*k/(N*dt);

Fw=dt*ft*exp(-j*t1'*w); subplot(222)

plot(w,abs(Fw)),grid on

axis([-10 10 -0.2 1.1*pi]) xlabel('\omega'),ylabel('F(w)') title('采样前信号频谱')

t2=-4:Ts:4;

fst=input('input fst=');

subplot(223)

plot(t1,ft,':'),hold on

stem(t2,fst),grid on

axis([-4 4 -0.1 1.1])

xlabel('Time(sec)'),ylabel('fs(t)') title('采样后的信号：'),hold off Fsw=Ts*fst*exp(-j*t2'*w); subplot(224)

plot(w,abs(Fsw)),grid on

axis([-10 10 -0.2 1.1*pi]) xlabel('\omega'),ylabel('Fs(w)') title('采样信号频谱：')

七、拉普拉斯变换

f = input('input f=')

F = laplace( f ) % Laplace变换

b = input('以z的幂的降序输入分子系数b:') ; % 以z的幂的降序输入分子系数a =input('以z的幂的降序输入分母系数a:') ; % 以z的幂的降序输入分母系数[ z , p , k ] = tf2zp( b , a ) ; % 计算零极点和增益

plot( real( z ) , imag( z ) , 'o' ) ; % 画零点

hold on

plot( real( p ) , imag( p ) , 'x' ) ; % 画极点

xlabel( 'Re-axis' ) ;

ylabel( 'Im-axis' ) ;

title( 'zero-pole plot' ) ;