《1、波形产生时间向量和正弦序列》由会员分享，可在线阅读，更多相关《1、波形产生时间向量和正弦序列（5页珍藏版）》请在金锄头文库上搜索。
1、1、波形产生：时间向量和正弦序列t = (0:0.001:1);generating data with a 1000 Hz sample frequencyy = sin(2*pi*50*t) + 2*sin(2*pi*120*t);randn(state,0);duplicate the results in the exampleyn = y + 0.5*randn(size(t);plot(t(1:50),yn(1:50)产生模拟序列产生离散序列（方法 1）n=1:100;y = sin(2*50*n) + 2*sin(2*120*n);y1=sin(2*50*n);stem(n,y)
2、figure(2);stem(n,y1)产生离散序列（方法 2）w=0:1:100*pi/100;%0,pi分成 101 份y1=sin(2*pi*1*w);可以显示高频信号改变 12500w1=0:1:500*pi/12500;y2=sin(2*50*w1) + 2*sin(2*120*w1);2、产生单位冲激函数(unit impulse)，单位阶跃函数(unit step)，单位斜坡函数（unit ramp）t = (0:0.001:1);y = 1; zeros(99,1); % impulsey = ones(100,1); % step (filter assumes 0 init
3、ial cond.)y = t; % rampy = t.2;y = square(4*t);All of these sequences are column vectors. The last three inherit their shapes from t用 plot(y)可以显示出上述模拟信号Generate discrete unit impulem=1:100;y=1 zeros(1,99);stem(m,y)generate discrete unit stepy1=ones(1,100);stem(m,y1)generate discrete unit ramp y=m;st
4、em(m,y)3、产生多通道函数t = (0:0.001:1);z = t t.2 square(4*t);均为列向量4、产生周期函数xtilde=x*ones(1,p);%x 是行向量， p 是周期xtilde=xtilde(:);xtilde=xtilde;5、Sinc 函数x = linspace(-5,5);y = sinc(x);plot(x,y)6、matlab 处理的数据的输入方式：（有待以后学习使用，暂时用不着）（1） 通过键盘直接输入（2） 通过 matlab 或工具箱的函数如 sin,cos,sawtooth,wquare 或 sinc。（3） 通过 matlab 的 lo
5、ad 命令导入 ASCII 或 MAT 文件。（4） 利用底层 I/O 命令 fopen、fread 和 fscanf 函数从 I/O 端口读入数据第二节 傅立叶操作和分析1、卷积（convolution ）（1）conv(1 1 1,1 1 1)ans =1 2 3 2 1（2）卷积有平滑波形的作用。相当于去除高频突变信号即实现了低通滤波器。由卷积得到的函数（f *g） （x ） ，一般要比 f，g 都光滑。特别当 g 为具有紧支集的光滑函数， f 为局部可积时，它们的卷积（f *g） （ x）也是光滑函数。利用这一性质，对于任意的可积函数 ，都可以简单地构造出一列逼近于 f 的光滑函数列
6、fs（x） ，这种方法称为函数的光滑化或正则化（http:/ 。时域的卷积对应于频域的乘积，所以卷积可以简化运算。（3）离散卷积 在数字信号处理和数字图像处理中常常用到离散卷积。 (f * g)(n) = sum(f(k)g(n k)（4）任意信号 与单位冲激信号 卷积的结果仍然是信号本身，即 ，任意信号 与一个延迟时间为 的单位冲激函数 相卷积的结果，相当于把信号本身延迟 ，即 2、传输函数H（z）是传输函数，b(i)和 a(i)是傅立叶系数，matlab 用行向量表示 a 和 b。3、 (1)filter 函数假设 a(1)1，滤波器在数字域可以表示为例子：b = 1; % Numerat
7、ora = 1 -0.9; % Denominatorx=1 zeros(1,99);y = filter(b,a,x);stem(y)(2)impz 函数b = 1; % Numeratora = 1 -0.9; % Denominatorh=impz(b,a,64);stem(h)4 、fvtool 工具的使用Filter Visualization tool(FVTool)，通过这个工具可以观察滤波器的频率响应，群延迟，冲激响应，阶跃响应，零极点分布和滤波器系数。语法：fvtool(b,a)fvtool(b1,a1,b2,a2,.bn,an) fvtool(Hd1,Hd2,.)h = fvtool(.)b=0.2 0.1;a=1 -0.4 -0.5;fvtool(b,a)可以看到幅度响应如下显示：冲激响应为阶跃响应为零极点分布为系数为
