《信号与系统》专题研究性学习实验报告.doc

信号与系统MATLAB专题研究性学习 专题一信号时域分析 1.基本信号的产生，语音的读取与播放 【研讨内容】 生成一个正弦信号，改变正弦信号的角频率和初始相位，观察波形变化； 生成一个幅度为1、基频为2Hz、占空比为50%的周期方波， 观察一定时期内的股票上证指数变化，生成模拟其变化的指数信号， 录制一段音频信号，进行音频信号的读取与播放 【题目分析】 (1) 正弦信号的形式为Acos（ω0t+ψ）或Asin（ω0t+ψ），分别用MATLAB的内部函数cos和sin表示，其调用形式为、 。生成正弦信号为y=5sin(t),再依次改变其角频率和初相，用matlab进行仿真。 (2) 幅度为1，则方波振幅为0.5，基频w0=2Hz，则周期T=pi，占空比为50%，因此正负脉冲宽度比为1。 (3) 将波形相似的某一段构造成一个指数函数，在一连续时间内构造不同的2~3个不同指数函数即可大致模拟出其变化。 (4) 录制后将文件格式转化为wav，再用wavread函数读取并播放，用plot函数绘制其时域波形。 【仿真】 正弦信号 正弦信号1： A=1;w0=1/4*pi;phi=pi/16; t=-8:0.001:8; xt1=A*sin(w0*t+phi); plot(t,xt1) title(xt1=sin(0.25*pi*t+pi/16)) 正弦信号2（改变1中频率） A=1;w1=1/4*pi;w2=1*pi;phi=pi/16; t=-8:0.001:8; xt1=A*sin(w1*t+phi); xt2=A*sin(w2*t+phi); plot(t,xt1,t,xt2) 正弦信号3（改变1中相位） A=1;w=1/4*pi;phi1=pi/16;phi2=pi/4; t=-8:0.001:8; xt1=A*sin(w*t+phi1); xt3=A*sin(w*t+phi2) plot(t,xt1,t,xt3) 方波信号 t=-100:0.01:100; T=0.5; f=1/T; y=square(2*pi*f*t,50); plot(t,y); axis([-2 2 -3 3]); 模拟股票上证指数变化的指数信号 x1=0:0.001:5; y1=2500+1.8*exp(x1); x2=5:0.001:10; y2=2847-1.5*exp(0.8*x2); x3=10:0.001:15; y3=2734+150*exp(-0.08*x3); x4=15:0.001:20; y4=2560-156*exp(-0.08*x4); x=[x1,x2,x3,x4]; y=[y1,y2,y3,y4]; plot(x,y); 音频信号的读取与播放 [x,Fs,Bits]=wavread(C:\Users\Ghb\Desktop\nansheng.wav) sound(x,Fs,Bits) plot(x) [x,Fs,Bits]=wavread(C:\Users\Ghb\Desktop\nvsheng.wav) sound(x,Fs,Bits) plot(x) 2.信号的基本运算（语音信号的翻转、展缩） 【研讨内容】 将原始音频信号在时域上进行延展、压缩， 将原始音频信号在频域上进行幅度放大与缩小， 将原始音频信号在时域上进行翻转， 【题目分析】 用matlab的wavread函数读取录制的音频，用length函数计算出音频文件的长度，最后计算出时间t，然后用plot函数输出录制的音频信号 （1） 延展与压缩分析 把时间t变为原来的一半，信号就被延展为原来的2倍，把时间他变为原来的2倍，信号就被压缩为原来的一半。 （2） 幅度放大与缩小 把信号的幅度值变为原来的2倍，信号的幅度就被放大了2倍，把信号的幅度变为原来的一半，信号就被缩小为原来的一半。 （3） 信号的翻转 把信号的t变为变为原来的相反数，就可以实现信号的翻转。 【仿真】 读取原始信号 fs=44100;bits=32; [x,fs,nbits]=wavread(C:\Users\Ghb\Desktop\nansheng.wav); plot(x);title(原始信号) wavplay(x,fs); 信号的延展 [x,fs,nbits]=wavread(C:\Users\Ghb\Desktop\nansheng.wav); x1=x(1:2:end) subplot(2,1,1);plot(x);title(原始信号) subplot(2,1,2);plot(x1);title(延展) wavplay(x1,fs); 信号的压缩 [x,fs,nbits]=wavread(C:\Users\Ghb\Desktop\