实验一：matlab实现序列卷积运算说课讲解

实验一：m a t l a b实现序列卷积运算

精品文档

收集于网络，如有侵权请联系管理员删除 ↑前言

MATLAB 是一套功能强大的工程计算及数据处理软件，广泛应用于工业，电子，医疗和建筑等众多领域。它是一种面向对象的，交互式程序设计语言，其结构完整又优良的可移植性。它在矩阵运

算，数字信号处理方面有强大的功能。另外，MATLAB 提供了方便的绘图功能，便于用户直观地输出处理结果。

本课程实验要求学生运用MATLAB 编程完成一些数字信号处理的基本功能，加深对教学内容的理解。

课程试验1——用MATLAB 实现序列卷积运算

一、实验目的

? 加深对常用离散信号的理解

? 熟悉并验证离散时间信号用数字序列表示的方法及序列的线性卷积运

算。

二、实验内容

1、（1）单位抽样序列

?

??=01)(n δ 00≠=n n 在MATLAB 中可以利用zeros()函数实现。

;

1)1();,1(==x N zeros x 如果)(n δ在时间轴上延迟了k 个单位，得到)(k n -δ即：

???=-01

)(k n δ 0

≠=n k n

实验四-使用matlab实现卷积的运算

一 实验目的 1、 学习MATLAB 语言的编程方法及熟悉MATLAB 指令； 2、 深刻理解卷积运算，利用离散卷积实现连续卷积运算； 二 实验内容 1、 完成)(1t f 与)(2t f 两函数的卷积运算 其中：)4()()(), ()(221--==-t u t u t f t u e t f t 在一个图形窗口中，画出)(1t f 、)(2t f 以 及卷积结果。要求每个坐标系有标题、坐标轴名称。 p = ; %定义时间间隔 t= 0:p:10; %定义时间向量 f1=exp(-2*t).*u(t); %将f （t ）表示出来 f2=u(t)-u(t-4); f=conv(f1,f2); subplot(1,2,1); plot(t,f1,t,f2); title('f1=e^-2t*u(t)'' / ''f2=u(t)-u(t-4)'); xlabel('t(sec)'); % 这行代码是给出x 坐标的标签 ylabel('f(t)'); grid on ； subplot(1,2,2); plot(f); title('f=f1*f2'); xlabel('t(sec)'); % 这行代码是给出x 坐标的标签 ylabel('f') grid on

2、 若系统模型为： )(3)()(4)(4)(' ' ' 't f t f t y t y t y +=++ 其中 )()(t u e t f t -= 求零状态响应，画出波形（函数本身画出一幅图，自己再画出一幅输入波形图）。 零状态响应: a= [1 4 4]; %将y （t ）各阶导数的系数放在向量a 中 b= [1 3]; %将f （t ）各阶导数的系数放在向量b 中 sys = tf(b, a); %求系统函数sys td = ; %定义时间间隔 t = 0 : td : 10; %定义时间向量 f = exp(-t).*u(t); %将f （t ）表示出来 y = lsim(sys, f, t); %求系统的零状态响应y plot(t, y); %绘出零状态响应的波形 xlabel('t(sec)'); % 这行代码是给出x 坐标的标签 ylabel('y(t)'); % 这行代码是给出y 坐标的标签 grid on

实验MATLAB符号运算

实验四 MATLAB 符号运算 一、实验目的 掌握符号变量和符号表达式的创建，掌握MATLAB 的symbol 工具箱的一些基本应用。 二、实验内容 (1) 符号变量、表达式、方程及函数的表示。 (2) 符号微积分运算。 (3) 符号表达式的操作和转换。 (4) 符号微分方程求解。 三、实验步骤 1. 符号运算的引入 在数值运算中如果求x x x πsin lim 0→，则可以不断地让x 接近于0，以求得表达式接近什么数，但是终究不能令0=x ，因为在数值运算中0是不能作除数的。MATLAB 的符号运算能解决这类问题。输入如下命令： >>f=sym('sin(pi*x)/x') >>limit(f,'x',0) >> f=sym('sin(pi*x)/x') f = sin(pi*x)/x >> limit(f,'x',0) ans = Pi 2. 符号常量、符号变量、符号表达式的创建 1) 使用sym( )创建 输入以下命令，观察Workspace 中A 、B 、f 是什么类型的数据，占用多少字节的内存空间。 >> A=sym('1') >> B=sym('x') >> f=sym('2*x^2+3*y-1') >> clear >> f1=sym('1+2') >> f2=sym(1+2) >> f3=sym('2*x+3') >> f4=sym(2*x+3) >> x=1 >> f4=sym(2*x+3) > A=sym('1') A = 1

>> B=sym('x') B = x >> f=sym('2*x^2+3*y-1') f = 2*x^2+3*y-1 >> clear >> f1=sym('1+2') f1 = 1+2 >> f2=sym(1+2) f2 = 3 >> f3=sym('2*x+3') f3 = 2*x+3 >> f4=sym(2*x+3) ??? Undefined function or variable 'x'. >> x=1 x = >> f4=sym(2*x+3) f4 =

matlab实现卷积运算

2、试求下列图片的卷积波形12()()f t f t * 2() f t t 1 -1 1() f t t 1 -1 列出编程步骤： p=0.01; k1=0:p:1; f1=ones(1,length(k1)); k2=-1:p:1; f2= (k2+1).*(k2<0)+(-k2+1).*(k2>=0); [f,k]=sconv(f1,f2,k1,k2,p) function [f,k]=sconv(f1,f2,k1,k2,p) 3、试求下列图片的卷积波形12()()f t f t *

1() f t t 1 0.5- 2() f t t 12 1 p=0.01; k1=-0.5:p:1; f1=ones(1,length(k1)); k2=0:p:2; f2= 0.5*k2; [f,k]=sconv(f1,f2,k1,k2,p) 4、试求下列图片的卷积波形12()()f t f t *

1() f t t 2 2 - 2() f t t 3-2 -3 21 p=0.01; k1=-2:p:2; f1= (k1==-2)+(k1==2); k2=-3:p:3; f2=(k2+3).*(k2<-2)+(-k2-1).*(k2>=-2).*(k2<=-1)+(k2-1).*(k2>=1).*(k2<=2)+(-k2+3).*(k2>2); [f,k]=sconv(f1,f2,k1,k2,p); 5、试求下列图片的卷积波形12()()f t f t *

1() f t t 5 -5 33 -2() f t t 3 -2 -3 21 p=0.01; k1=-10:p:10; f1=(k1>=-5).*(k1<=-3)+(k1>=3).*(k1<=5); k2=-3:p:3; f2=(k2+3).*(k2<-2)+(-k2-1).*(k2>=-2).*(k2<=-1)+(k2-1).*(k2>=1).*(k2<=2)+(-k2+3).*(k2>2); [f,k]=sconv(f1,f2,k1,k2,p);

利用MATLAB实现循环卷积.doc

一、实验目的 1.利用MATLAB 实现循环卷积。 2.比较循环卷积与线性卷积的区别。 二、实验条件 PC 机，MATLAB7.0 三、实验内容 1）循环卷积的定义：两个序列的N 点循环卷积定义为： )0()()()]()([1 0N n m n x m h n x n h N k N N <≤-=?∑-= 利用MATLAB 实现两个序列的循环卷积可以分三个步骤完成： (1)初始化：确定循环点数N ，测量输入2个序列的长度。 (2)循环右移函数：将序列x(n)循环右移，一共移N 次（N 为循环卷积的循环次数），最后将每次循环成的新序列组成一个矩阵V 。 (3)相乘：将x(n)移位后组成的矩阵V 与第二个序列h(n)对应相乘，即得循环卷积结果。程序如下： 程序一： clear;close all ; N=10; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1); stem(xxn1,x1); subplot(3,1,2); stem(xxn2,x2); x11=fft(x1,N);

x12=fft(x2,N); y11=x11.*x12; y1=ifft(y11,N); subplot(3,1,3); n=0:length(y1)-1; stem(n,y1,'.'); title('循环卷积的结果'); xlabel('n');ylabel('y1(n)'); 运行后所得图形如下： 观察所得的循环卷积结果发现并没有呈现周期性的序列，因此将程序做下列改变。程序二： clear;close all; N=40; x1=[6 15 -6 3 5 7 0 1]; x2=[7 1 2 9 4 3 20 6]; x2=[x2,x2,x2,x2]; xn1=length(x1); xxn1=0:xn1-1; xn2=length(x2); xxn2=0:xn2-1; subplot(3,1,1);

实验MATLAB符号计算

实验四符号计算 符号计算的特点：一，运算以推理解析的方式进行，因此不受计算误差积累问题困扰；二，符号计算，或给出完全正确的封闭解，或给出任意精度的数值解（当封闭解不存在时）；三，符号计算指令的调用比较简单，经典教科书公式相近；四，计算所需时间较长，有时难以忍受。 在MATLAB中，符号计算虽以数值计算的补充身份出现，但涉及符号计算的指令使用、运算符操作、计算结果可视化、程序编制以及在线帮助系统都是十分完整、便捷的。 MATLAB的升级和符号计算内核Maple的升级，决定着符号计算工具包的升级。但从用户使用角度看，这些升级所引起的变化相当细微。即使这样，本章还是及时作了相应的更新和说明。如MATLAB 6.5+ 版开始启用Maple VIII的计算引擎，从而克服了Maple V计算“广义Fourier变换”时的错误（详见第5.4.1节）。 5.1符号对象和符号表达式 5.1.1符号对象的生成和使用 【例5.1.1-1】符号常数形成中的差异 a1=[1/3,pi/7,sqrt(5),pi+sqrt(5)] % <1> a2=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) % <2> a3=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)],'e') % <3> a4=sym('[1/3,pi/7,sqrt(5),pi+sqrt(5)]') % <4> a24=a2-a4 a1 = 0.3333 0.4488 2.2361 5.3777 a2 = [ 1/3, pi/7, sqrt(5), 6054707603575008*2^(-50)] a3 = [ 1/3-eps/12, pi/7-13*eps/165, sqrt(5)+137*eps/280, 6054707603575008*2^(-50)] a4 = [ 1/3, pi/7, sqrt(5), pi+sqrt(5)] a24 = [ 0, 0, 0, 189209612611719/35184372088832-pi-5^(1/2)] 【例5.1.1-2】演示：几种输入下产生矩阵的异同。 a1=sym([1/3,0.2+sqrt(2),pi]) % <1> a2=sym('[1/3,0.2+sqrt(2),pi]') % <2> a3=sym('[1/3 0.2+sqrt(2) pi]') % <3> a1_a2=a1-a2 % a1 = [ 1/3, 7269771597999872*2^(-52), pi] a2 = [ 1/3, 0.2+sqrt(2), pi] a3 = [ 1/3, 0.2+sqrt(2), pi] a1_a2 = [ 0, 1.4142135623730951010657008737326-2^(1/2), 0]

基于Matlab实现线性卷积等

线性卷积与循环卷积 一、作品目的 通过matlab的强大功能展示线性卷积和循环卷积过程中方方面面的计算和变化，让大家对这两种卷积有一个更加完美的认识。 二、概念简介 卷积是一种典型的乘累加运算。 1.线性卷积 线性卷积是对线性移不变（LSI）系统的输入输出关系的描述，体现系统的特性。 线性卷积的表达式为 一般情况，现实的系统为因果系统，有k<0时，恒有h(k)=0，则 若x(n)是一个N点序列，h(n)是一个m点序列，则卷积的结果y(n)将是L=N+M-1点的序列。 2.循环卷积

设x1(n) 和x2(n) 是两个长度为L、M的有限长序列，它们的N 点循环卷积x3(n) 定义为： 注意：其中N>=Max{L,M}如果其中一个序列(或者两个序列）的长度没有所求N点循环卷积的长度长，那在该序列后面补零，直到长度达到N。 三、设计思路及程序 1. 线性卷积： （1）以输入序列x（n）=[5,4,3,2,1]，脉冲响应h（n）=[1,1,1,1]为列进行演示。 （2）计算输入序列和脉冲响应的长度。 （3）画出补零后的输入序列和脉冲响应 （4）设计一个循环，在循环中实现反转、位移和计算。并画出反转后的图像变化和卷积图像，将每一次移位结果保存为fig图。（5）最后将上一步所生成的所有fig图合起来生成一张gif图 程序展示： clear; clc; close all; (1)(2)

xn=[5,4,3,2,1]; M=length(xn);%输入任意序列并计算长度M hn=[1,1,1,1]; N=length(hn);%输入任意脉冲响应并计算长度N m=[-(M-1):M+N-2];%设置代换变量的范围以便x(m)翻转和移位(3) xm=[zeros(1,M-1),xn,zeros(1,N-1)];%补零以便与m对应绘图 subplot(2,2,1);stem(m,xm,'r.');%%绘输入序列x(m) ylabel('x(m)'); grid on; title('(a)输入序列x(m)'); hm=[zeros(1,M-1),hn,zeros(1,M-1)];%补零以便与m对应绘图 subplot(2,2,2);stem(m,hm,'r.');%绘脉冲响应 ylabel('h(m)'),grid,title('(b)脉冲响应h(m)');%%加标签网格和标题 yn=zeros(1,2*M+N-2);%卷积输出初始化 (4) for n=0:M+N-2;%逐个计算卷积输出 if n==0; xmfy=[fliplr(xn),zeros(1,M+N-2)];%实现翻转 else for k=M:-1:1;

用matlab实现两个离散序列的卷积(不使用conv函数)

作业2.用matlab实现离散序列的卷积. N=14; n=[1:N-1]; f=1/16; signal1=5*sin(2*pi*n/8); figure(1); subplot(3,1,1) stem(n,signal1);title( ' 信号1' );xlabel( 'n' );ylabel( axis([0 15 -6 6]) long_M=5; signal2=ones(1,long_M); subplot(3,1,2) stem(signal2);title( ' 信号2' );xlabel( 'n' );ylabel( axis([0 6 -2 2]); grid on; long_N=length(signal1); fk=zeros(0,long_N+long_M+10); if (long_N>long_M) for k=1:1:long_N+long_M-1 a=0; if (k<=long_N) for i=1:1:k if (i>long_M) fk(k)=a; else fk(k)=a+signal2(i)*signal1(k-i+1); a=fk(k); end end else for i=1:1:k if (k-long_N+i>long_M) fk(k)=a; else fk(k)=a+signal2(k-long_N+i)*signal1(long_N-i+1); a=fk(k); end end end end end subplot(3,1,3) stem(fk);title( ' 卷积函数的实现' );xlabel( 'n' );ylabel( 'y(n)' ); 'y(n)' ); 幅度' );

MATLAB符号计算实验报告

实验六符号计算 学院：数计学院班级：1003班姓名：黄晓丹学号：1051020144 一、实验目的 1、了解富符号对象和数值对象之间的差别，以及它们之间的互相转换 2、了解符号运算和数值运算的特点、区别和优缺点 3、掌握符号对象的基本操作和运算，以及符号运算的基本运用 二、实验内容 1、符号常数形成和使用 （1）符号常数形成中的差异 >> a1=[1/3,pi/7,sqrt(5),pi+sqrt(5)] a1 = 0.3333 0.4488 2.2361 5.3777 >> a2=sym([1/3,pi/7,sqrt(5),pi+sqrt(5)]) a2 = [ 1/3, pi/7, sqrt(5),

6054707603575008*2^(-50)] >> a3=sym('[1/3,pi/7,sqrt(5),pi+sqrt(5)]') a3 = [ 1/3, pi/7, sqrt(5), pi+sqrt(5)] >> a24=a2-a3 a24 = [ 0, 0, 0, 189209612611719/35184372088832-pi-5^(1/2)] （2）把字符表达式转化为符号变量 >> y=sym('2*sin(x)*cos(x)') y = 2*sin(x)*cos(x) >> y=simple(y)

y = sin(2*x) （3）用符号计算验证三角等式 >> syms fai1 fai2;y=simple(sin(fai1)*cos(fai2)-cos(fai1)*sin(fai2)) y = sin(fai1-fai2) （4）求矩阵的行列式值、逆和特征值 >> syms a11 a12 a21 a22;A=[a11,a12;a21,a22] A = [ a11, a12] [ a21, a22] >> DA=det(A),IA=inv(A),EA=eig(A) DA =

MATLAB实现卷积码编译码-

本科生毕业论文（设计） 题目：MATLAB实现卷积码编译码 专业代码： 作者姓名： 学号： 单位： 指导教师： 年月日

目录 前言----------------------------------------------------- 1 1. 纠错码基本理论---------------------------------------- 2 1.1纠错码基本理论 ----------------------------------------------- 2 1.1.1纠错码概念 ------------------------------------------------- 2 1.1.2基本原理和性能参数 ----------------------------------------- 2 1.2几种常用的纠错码 --------------------------------------------- 6 2. 卷积码的基本理论-------------------------------------- 8 2.1卷积码介绍 --------------------------------------------------- 8 2.1.1卷积码的差错控制原理----------------------------------- 8 2.2卷积码编码原理 ---------------------------------------------- 10 2.2.1卷积码解析表示法-------------------------------------- 10 2.2.2卷积码图形表示法-------------------------------------- 11 2.3卷积码译码原理---------------------------------------------- 15 2.3.1卷积码三种译码方式------------------------------------ 15 2.3.2V ITERBI译码原理---------------------------------------- 16 3. 卷积码编译码及MATLAB仿真---------------------------- 18 3.1M ATLAB概述-------------------------------------------------- 18 3.1.1M ATLAB的特点------------------------------------------ 19 3.1.2M ATLAB工具箱和内容------------------------------------ 19 3.2卷积码编码及仿真 -------------------------------------------- 20 3.2.1编码程序 ---------------------------------------------- 20 3.3信道传输过程仿真-------------------------------------------- 21 3.4维特比译码程序及仿真 ---------------------------------------- 22 3.4.1维特比译码算法解析------------------------------------ 23 3.4.2V ITERBI译码程序--------------------------------------- 25 3.4.3 VITERBI译码MATLAB仿真----------------------------------- 28 3.4.4信噪比对卷积码译码性能的影响 -------------------------- 28

用MATLAB实现序列圆周卷积

数字信号处理实验报告 实验项目名称：用MATLAB实现序列的圆周卷积 实验日期： 2012-11-28 实验成绩： 实验评定标准： 一、实验目的 通过本实验，掌握一些基本而且重要的离散时间信号，熟悉基本离散时间信号的MATLAB实现方法。 二、实验器材 PC机，MATLAB软件。 三、实验内容 计算两序列x1（n）={1,2,3,4,5}，x2（n）={1,2,3,4,5,4,3,2,1}的圆周卷积。 四、实验结果 实验代码： clear all close all clc x1=[1,2,3,4,5,6,7,8]; x2=[1,2,3,4,5,6,7,8,7,6,5,4,3,2, 1]; N=length(x1)+length(x2); n=0:N-1 n1=0:N-2; n2=0:N-3; y1=circonvt(x1,x2,N); y2=circonvt(x1,x2,N-1);

y3=circonvt(x1,x2,N-2); x1=[x1 zeros(1,N-length(x1))]; x2=[x2 zeros(1,N-length(x2))]; Xf1=dft(x1,N); Xf2=dft(x2,N); Xf=Xf1.*Xf2; x=idft(Xf,N); x=real(x); subplot(2,3,1) stem(n,x1); title('x1(n)'); subplot(2,3,2) stem(n,x2); title('x2(n)') subplot(2,3,3); stem(n,x); title('x(n)=IDFT(X(k))'); subplot(2,3,4); stem(n,y1); title('N点圆周卷积'); subplot(2,3,5); stem(n1,y2); title('N-1点圆周卷积'); subplot(2,3,6); stem(n2,y3); title('N-2点圆周卷积'); function y=circonvt(x1,x2,N) if length(x1)>N error('N 必须 >= x1的长度') end if length(x2)>N error('N 必须 >= x2的长度') end x1=[x1 zeros(1,N-length(x1))]; x2=[x2 zeros(1,N-length(x2))]; m=[0:1:N-1]; x2=x2(mod(-m,N)+1); H=zeros(N,N); for n=1:1:N H(n,:)=cirshift(x2,n-1,N); end y=x1*H; function y=cirshift(x,m,N) if length(x)>N error('N 必须 >= x的长度') end x=[x zeros(1,N-length(x))]; n=[0:1:N-1]; n=mod(n-m,N); y=x(n+1); function [Xk]=dft(xn,N) n=[0:1:N-1];k=[0:1:N-1]; WN=exp(-j*2*pi/N); nk=n'*k; WNnk=WN.^nk; Xk= xn * WNnk; function [xn]=idft(Xk,N) %计算逆离散傅里叶变换 %[xn]=idft(Xk,N) n=[0:1:N-1];

matlab实验五多项式和符号运算

实验五：Matlab多项式和符号运算 一、实验目的 1．掌握Matlab多项式的运算。 2．了解符号运算。 二、实验内容 1．将多项式()(2)(3)(7)(1) =-+-+化为x的降幂排列。 P x x x x x syms x; y=(x-2)*(x+3)*(x-7)*(x+1); expand(y) ans = x^4-5*x^3-19*x^2+29*x+42 2．求一元高次方程的根。 98765432 --++--++= 53015027313658204100576-28800 x x x x x x x x x syms x y; y=x^9-5*x^8-30*x^7+150*x^6-1365*x^4-820*x^3+410 0*x^2+576*x-2880; solve(y,x) ans = 6.81947687944124431946 1.42761488953013276419+.8192491831*i 2.865487219+2.49263348244446271927*i

-1.887673354+1.812452594*i -.9583509633 -5.922730991 -1.887673354-1.812452594*i 2.865487219-2.49263348244446271927*i 1.42761488953013276419-.8192491831*i 3．求一元高次方程的根，并画出左边多项式函数在[2,2] x∈-区间内的曲线。 42 -+= x x 210 a=[1 0 -2 0 1]; r=roots(a) syms x; x=-2:2; y=[1 0 -2 0 1]; plot(x,y) r = 1.0000 + 0.0000i 1.0000 - 0.0000i -1.0000 -1.0000

用MATLAB实现线性卷积运算

北京邮电大学 实验报告 实验名称：用MATLAB实现线性卷积运算学院：信息与通信工程学院 班级： 姓名： 学号： 日期：2012年5月

一、实验原理 1、算法产生背景 DFT 是连续傅里叶变换在时域和频域上都离散的形式，将时域信号的采样变换为在离散时间傅里叶变换频域的采样。在形式上，变换两端（时域和频域上）的序列是有限长的。DFT 具备明确且合理的物理含义，适合应用于数字系统，同时可以方便地由计算机进行运算。对于线性非移变离散系统，可由线性卷积表示时域输入输出关系，即 () ()*() ynxnhn 通常采用循环卷积降低运算量，但实际中往往无法满足对信号处理的实时性要求。因此，产生了重叠相加法和重叠保留法两种典型的算法，用以快速计算线性卷积，成为了DFT 的一个重要应用。 2、算法基本思想 1）重叠相加法 重叠相加法是将待过滤的信号分割成长为N 的若干段，如图1 所示，每一段都可以和有限时宽单位取样响应作卷积，再将过滤后的各段重叠相加。 具体算法实现原理如图2 所示，建立缓存序列，每次输入N 点序列，通过计算x(n) 和h(n) 的循环卷积实现线性卷积运算，将缓存的M-1 点序列和卷积结果相加，并输出前N 点作为计算结果，同时缓存后M-1 点，如此循环，直至所有分段计算完毕，则输出序列y(n)为最终计算结果。

2）重叠保留法 重叠保留法相当于将x l(n)和h(n)作循环卷积，然后找出循环卷积中相当于线性卷积的部分。在这种情况下，将序列y(n)分为长为N的若干段(如图3所示)，每个输入段和前一段有M-1个重叠点。此时只需要将发生重叠的前M-1个点舍去，保留重叠的部分并输出，则可获得序列y(n)，算法如图4所示。

MATLAB计算卷积

本科实验报告 实验名称：MATLAB计算卷积 学员：学号： 培养类型：技术类年级： 2008级 专业：电子工程所属学院：电子科学与工程学院指导教员：职称：教授 实验室：实验日期： 2010年12月23日 国防科学技术大学训练部制

实验一 MATLAB 计算卷积 一. 实验目的 通过MATLAB 卷积程序的编写, 学会使用计算机软件编程来处理信号(即信号的采样与计算)并进一步了解信号卷积计算的方法与步骤,以更好地理解信号与系统课程知识理论，加强理论与实践的结合。 二. 实验要求 使用MATLAB 编写主程序和卷积函数，实现两个信号的卷积运算。要求： 1. 程序必须对任意两个函数信号都适用； 2. 结果必须是至少两个信号的卷积。 三. 实验平台 笔记本电脑 MATLAB 软件平台 四. 实验原理 信号卷积的计算公式为 可通过图形变换的方法来计算两个信号的卷积。 用图解法求解卷积的步骤是：翻转、滑动、相乘、积分。 a. 确定卷积结果的分段时限； b. 确定每段中积分的上下限； c. 确定每段中积分函数的表达式具体实现如下图 ()()()()()c t f t g t f g t d τττ ∞-∞ =*= -?

五. 实验内容 对于信号 计算其卷积 六. 实验分析 （1）根据题目要求输入信号，观察信号和卷积结果的波形表示。 信号f 的波形为 2000007 6 6 02 0 0 other Pair 1: = 15 w here 1010 1010 2010ππτ τττ-?+≤≤=? ?-+-==?=?=?* exp()()()()(.)()() .,,./f t kt t s t s t s t s t g t s t f H z s k M H z s

matlab符号计算实验报告

1. 已知x=6,y=5, 利用符号表达式求z =>> syms x >> z=(x+1)/(sqrt(x+3)-sqrt(y)); >> subs(z,x,5) ans =6/(8^(1/2)-y^(1/2)) >> subs(ans,6) ans = 15.8338 2. 分解因式。 （1）x y -44； >> syms x y >> factor(x^4-y^4) ans =(x-y)*(x+y)*(x^2+y^2) （2）x x x +++642 12575151 >> syms x >> factor(125*x^6+75*x^4+15*x^2+1) ans =(5*x^2+1)^3 3. 化简表达式 （1）sin cos cos sin ββββ-1212； >> syms x y >> f=sin(x).*cos(y)-cos(x).*sin(y)； >> sfy1=simple(f) 结果：sfy1 =sin(x-y) （2）x x x +++248321 >> syms x >> f=(4*x^2+8*x+3)/(2*x+1);sfy1=simplify(f) sfy1 =2*x+3 4、求下列极限，将完成实验的程序写到文件sy1.m 中： （1） （2） （3） （4） （5） （1）>> syms x >> F1=atan(x)/(x); >> w=limit(F1) w =1 （2）>> syms x F2=((1+x)/(1-x))^(1/x); >> w=limit(F2) w =exp(2) （3）>> syms x F3=(x.*log(1+x))/(sin(x^2)); >> w=limit(F3) w =1 （4）>> syms x F4=atan(x)/(x); >> w=limit(F4,x,inf) w =0 （5）>> syms x F5=(1/(1-x)-1/(1-x^3)); >> w=limit(F5,x,1) w =NaN 5、求下列函数的导数，将完成实验的程序写到文件sy2.m 中： 1、 >> x = sym('x'); >> y1=(cos(x))^3-cos(3*x); >> diff(y1)ans =-3*cos(x)^2*sin(x)+3*sin(3*x) 2、 >> x = sym('x'); >> y2=x.*sin(x).*(log(x)); >> diff(y2)ans =sin(x)*log(x)+x*cos(x)*log(x)+sin(x) 3、 >> x = sym('x'); >> y3=(x.*exp(x)-1)/sin(x); >> diff(y3) ans =(exp(x)+x*exp(x))/sin(x)-(x*exp(x)-1)/sin(x)^2*cos(x) 4、 x x x x F 1011lim 2??? ??-+=→3 1115lim()11x F x x →=---20sin )1ln(lim 3x x x F x +=→x x F x arctan lim 10→=arctan 4lim x x F x →∞=x x y 3cos cos 13-=x x x y ln sin 2=x xe y x sin 13-=cos x y e x =

实验二优选资料连续时间信号卷积运算的MATLAB实现

实验二连续时间信号卷积运算的M A T L A B实现 一．实验目的 （1）理解掌握卷积的概念及物理意义。 （2）理解单位冲激响应的概念及物理意义。 二．实验原理 三.实验参考程序 用MATLAB实现连续信号f1（t）和f2（t）卷积。首先利用MATLAB实现连续信号卷积的通用函数sconv(): function[f,k]=sconv(f1,f2,k1,k2,p); f=conv(f1,f2); f=f*p; k0=k1(1)+k2(1); k3=length(f1)+length(f2)-2; k=k0:p:k3*p; subplot(2,2,1); plot(k1,f1); title('f1(t)'); xlabel('t'); ylabel('f1(t)'); subplot(2,2,2); plot(k2,f2); title('f2(t)'); xlabel('t'); ylabel('f2(t)'); subplot(2,2,3); plot(k,f); h=get(gca,'position'); h(3)=*h(3); set(gca,'position',h); title('f(t)=f1(t)*f2(t)'); xlabel('t'); ylabel('f(t)'); 例2-1 已知两连续时间信号，试用MATLAB求f（t）=f1（t）*f2（t），并绘出f（t）的时域波形图。 实现上述过程的MATLAB命令如下： p=; k1=0:p:2; f1=*k1;

k2=k1; f2=f1; [f,k]=sconv(f1,f2,k1,k2,p); 运行程序结果图： 而当p=时连续时间信号的波形图为 对比上面两幅图可见,当抽样时间p 足够小时,函数sconv()的计算结果就是连续时间卷积 )(*)()(21t f t f t f =的较好的数值近似。 四.实验要求 在课本卷积部分课后习题中任选两题，完成信号的卷积运算，绘出波形。 已知两连续时间信号如图所示，使用MATLAB 求()()()12f t f t f t =*，并绘出()f t 的时域波形图。 五．实验心得 通过本次实验我学会使用MATLAB 软件，并用软件实现了的连续时间信号的卷积运算，且进一步理解了卷积的概念和物理意义，对我以后学习卷积知识有了很大的帮助。

实验五 使用matlab实现卷积的运算

实验五 使用matlab 实现卷积的运算 一 实验目的 1、 学习MATLAB 语言的编程方法及熟悉MA TLAB 指令； 2、 深刻理解卷积运算，利用离散卷积实现连续卷积运算； 二 实验内容 1、 完成)(1t f 与)(2t f 两函数的卷积运算 其中：)4()()(), ()(221--==-t u t u t f t u e t f t 在一个图形窗口中，画出)(1t f 、 )(2t f 以及卷积结果。要求每个坐标系有标题、坐标轴名称。 >> p=0.1; t=0:p:10; f1=exp(-2*t).*u(t); f2=u(t)-u(t-4); f=conv(f1,f2); subplot(1,3,1); plot(t,f1,'r'); title('f1(t)=e^-2*t*u(t)'); xlabel('t(sec)'); ylabel('f1(t)'); subplot(1,3,2); plot(t,f2,'g'); title('f2(t)=u(t)-u(t-4)'); xlabel('t(sec)'); ylabel('f2(t)'); subplot(1,3,3); plot(f); title('f(t)=f1(t)*f2(t)'); xlabel('t(sec)'); ylabel('f(t)');

05100 0.10.20.30.40.5 0.6 0.70.8 0.9f1(t)=e - 2*t*u(t)t(sec) f 1(t ) 510 0.10.20.30.4 0.50.6 0.70.80.9 1f2(t)=u(t)-u(t-4) t(sec) f 2(t ) 0200400 0.511.522.53 3.54 4.5 5f(t)=f1(t)*f2(t) t(sec) f (t ) 2、 若系统模型为： )(3)()(4)(4)(' ' ' 't f t f t y t y t y +=++ 其中 )()(t u e t f t -= 求零状态响应，画出波形（函数本身画出一幅图，自己再画出一幅输入波形图）。 a=[1 4 4]; b=[1 3]; sys=tf(b,a); td=0.01; t=0:td:10; f=exp(-t).*u(t); y=lsim(sys,f,t); plot(t,y); xlabel('t(sec)'); ylabel('y(t)');

实验3 Matlab 符号运算及求函数极值

实验3 Matlab 符号运算及求函数极值一、实验目的和要求 掌握用Matlab软件进行符号运算以及求函数的极值。 二、实验环境 Windows系列操作系统，Matlab软件。 三、实验内容 1.用MATLAB进行符号运算； 2.编程求函数的极值。 四、实验步骤 3.开启软件平台——Matlab，开启Matlab编辑窗口； 4.根据求解步骤编写M文件； 5.保存文件并运行； 6.观察运行结果(数值或图形)； 7.根据观察到的结果和体会写出实验报告。 五、示例 1．计算一元函数的极值 例1求 2 2 344 1 x x y x x ++ = ++ 的极值 解首先建立函数关系： s yms x y=(3*x^2+4*x+4)/( x^2+x+1); 然后求函数的驻点： dy=diff(y); xz=solve(dy) xz= [0] [-2] 知道函数有两个驻点x 1=0和x 2 =-2， 接下来我们通过考察函数的图形，则它的极值情况和许多其它特性是一目了然的。而借助MATLAB的作图功能，我们很容易做到这一点。 例2 画出上例中函数的图形

解 syms x y=(3*x^2+4*x+4)/( x^2+x+1); 得到如下图形 ezplot(y) 2．计算二元函数的极值 MATLAB 中主要用diff 求函数的偏导数,用jacobian 求Jacobian 矩阵。 例1 求函数42823z x xy y =-+-的极值点和极值. 首先用diff 命令求z 关于x,y 的偏导数 >>clear; syms x y; >>z=x^4-8*x*y+2*y^2-3; >>diff(z,x) >>diff(z,y) 结果为 ans =4*x^3-8*y ans =-8*x+4*y 即348,84z z x y x y x y ??=-=-+??再求解方程，求得各驻点的坐标。一般方程组的符号解用solve 命令，当方程组不存在符号解时，solve 将给出数值解。求解方程的MA TLAB 代码为：

卷积神经网络CNN代码解析-matlab

卷积神经网络CNN代码解析 deepLearnToolbox-master是一个深度学习matlab包，里面含有很多机器学习算法，如卷积神经网络CNN，深度信念网络DBN，自动编码AutoEncoder（堆栈SAE，卷积CAE）的作者是Rasmus Berg Palm (rasmusbergpalm@https://www.wendangku.net/doc/0617214987.html,) 代码下载：https://https://www.wendangku.net/doc/0617214987.html,/rasmusbergpalm/DeepLearnToolbox 这里我们介绍deepLearnToolbox-master中的CNN部分。 DeepLearnToolbox-master中CNN内的函数： 调用关系为： 该模型使用了mnist的数字mnist_uint8.mat作为训练样本，作为cnn的一个使用样例， 每个样本特征为一个28*28=的向量。

网络结构为： 让我们来看看各个函数： 一、Test_example_CNN: (2) 三、cnntrain.m (5) 四、cnnff.m (6) 五、cnnbp.m (7) 五、cnnapplygrads.m (10) 六、cnntest.m (11) 一、Test_example_CNN: Test_example_CNN: 1设置CNN的基本参数规格，如卷积、降采样层的数量，卷积核的大小、降采样的降幅 2 cnnsetup函数初始化卷积核、偏置等 3 cnntrain函数训练cnn，把训练数据分成batch，然后调用 3.1 cnnff 完成训练的前向过程，

3.2 cnnbp计算并传递神经网络的error，并计算梯度（权重的修改量） 3.3 cnnapplygrads 把计算出来的梯度加到原始模型上去 4 cnntest 函数，测试当前模型的准确率 该模型采用的数据为mnist_uint8.mat， 含有70000个手写数字样本其中60000作为训练样本，10000作为测试样本。把数据转成相应的格式，并归一化。 设置网络结构及训练参数 初始化网络，对数据进行批训练，验证模型准确率 绘制均方误差曲线 二、Cnnsetup.m 该函数你用于初始化CNN的参数。 设置各层的mapsize大小， 初始化卷积层的卷积核、bias 尾部单层感知机的参数设置 * bias统一设置为0