Matlab信号与图像处理实验-模板.docVIP

本 科 实 验 报 告 实验名称： 小组名称： 长： 学 号： 组 员： 学 号： 组 员： 学 号： 组 员： 学 号： 组 员： 学 号： 专 业： 所属学院： 教 师： 实验日期： 2013-10-15 大学制 确定卷积结果的分段时限； 确定每段中积分的上下限； 确定每段中积分函数的表达式具体实现如下图 实验内容 对于信号 计算其卷积 实验分析 （1）根据题目要求输入信号，观察信号和卷积结果的波形表示。 信号f的波形为 信号g的波形为 卷积结果的波形为 由图形可知，卷积后结果信号的长度为两输入信号长度之和；且卷积结果（经db转换）每一点都是两输入信号翻转平移后重叠区域的面积。 说明卷积函数编写正确，卷积结果波形符合理论。 学习心得 通过编写MATLAB程序计算卷积，使我对信号之间计算卷积的方法和步骤有了更深的理解。同时这次实验使我对MATLAB这一数学软件有了进一步的理解，它除了功能强大且语言简单易懂外，尤其适合对信号的处理和计算。此次实验告诉我们在动手计算解决问题的同时要学会运用功能强大的计算机软件来解决问题并做仿真验证以进一步研究。 程序源码 主程序：%UNTITLED Summary of this function goes here % Detailed explanation goes here clear all % Parameters definition Tp=10.e-6; % transmitted pulse width fc=1.e+7; % carrier frequency Fs=200.e+6; % sampling frequency kr=2.e+12; % chirp rate tr=0:1/Fs:Tp; % discrete time array %Functions definition s=exp(j*2*pi*fc*tr+j*pi*kr*(tr).^2); f=exp(j*2*pi*fc*(tr-Tp)+j*pi*kr*(tr-Tp).^2)+exp(j*2*pi*fc*(tr-1.5*Tp)+j*pi*kr*(tr-1.5*Tp).^2); g=conj(s); % Processing y=conv_myself(f,g); y=y/max(y); y=20*log10(abs(y)); % plot result figure(1); plot(tr,f); title(function f); figure(2); plot(tr,g); title(function g); tc=0:1/Fs:2*Tp; figure(3); plot(tc,y);title(the result of conv); grid; 卷积程序 function w=conv_myself(f,g) len=length(f)+length(g)-1; fn=zeros(1,len); gn=zeros(1,len); fn(length(g):end)=f; for m=1:length(g) g1(m)=g(length(g)-m+1); end gn(1:length(g))=g1; for n=1:len gn1=circshift(gn,[0,n-1]); w(n)=sum(fn.*gn1); end