数字信号处理课程设计实验报告
基础实验
实验一离散时间系统及离散卷积
一、实验目的
（1）熟悉MATLAB软件的使用方法。
（2）熟悉系统函数的零极点分布、单位脉冲响应和系统频率响应等概念。（3）利用MATLAB绘制系统函数的零极点分布图、系统频率响应和单位脉冲响应。
（4）熟悉离散卷积的概念，并利用MATLAB计算离散卷积。
二、实验容
1、离散时间系统的单位脉冲响应
（1）选择一个离散时间系统；
（2）用笔进行差分方程的递推计算；
（3）编制差分方程的递推计算程序；
（4）在计算机上实现递推运算；
（5）将程序计算结果与笔算的计算结果进行比较，验证程序运行的正确性；2.离散系统的幅频、相频的分析方法
（1）给定一个系统的差分方程或单位取样响应；
（2）用笔计算几个特殊的幅频、相频的值，画出示意曲线图；
（3）编制离散系统的幅频、相频的分析程序；
（4）在计算机上进行离散系统的幅频、相频特性，并画出曲线；
（5）通过比较，验证程序的正确性；
3. 离散卷积的计算
（1）选择两个有限长序列，用笔计算其线性卷积；
（2）编制有限长序列线性卷积程序；
（3）利用计算程序对（1）选择的有限长序列进行卷积运算；
（4）比较结果验证程序的正确性。
三、实验要求
（1）编制实验程序，并给编制程序加注释；
（2）按照实验容项要求完成笔算结果；
（3）验证编制程序的正确性，记录实验结果。
（4）至少要求一个除参考实例以外的实例，在实验报告中，要描述清楚实例中
的系统，并对实验结果进行解释说明。
四、实验程序及其结果
1、离散时间系统的脉冲响应
% y(n)-0.9y(n-1)+ 0.5y(n-2)=x(n)+0.5x(n-1) 设x(n)=2^n; 0<=n<=15 clear all;
b=[1,0.5];a=[1,-0.9,0.5];
n=0:15; x=impseq(0,0,15)
h=filter(b,a,x); % 系统冲击响应
subplot(2,1,1);stem(n,h);title('系统冲激响应');
xlabel('n');ylabel('h');
n=0:15; x=2.*n; nx=0:15;nh=0:15;
y=conv_m(x,nx,h,nh);ny=length(y);
n=0:ny-1;
subplot(2,1,2);stem(n,y);title('系统对x(n)响应');
xlabel('n');ylabel('y');
B=roots(b);A=roots(a);
figure; zplane(B,A);
2、离散系统的幅频、相频的分析方法
% 差分方程为y(n)-1.76y(n-1)+1.1829y(n-2)-0.2781y(n-3)
% =0.0181x(n)+0.0543x(n-1)+0.0543x(n-2)+0.0181x(n-3) b=[0.0181,0.0543,0.0543,0.0181];
a=[1.000,-1.76,1.1829,-0.2781];
m=0:length(b)-1;
l=0:length(a)-1;
K=500;
k=0:1:K;
w=pi*k/K;
H=(b*exp(-j*m'*w))./(a*exp(-j*l'*w)); %计算频率相应
magH=abs(H); % magH为幅度
angH=angle(H); %a ngH为相位
subplot(2,1,1); plot(w/pi,magH);grid;
xlabel('以\pi为单位的频率'); ylabel('幅度');
title('幅度响应');
subplot(2,1,2); plot(w/pi,angH); grid;
xlabel('以\pi为单位的频率'); ylabel('相位');
title('相位响应');
3、离散卷积的计算
% x=[1,4,3,5,3,6,5] , -4<=n<=2
% h=[3,2,4,1,5,3], -2<=n<=3
% 求两序列的卷积
clear all;
x=[1,4,3,5,3,6,5]; nx=-4:2;
h=[3,2,4,1,5,3];nh=-2:3;
ny=(nx(1)+nh(1)):(nx(length(x))+nh(length(h)));
y=conv(x,h);
n=length(ny);
subplot(3,1,1);stem(nx,x);xlabel('nx');ylabel('x');
subplot(3,1,2);stem(nh,h);xlabel('nh');ylabel('h');
subplot(3,1,3);stem(ny,y);xlabel('n');ylabel('x和h的卷积');
实验二 离散傅里叶变换与快速傅里叶变换
一、实验目的
（1）加深理解离散傅立叶变换及快速傅立叶变换概念； （2）学会应用FFT 对典型信号进行频谱分析的方法； （3）研究如何利用FFT 程序分析确定性时间连续信号； （4）熟悉应用FFT 实现两个序列的线性卷积的方法。 二、实验容
1、用离散傅立叶变换程序处理时间抽样信号，并根据实序列离散傅立叶变换的对称性，初步判定程序的正确性；
2.观察三角波和反三角波序列的时域和幅频特性，用N=8点FFT 分析信号序列的幅频特性，观察两者的序列形状和频谱曲线有什么异同？绘出两序列及其幅频特性曲线。
三角波序列
反三角波序列
??
?
??≤-≤≤+=else
n n
n n n x c 048301)(??
?
??≤-≤≤-=else
n n n n n x d 043
304)(
3.已知余弦信号如下
当信号频率 F=50Hz 采样间隔 T=0.000625s ，采样长度 N=64时，对该信号进行傅里叶变换。用FFT 程序分析正弦信号，分别在以下情况进行分析。
（1）F=50, N=32,T=0.000625； （2）F=50, N=32,T=0.005； （3）F=50, N=32,T=0.0046875； （4）F=50, N=32,T=0.004； （5）F=50, N=64=0.000625
4.选定某一时间信号进行N=64点离散傅立叶变换，并且，对同一信号进行快速傅立叶变换，并比较它们的速度。 三、实验要求
1、编制DFT 程序及FFT 程序，并比较DFT 程序与FFT 程序的运行时间。给编制的程序加注释；
2.完成实验容2，并对结果进行分析。在单位圆Z 上的变换频谱会相同吗？如果不同，你能说出那个低频分量更多吗？为什么？
3.完成实验容3，并对结果进行分析；
4.利用编制的计算卷积的计算程序，给出一下三组函数的卷积结果。 （1）
)
2cos()(Ft t x π=??
?≤≤=其它
14
01)(n n x 14
054)(≤≤?
?
? ??=n n h n
（2） （3）
四、实验程序及其结果分析 1、离散傅里叶变换
function [ Xk] = dft( xn,N ) %计算N 点的DFT 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; End
2、（1）三角波 clear all; n=0:1:7;
x1=(n+1).*(n>=0).*(n<=3)+(8-n).*(n>=4); % 三角波序列 subplot(3,1,1);stem(n,x1); title('三角波序列');xlabel('n'); N=8; k=0:N-1; X1=fft(x1,N);
magX1=abs(X1);phaX1=angle(X1); subplot(3,1,2);stem(k,magX1); xlabel('k');ylabel('三角波DFT 的幅度'); subplot(3,1,3);stem(k,phaX1);
xlabel('k');ylabel('三角波DFT 的相位');
??
?≤≤=其它
9
01)(n n x 19
0)5.0sin(5.0)(≤≤=n n n h 9
0)
1.01()(≤≤-=n n n x 9
01.0)(≤≤=n n
n h
（2）反三角波
clear all;
n=0:1:7;
x2=(4-n).*(n>=0).*(n<=3)+(n-3).*(n>=4); % 反三角波序列subplot(3,1,1);stem(n,x2);
title('反三角波序列');xlabel('n');
N=8; k=0:N-1;
X2=fft(x2,N);
magX2=abs(X2);phaX2=angle(X2);
subplot(3,1,2);stem(k,magX2);
xlabel('k');ylabel('幅度');
subplot(3,1,3);stem(k,phaX2);
xlabel('k');ylabel('相位');
结果分析：由图知：三角波和反三角波序列的波形不同，当N=8时，正反三角波的幅频特性相同，因为两者的时域只差一个相位
3、余弦信号的FFT分析
clea all;
F=50; N=32;T=0.000625; % (1) F=50,N=32,T=0.000625
n=1:N;
x=cos(2*pi*F*n*T);
figure(1); subplot(2,1,1); plot(n,x);
ylabel('x(n)'); xlabel('n');title('(1) F=50,N=32,T=0.000625');
X=fft(x); magX=abs(X);
subplot(2,1,2);plot(n,X);
ylabel('FFT|X|'); xlabel('f(pi)');
F=50; N=32;T=0.005; % (2) F=50,N=32,T=0.005
n=1:N;
x=cos(2*pi*F*n*T);
figure(2);
subplot(2,1,1);plot(n,x);ylabel('x(n)');
xlabel('n');title('(2) F=50,N=32,T=0.005');
X=fft(x); magX=abs(X);
subplot(2,1,2); plot(n,X); ylabel('FFT|X|'); xlabel('f(pi)');
F=50,N=32,T=0.0046875; % (3) F=50,N=32,T=0.0046875 n=1:N;
x=cos(2*pi*F*n*T);
figure(3);
subplot(2,1,1); plot(n,x); ylabel('x(n)');
xlabel('n');title('(3) F=50,N=32,T=0.0046875');
X=fft(x);
magX=abs(X);
subplot(2,1,2);
plot(n,X); ylabel('FFT|X|'); xlabel('f(pi)');
F=50,N=32,T=0.004; % (4) F=50,N=32,T=0.004
n=1:N; x=cos(2*pi*F*n*T);
figure(4);
subplot(2,1,1);
plot(n,x); ylabel('x(n)');
xlabel('n');title('(4) F=50,N=32,T=0.004');
X=fft(x); magX=abs(X);
subplot(2,1,2);
plot(n,X);ylabel('FFT|X|'); xlabel('f(pi)');
F=50,N=64,T=0.000625; % (5) F=50,N=64,T=0.000625
n=1:N;
x=cos(2*pi*F*n*T);
figure(5); subplot(2,1,1); plot(n,x); ylabel('x(n)'); xlabel('n'); title('(5) F=50,N=64,T=0.000625');
X=fft(x); magX=abs(X);
subplot(2,1,2); plot(n,X); ylabel('FFT|X|'); xlabel('f(pi)');
结果分析；不同的采样间隔会产生不同的栅栏效应，相当于透过栅栏观赏风景,只能看到频谱的一部分,而其它频率点看不见,因此很可能使一部分有用的频率成分被漏掉，从而产生不同的频谱图，减小栅栏效应可用提高采样间隔也就是频率分辨力的方法来解决。间隔小，频率分辨力高，被“挡住”或丢失的频率成分就会越少。但会增加采样点数，使计算工作量增加。
4、卷积计算
clear all;
n=0:14;x1=1.*(n>=0).*(n<=14); % 第一组函数
h1=(4/5).^n;
N=27; y=circonvt(x1,h1,N);
disp(y);subplot(3,1,1);
i=0:N-1; stem(i,y); ylabel('式（1）卷积');
n=0:9;x2=1.*(n>=0).*(n<=9); % 第二组函数
n=0:19;h2=(0.8).*sin(0.5.*n);
N=27; y=circonvt(x2,h2,N);
disp(y); subplot(3,1,2);
i=0:N-1; stem(i,y); ylabel('式（2）卷积');
n=0:9; x3=(1-0.1.*n); h3=0.1.*n; % 第三组函数
N=17; y=circonvt(x3,h3,N);
disp(y); subplot(3,1,3);
i=0:N-1; stem(i,y); ylabel('式（3）卷积');
实验三IIR数字滤波器设计
一、试验目的
1、学习模拟－数字变换滤波器的设计方法；
2.掌握双线性变换数字滤波器设计方法；
3.掌握实现数字滤波器的具体方法。
二、实验容
1、设计一个巴特沃思数字低通滤波器，设计指标如下：
通带w 阻带w>ws=0.35pi 幅度衰减不小于15dB； 2、编制计算设计的数字滤波器幅度特性和相位特性的程序，并进行实验验证。 3、编制实现该数字滤波器程序并且实现数字滤波 （1）分别让满足所设计的滤波器的通带、过渡带、阻带频率特性的正弦波通过滤波器，验证滤波器性能； （2）改变正弦抽样时间，验证数字低通滤波器的模拟截止频率实抽样时间的函数。 三、实验要求 1、编制实验容要求的程序，并给程序加注释； 2、根据实验结果，给出自己设计的数字滤波器的幅度特性和相位特性； 3.用所设计的滤波器对不同频率的正弦波信号进 行滤波，以说明其特性； 4、fp=0.2KHz,Rp=1dB,fs=0.3KHz,As=25dB,T=1ms；分别用脉冲响应不变法及双线性变换法设计一Butterworth数字低通滤波器，观察所设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。 四、实验程序及结果分析 1、巴特沃斯低通滤波器 clear all; wp=0.2*pi;ws=0.3*pi;Rp=1;As=15; R=10^(-Rp/20);A=10^(-As/20); [b,a]=afd_buttap(wp,ws,Rp,As); [ C,B,A] = sdir2cas(b,a); %将直接型装换成级联型 [db,mag,pha,w]=freqs_m(b,a,0.5*pi); %计算系统频率响应 [ha,x,t]=impulse(b,a); subplot(2,1,1);plot(w/pi,mag);title('幅度响应');grid; ylabel('幅度');xlabel('以\pi为单位的频率'); subplot(2,1,2);plot(w/pi,pha);title('相位响应');grid; ylabel('相位');xlabel('以\pi为单位的频率'); 2、脉冲响应不变法设计巴特沃斯滤波器 % 冲击响应不变法设计数字巴特沃斯低通滤波器 clear all; T=0.0001; wp=200*2*pi*T;Rp=1; ws=300*2*pi*T;As=25; %数字滤波器指标T=1; omegap=wp/T; omegas=ws/T; %转换为模拟域指标 [cs,ds]=afd_buttap(omegap,omegas,Rp,As); %模拟巴特沃斯滤波器的计算% 冲击响应不变法 [b,a]=imp_invr(cs,ds,T); [C,B,A]=dir2par(b,a); % 将IIR滤波器直接型转换为并联型的函数 % 计算数字滤波器频率响应 [db,mag,pha,grd,w]=freqz_m(b,a); subplot(3,1,1);plot(w/pi,mag); title('冲击响应不变法设计巴特沃斯滤波器'); axis([0 0.5 0 1]);grid; ylabel('|幅度|');xlabel('以\pi为单位的频率'); subplot(3,1,2);plot(w/pi,db);axis([0 0.8 -400 0]);grid; ylabel('对数幅度（db）');xlabel('以\pi为单位的频率'); 3、双线性变换法设计巴特沃斯滤波器 clear all; T=0.0001;Fs=1/T; wp=200*2*pi*T;Rp=1; ws=300*2*pi*T;As=25; %数字滤波器指标 % 转换成模拟域指标 omegap=(2/T)*tan(wp/2); omegas=(2/T)*tan(ws/2); %模拟巴特沃斯滤波器的计算 [cs,ds]=afd_buttap(omegap,omegas,Rp,As); %双线性变换 [b,a]=bilinear(cs,ds,Fs); [C,B,A]=sdir2cas(b,a); [db,mag,pha,grd,w]=freqz_m(b,a); subplot(3,1,1);plot(w/pi,mag);ylabel('|幅度|');title('双线性变换法设计巴特沃斯滤波器'); xlabel('以\pi为单位的频率');axis([0 0.5 0 1]);grid; subplot(3,1,2);plot(w/pi,db);axis([0 0.8 -400 0]);grid; ylabel('对数幅度（db）');xlabel('以\pi为单位的频率'); subplot(3,1,3);plot(w/pi,pha);axis([0 0.4 -4 4]);grid; ylabel('相位');xlabel('以\pi为单位的频率'); 4、数字滤波 clear all; F=0.1;N1=512;T1=0.1; n=0:5:N1; x=sin(2*pi*F*n*T1); wp=0.2*pi;ws=0.35*pi; Rp=1;As=15; T=1;Fs=1/T; OmegaP=(2/T)*tan(wp/2); %预畸变 OmegaS=(2/T)*tan(ws/2); ep=sqrt(10^(Rp/10)-1); %求平方根 Ripple=sqrt(1/(1+ep.^2));Attn=1/10^(As/20); N=ceil((log10((10^(Rp/10)-1)/(10^(As/10)-1)))/(2*log10(OmegaP/OmegaS))); OmegaC=OmegaP/((10.^(Rp/10)-1).^(1/(2*N))); [cs,ds]=u_buttap(N,OmegaC); [b,a]=bilinear(cs,ds,Fs); %括号为模拟滤波器系数，中括号中的为数字滤波器系数 [mag,db,pha,w]=freqz_m(b,a); y=filter(b,a,x); subplot(2,1,1);stem(x); axis([0 100 -1 1]);title('输入正弦序列'); subplot(2,1,2);stem(y); axis([0 100 -1 1]); title('滤波输出序列'); 结果分析：滤波输出序列前有一段是由于延时造成的。 实验四FIR数字滤波器 1、FIR滤波器的设计 clear all; % 用哈明窗函数法来设计FIR滤波器选N=51 wp=0.2*pi; ws=0.35*pi; N=51; n=[0:1:N-1]; wc=(ws+wp)/2; % 理想低通的截止频率 hd=ideal_lp(wc,N); % 理想低通的冲激响应 w_ham=(hamming(N))'; h=hd.*w_ham % FIR滤波器的冲激响应 [db,mag,pha,grd,w]=freqz_m(h,[1]); delta_w=2*pi/1000; Rp=-(min(db(1:1:wp/delta_w+1))); % 实际的通带衰减 As=-round(max(db(ws/delta_w+1:1:501))); % 实际的阻带衰减 subplot(2,2,1);stem(n,hd);title('理想冲激响应'); axis([0 N-1 -0.1 0.3]);xlabel('n'); ylabel('hd(n)'); subplot(2,2,2);stem(n,h);title('实际冲激响应'); axis([0 N-1 -0.1 0.3]);xlabel('n'); ylabel('h(n)'); subplot(2,2,3);plot(w/pi,pha);title('相位响应');grid; xlabel('以\pi为单位的频率'); ylabel('相位'); subplot(2,2,4);plot(w/pi,db);title('幅度响应（db）');grid; xlabel('以\pi为单位的频率'); ylabel('H(db)'); 2、分别用矩形窗、汉宁窗、海明窗、三角窗和Blackman窗设计一个线性相位 的带通滤波器。 clear all; w1=0.3*pi; w2=0.5*pi; N=51; n=[0:1:N-1]; hd=ideal_lp(w2,N)-ideal_lp(w1,N); % 理想低通的冲激响应 w_rec=(rectwin(N))'; h1=hd.*w_rec ; % 矩形窗FIR滤波器的冲激响应 [db,mag,pha,grd,w]=freqz_m(h1,[1]); figure;plot(w/pi,db);title('矩形窗幅度响应（db）'); xlabel('以\pi为单位的频率'); ylabel('H1(db)'); w_hann=(hann(N))'; h2=hd.*w_hann % 汉宁窗FIR滤波器的冲激响应 [db,mag,pha,grd,w]=freqz_m(h2,[1]); figure;plot(w/pi,db);title('汉宁窗幅度响应（db）'); xlabel('以\pi为单位的频率'); ylabel('H2(db)'); w_ham=(hamming(N))'; h3=hd.*w_ham % 哈明窗FIR滤波器的冲激响应 吉林大学本科“十二五”规划教材修订教材选题立项名单 序号教材名称版次主编主编单位方式出版单位 1 田野考古学一版冯恩学文学院考古系修订吉林大学出版社 2 藏品管理学二版吕军文学院博物馆系修订高等教育出版社 3 语言国情学一版刘佐艳外国语学院修订吉林大学出版社 4 新世纪实用英语写作三版王青华公共外语教育学院修订外语教研出版社 5 国际贸易一版李俊江经济学院修订高等教育出版社 6 发展经济学三版马春文经济学院修订高等教育出版社 7 保险学教程一版池晶经济学院修订科学出版社 8 刑法学（上、下册）一版李洁法学院修订中人大出版社 9 地缘政治学一版刘雪莲行政学院修订高等教育出版社 10 当代国际关系一版黄凤志行政学院修订吉林大学出版社 11 会计信息系统原理与实验教程二版李清商学院修订清华大学出版社 12 电子档案管理基础一版王萍管理学院修订清华大学出版社 13 大学数学系列教材三版李辉来数学中心: 修订高等教育出版社 14 经济数学基础系列教材二版孙毅数学中心: 修订清华大学出版社 15 工程数学系列教材二版袁洪军数学中心: 修订高等教育出版社 16 常微分方程一版伍卓群、李勇、 史少云 数学学院修订高等教育出版社 17 高等代数一版杜现昆、原永久、 牛凤文 数学学院修订高等教育出版社 18 近代物理实验（一）二版韩炜物理学院修订科学出版社 19 热力学与统计物理学一版崔海宁物理学院修订吉大出版社 谢谢你的阅读 20 医用物理学二版梁路光物理教学中心修订高教出版社 21 大学物理学（上下）一版张铁强物理教学中心修订高教出版社 22 医用物理实验一版付研物理教学中心修订高教出版社 23 无机化学（下）二版宋天佑化学学院修订高等教育出版社 24 有机化学实验二版丁长江化学学院修订科学出版社 25 精细化工概要一版高大维、张玉敏化学学院修订吉林大学出版社 26 综合和设计化学一版张寒琦化学学院修订高等教育出版社 27 生命科学仪器使用技术教程一版滕利荣生命科学学院修订科学出版社 28 电视原理一版李秀英电子科学与工程学院修订西安电子科技大学出版社 29 计算机操作系统教程三版左万历计算机科学与技术学院修订高等教育出版社 30 大学计算机程序设计基础一版陈娟计算机科学与技术学院修订清华大学出版社 31 数字电子技术二版杨永健软件学院修订人民邮电出版社 32 计算机网络三版刘衍珩软件学院修订科学出版社 33 机械制造技术基础二版邹青机械科学与工程学院修订机械工业出版社 34 机械设计二版谭庆昌机械科学与工程学院修订高等教育出版社 35 工程设计制图一版潘淑璋机械科学与工程学院修订科学出版社 36 工程流体力学一版于萍机械科学与工程学院修订科学出版社 37 工程机械设计一版秦四成机械科学与工程学院修订科学出版社 38 CATIA V5R21实用基础教程一版潘志刚机械科学与工程学院修订科学出版社 39 AutoCAD2013实用教程一版张云辉机械科学与工程学院修订科学出版社 40 现代工程制图一版张秀芝机械科学与工程学院修订科学出版社 41 机械设计课程设计一版寇尊权机械科学与工程学院修订机械工业出版社 42 机械设计基础一版贾艳辉机械科学与工程学院修订高等教育出版社 43 机械精度设计与检测二版陈晓华、侯磊机械科学与工程学院修订中国质检出版社 44 汽车设计四版宋传学汽车工程学院修订机械工业出版社 谢谢你的阅读 上海杉达学院 商务流程综合实训总结 单位名称：微科电子有限公司（加盖公章）姓名：陈恩娜 学院：胜祥商学院 专业：国际经济与贸易 班级： f130219 时间： 2016.11.14 工作总结历时10天的跨专业实训圆满落幕，作为国贸专业的我们参加了企业经营的模拟实训，通过这几天的实训，我也算是了解到了企业经营的基本流程和一般模式。前三天，我们基本在忙碌企业创立的事情。我们首先要做的是公司注册登记。公司注册流程共包括名称预先登记、设立登记申请书、准备申请材料、银行开户入资、验资、前置审批、报送申请材料、工商审批发照、刻制公章、开设银行帐户、办理各项登记、股东入资证明、企业机构代码、企业税务登记、企业劳动备案、社会保险登记、企业统计备案、特殊行业备案18类。我们实训时规定了公司类型为制造型企业，开始注册资金为500万元。在实训中，我们注册的公司为微科电子有限公司，地址位于上海市陆家嘴，股东为3人。申请表提交上去后，我们得到了审批，就开始了下面的企业经营规划。我们企业经营过程共分为九个相互联系又相互独立的部门。这十个部门分别为总经办、人力资源部门、物流部门、质检部门、行政部门、营销部门、生产部门、财务部门、采购部门。总经办为CEO代理，也就是我们常说的总经理，其他分别为营销总监、生产总监、物流总监、采购总监、财务总监、采购经理、人事经理、行政主管、质检经理。这九个职位分别为十一个同学完成，我作为采购部经理也参与其中。采购部的岗位职责可归纳为：依照公司生产需要及物资采购计划，全面负责公司的物料采购和供应工作；对初次进行合作的供应商进行调查了解，根据实际掌握的资料与信息做出相应的评价判定。对于符合公司要求的，方可与之开展业务往来与合作；采购工作的开展应当做到精打细算，尽力降低公司的采购成本，避免出现浪费公司资金的情况；负责对所有与公司有业务往来的供应商或供货企业进行定期的资质评价并给出明确的评价等级。针对不同的评价等级给出相应的处理意见；完成公司领导布置的其它各项工作。各岗位要各司其职，填写经营流程表，有序的完成一年的经营。 数字信号处理作业提交日期：2016年7月15日 实验一 维纳滤波器的设计 第一部分 设计一维纳滤波器。 (1)产生三组观测数据，首先根据()(1)()s n as n w n =-+产生信号()s n ，将其加噪(信噪比分别为20,10,6dB dB dB )，得到观测数据123(),(),()x n x n x n 。 (2)估计()i x n ，1,2,3i =的AR 模型参数。假设信号长度为L ，AR 模型阶数为N ，分析实验结果，并讨论改变L ，N 对实验结果的影响。 1 实验原理 滤波技术是信号分析、处理技术的重要分支，无论是信号的获取、传输，还是信号的处理和交换都离不开滤波技术，它对信号安全可靠和有效灵活地传递是至关重要的。信号分析检测与处理的一个十分重要的内容就是从噪声中提取信号，实现这种功能的有效手段之一是设计一种具有最佳线性过滤特性的滤波器，当伴有噪声的信号通过这种滤波器的时候，它可以将信号尽可能精确地重现或对信号做出尽可能精确的估计，而对所伴随噪声进行最大限度地抑制。维纳滤波器就是这种滤波器的典型代表之一。 维纳（Wiener ）是用来解决从噪声中提取信号的一种过滤（或滤波）方法。这种线性滤波问题，可以看做是一种估计问题或一种线性估计问题。 设一线性系统的单位样本响应为()h n ，当输入以随机信号()x n ，且 ()() () x n s n v n =+，其中()s n 表示原始信号，即期望信号。()v n 表示噪声，则输出()y n 为()=()()m y n h m x n m -∑，我们希望信号()x n 经过线性系统()h n 后得到的()y n 尽可能接近 于()s n ，因此称()y n 为估计值，用?()s n 表示。 则维纳滤波器的输入-输出关系可用下面表示。 设误差信号为()e n ,则?()()()e n s n s n =-,显然)(n e 可能是正值，也可能是负值，并且它是一个随机变量。因此，用它的均方误差来表达误差是合理的，所谓均方误差最小即 它的平方的统计期望最小：222?[|()|][|()()|][|()()|]E e n E s n s n E s n y n =-=-=min 。而要使均方误差最小，则需要满足2[|()|]j E e n h ?=0. 进一步导出维纳-霍夫方程为：()()()()*(),0,1,2...xs xx xx i R m h i R m i R m h m m =-==∑ 写成矩阵形式为：xs xx R R h =，可知：1xs xx h R R -=。表明已知期望信号与观测数据的互相关函数以及观测信号的自相关函数时，可以通过矩阵求逆运算，得到维纳滤波器的 DSP 课程设计实验 一、语音信号的频谱分析： 要求首先画出语音信号的时域波形，然后对语音信号进行频谱分析。在MATLAB 中，可以利用函数fft 对信号进行快速傅立叶变换，得到信号的频谱特性，从而加深对频谱特性的理解。 其程序为： >> [y,fs,bits]=wavread('I:\',[1024 5120]); >> sound(y,fs,bits); >> Y=fft(y,4096); >> subplot(221);plot(y);title('原始信号波形'); | >> subplot(212);plot(abs(Y));title('原始信号频谱'); 程序运行结果为： 二、设计数字滤波器和画出频率响应： 根据语音信号的特点给出有关滤波器的性能指标： 低通滤波器性能指标，p f =1000Hz ，c f =1200Hz ，s A =100dB ，p A =1dB ； 高通滤波器性能指标，c f =4800Hz ，p f =5000Hz ，s A =100dB ，p A =1dB ； 带通滤波器性能指标，1p f =1200Hz ，2p f =3000Hz ，1c f =1000Hz ，2c f =3200Hz ，s A =100dB ， p A =1dB ； 】 要求学生首先用窗函数法设计上面要求的三种滤波器，在MATLAB中，可以利用函数firl 设计FIR滤波器；然后再用双线性变换法设计上面要求的三种滤波器，在MATLAB中，可以利用函数butte、cheby1和ellip设计IIR滤波器；最后，利用MATLAB中的函数freqz画出各种滤波器的频率响应，这里以低通滤波器为例来说明设计过程。 低通： 用窗函数法设计的低通滤波器的程序如下： >> fp=1000;fc=1200;As=100;Ap=1;fs=22050; >> wc=2*fc/fs;wp=2*fp/fs; >> N=ceil(/*(wc-wp)/2))+1; >> beta=*; >> Win=Kaiser(N+1,beta); 、 >>b=firl(N,wc,Win); >>freqz(b,1,512,fs); 程序运行结果： 这里选用凯泽窗设计，滤波器的幅度和相位响应满足设计指标，但滤波器长度（N=708）太长，实现起来很困难，主要原因是滤波器指标太苛刻，因此，一般不用窗函数法设计这种类型的滤波器。 用双线性变换法设计的低通滤波器的程序如下： >> fp=1000;fc=1200;As=100;Ap=1;fs=22050; >> wc=2*fc/fs;wp=2*fp/fs; 》 >> [n,wn]=ellipord(wp,wc,Ap,As); >> [b,a]=ellip(n,Ap,As,wn); >> freqz(b,a,512,fs); ^ 实验一MATLAB语言的基本使用方法 实验类别：基础性实验 实验目的： （1）了解MATLAB程序设计语言的基本方法，熟悉MATLAB软件运行环境。 （2）掌握创建、保存、打开m文件的方法，掌握设置文件路径的方法。 （3）掌握变量、函数等有关概念，具备初步的将一般数学问题转化为对应计算机模型并进行处理的能力。 （4）掌握二维平面图形的绘制方法，能够使用这些方法进行常用的数据可视化处理。 实验内容和步骤： 1、打开MATLAB，熟悉MATLAB环境。 2、在命令窗口中分别产生3*3全零矩阵，单位矩阵，全1矩阵。 3、学习m文件的建立、保存、打开、运行方法。 4、设有一模拟信号f(t)=1.5sin60πt,取?t=0.001,n=0,1,2,…,N-1进行抽样，得到 序列f(n),编写一个m文件sy1_1.m,分别用stem,plot,subplot等命令绘制32 点序列f(n)(N=32)的图形，给图形加入标注，图注，图例。 5、学习如何利用MATLAB帮助信息。 实验结果及分析： 1）全零矩阵 >> A=zeros(3,3) A = 0 0 0 0 0 0 0 0 0 2）单位矩阵 >> B=eye(3) B = 1 0 0 0 1 0 0 0 1 3）全1矩阵 >> C=ones(3) C = 1 1 1 1 1 1 1 1 1 4)sy1_1.m N=32; n=0:N-1; dt=0.001; t=n*dt; y=1.5*sin(60*pi*t); subplot(2,1,1), plot(t,y); xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('正弦函数'); title('二维图形'); subplot(2,1,2), stem(t,y) xlabel('t'); ylabel('y=1.5*sin(60*pi*t)'); legend('序列函数'); title('条状图形'); 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 二维图形 00.0050.010.0150.020.0250.030.035 t y = 1 . 5 * s i n ( 6 * p i * t ) 条状图形 DSP课程设计总结（2013-2014学年第2学期） 题目： 专业班级：电子1103 学生姓名：万蒙 学号：11052304 指导教师： 设计成绩： 2014 年6 月 目录 一设计目的----------------------------------------------------------------------3 二系统分析----------------------------------------------------------------------3 三硬件设计 3.1 硬件总体结构-----------------------------------------------------------3 3.2 DSP模块设计-----------------------------------------------------------4 3.3 电源模块设计----------------------------------------------------------4 3.4 时钟模块设计----------------------------------------------------------5 3.5 存储器模块设计--------------------------------------------------------6 3.6 复位模块设计----------------------------------------------------------6 3.7 JTAG模块设计--------------------------------------------------------7 四软件设计 4.1 软件总体流程-----------------------------------------------------7 4.2 核心模块及实现代码---------------------------------------8 五课程设计总结-----------------------------------------------------14 机械工程及自动化专业本科培养方案2009版 一、培养目标 培养适应现代化建设和未来社会与科技发展需要，立志为国家富强、民族振兴和人类文明进步而奋斗，德智体美全面发展与健康个性和谐统一，富有创新精神和实践能力的高级机械工程专业人才。 学生毕业后可在有关企业、科研单位、国家机关和高等院校从事设计、制造、管理、营销、科研和教学等工作。 二、业务培养要求 本专业培养的学生具有坚实的数学、力学、外语和计算机基础,掌握机械科学的基础理论和基本技能,了解市场经济的基本知识,获得初步的科学研究、科技开发、组织管理和社会活动能力训练。本专业的毕业生应能在以下几方面获得相应的知识和能力： 1．掌握现代设计的基本理论和技能, 有较强的自动化技术应用能力，能够较熟练地运用计算机辅助设计技术，初步具有综合运用机械、电子、液压等知识进行机械产品设计的能力。 2．熟悉现代机械制造的基本理论、技术和装备，能够制定机械产品的加工和装配工艺规程,正确选择和设计工艺装备，具有加工质量及产品性能的检测、分析与控制的基本知识和能力。 3. 掌握自动控制的基本理论，熟悉数控技术，初步具有现代机械制造系统如数控机床、加工中心、柔性制造系统的运行和维护能力。 4.初步具有新工艺、新技术、新设备的研究与开发能力，懂得机械产品及其制造过程的技术经济分析与生产组织管理。 5．具有较强的自学能力，掌握独立获取、消化和应用新知识的能力和方法；基本掌握一门外国语，能顺利阅读本专业的外文资料，具有一定的国际学术交流能力。 6．具有一定的市场经济知识与管理知识；懂得一定的法律知识和国防知识。 三、主干学科及主要课程 主干学科：机械设计及理论、机械制造及自动化、机械电子工程。 主要课程：理论力学、材料力学、工程图学、工程材料、机械原理、机械设计、制造技术基础、电工学、控制工程基础、微机原理与接口技术、测试与传感技术、机电传动控制、液压与气压传动、机械制造装备设计、机电控制系统分析与设计、工程机械设计等。 四、专业特色及专业方向 机械工程及自动化专业是集机械、电子、信息技术为一体的综合性专业，知识结构先进，综合性强。培养的学生知识面宽，思路开阔，具有较强的创新意识和实践能力，适应面广阔，发展迅速，需求巨大。 五、学制 一般为四年 六、学位授予 工学学士 七、毕业合格标准 1．具有良好的思想和身体素质，符合学校规定的德育和体育标准。 2．通过培养方案规定的全部教学环节，按照规定内容取得212学分。 宁波职业技术学院 外贸实务 II - 提高 实训报告 指导老师：江彬 班级：国贸3141 学生姓名：肖思洁 学号： 1426263133 日期： 2016-04-25 课程名称：外贸实务II-提高 1 实训目的及要求 1．1 1．2 2 实训内容及步骤（包含简要的实训步骤流程） 2．1 本人所扮演的角色 2．2 贸易资料及实训步骤 2．2．1 2．2．2 2．2．3 2．2．4 3 实训结果（包括实训项目的完成情况，代表性邮件，单据，程序或图表、结论陈述、核算表数据记录及分析等） 3．1 3．1．1 3．1．2 3．1．3 3．2 4 实训中遇到的问题及其解决方法 5 实训总结（包括心得体会、对SimTrade系统的评价、对自己实训效果的评价如实训收获不足及实训改进意见等） 6 实训评价 外贸实训报告 外贸实务实训体会总结，这次是项目过程的完成让我感觉很深刻。短短的32学时已经结束，静下心来回想这次操作模拟学习真是感受颇深。我们知道动手操作是大学教育中一个极为重要的实践性环节，通过实习，可以使我们在实践中接触与本专业相关的一些实际工作，培养和锻炼我们综合运用所学的基础理论、基本技能和专业知识，去独立分析和解决实际问题的能力，把理论和实践结合起来，提高我们的实际动手能力，为将来我们毕业后走上工作岗位打下一定的基础。通过这段时间的学习，从无知到认知，到深入了解，渐渐地我喜欢上这个专业，让我深刻的体会到学习的过程是最美的，在整个实习过程中，我每天都有很多的新的体会，新的想法。 回顾我的实习生活，感触是很深的，收获也是丰硕的。可以模拟出口商与非洲、中东等地方的外国商人做生意，他们在我公司下订单，我们再把订单下到厂里，从中赚取差额，或者作为进口商与出口商讨价还价，去除运费保险费等来赚取利益。当然对工厂的基本流程也有一定的了解。这次的实习经历我对外贸这个专业有了更加理性的认识和更深刻的体会。在这次是学习中，我学到了很多过去两年没有体会到的东西，这不仅仅只是上课模拟，也是一次对真实工作流程经历。 对实践的看法： 在操作过程中，根据本门课程的内容、特点，通过走出去、请进来等方式，精心组织方案。通过听、看、做使一些看起来繁杂的专业知识很快被我们理解和掌握。只有将理论联系实际，教学与实际相结合，才是培养我们能力的一种有效形式。 出口商+进口商+工厂，一共写了78封邮件。一共发布了8封广告和信息。 查了B2B里的多数产品信息。 银行汇率：欧元大多汇率为8.6402，美元大多为6.1463等。保险费：一切险（ALL RISKS）为0.8%，战争险（W AR RISKS）为0.08% 保险费计算方式为 (1)按CIF进口时：保险金额＝CIF货价×1.1 (2)按CFR进口时：保险金额＝CFR货价×1.1 / (1 - 1.1 ×r)，其中r为保险费率，请在"淘金网"的"保险费"页面查找，将所投险别的保险费率相加即可。 (3)按FOB进口时：保险金额＝(FOB货价+ 海运费)×1.1 / (1 - 1.1 ×r)，其中FOB 货价就是合同金额，海运费请在装船通知中查找，由出口商根据配舱通知填写，如果出口商填写错误，请其查看配舱通知。 实训的基本流程： 第一周完成了进口商，出口商，工厂，进口、出口地银行的基本资料。然后熟悉了下系统的基本轮廓，如B2B里面可以查询写什么，市场，海关等在哪个位置。根据老师的知道，试着去发广告与写邮件。 第二周确定角色，开始寻找有利信息，搜索信息，同业务伙伴建合作关系。 我先进行成本、费用、利润等的核算，若有盈利则进一步磋商合作，若亏损就跟对方进行讨价还价。过程为询盘——发盘——还盘——接受。 第三周之后进入交易准备阶段——交易磋商阶段——签订合同(T/T+FOB)与接收信用证（L/C+CIF）——履行合同阶段。 签订合同之后进行合同履行阶段。 首先作为出口商，与进口商进行磋商商定后确定的价格，之后跟工厂进行合作并进一步签订SALES CONFIRMATON。等工厂交货物发过来后，与工厂的业务就能完成。 武汉工程大学 数字信号处理实验报告 姓名：周权 学号：1204140228 班级：通信工程02 一、实验设备 计算机，MATLAB语言环境。 二、实验基础理论 1.序列的相关概念 2.常见序列 3.序列的基本运算 4.离散傅里叶变换的相关概念 5.Z变换的相关概念 三、实验内容与步骤 1.离散时间信号（序列）的产生 利用MATLAB语言编程产生和绘制单位样值信号、单位阶跃序列、指数序列、正弦序列及随机离散信号的波形表示。 四实验目的 认识常用的各种信号，理解其数字表达式和波形表示，掌握在计算机中生成及绘制数字信号波形的方法，掌握序列的简单运算及计算机实现与作用，理解离散时间傅里叶变换，Z变换及它们的性质和信号的频域分 实验一离散时间信号（序列）的产生 代码一 单位样值 x=2; y=1; stem(x,y); title('单位样值 ') 单位阶跃序列 n0=0; n1=-10; n2=10; n=[n1:n2]; x=[(n-n0)>=0]; stem(n,x); xlabel('n'); ylabel('x{n}'); title('单位阶跃序列'); 实指数序列 n=[0:10]; x=(0.5).^n; stem(n,x); xlabel('n'); ylabel('x{n}'); title('实指数序列'); 正弦序列 n=[-100:100]; x=2*sin(0.05*pi*n); stem(n,x); xlabel('n'); ylabel('x{n}'); title('正弦序列'); 随机序列 n=[1:10]; x=rand(1,10); subplot(221); stem(n,x); xlabel('n'); ylabel('x{n}'); title('随机序列'); 电气信息工程学院 D S P技术与综合训练 实验报告 班级 08通信1W 姓名丁安华 学号 08313115 指导老师倪福银刘舒淇 2011年09 月 目录 实验一 LED演示 1.1.实验目的 -------------------------------------------------P2 1. 2.实验设备-------------------------------------------------P2 1. 3.实验原理-------------------------------------------------P2 1. 4.实验程序设计流程------------------------------------------P3 1. 5.实验程序编写----------------------------------------------P4 1. 6.实验步骤-------------------------------------------------P7 1. 7.实验结果与分析--------------------------------------------P7实验二键盘输入 2.1.实验目的 -------------------------------------------------P8 2.2.实验设备-------------------------------------------------P8 2. 3.实验原理-------------------------------------------------P8 2. 4.实验程序设计流程------------------------------------------P9 2. 5.实验程序编写----------------------------------------------P10 2. 6.实验步骤-------------------------------------------------P14 2. 7.实验结果与分析--------------------------------------------P14实验三液晶显示器控制显示 3.1.实验目的 -------------------------------------------------P15 3.2.实验设备-------------------------------------------------P15 3.3.实验原理-------------------------------------------------P15 3. 4.实验程序设计流程------------------------------------------P17 3. 5.实验程序编写----------------------------------------------P18 3. 6.实验步骤-------------------------------------------------P22 3. 7.实验结果与分析--------------------------------------------P23实验四有限冲激响应滤波器（FIR）算法 4.1.实验目的 -------------------------------------------------P23 4.2.实验设备-------------------------------------------------P23 4.3.实验原理-------------------------------------------------P24 4.4.实验程序设计流程------------------------------------------P25 4. 5.实验程序编写----------------------------------------------P25 4. 6.实验步骤-------------------------------------------------P27 4. 7.实验结果与分析--------------------------------------------P28 吉大《数字信号处理》（二） 第二章 离散时间信号与系统的z 域分析和频域分析 一、什么是z 变换？ 在时域离散信号和系统中，用序列的傅里叶变换进行频域分析，z 变换则是其推广，用以 对序列进行复频域分析。 序列x(n)的z 变换定义为 式中z 是一个复变量，它所在的复平面称为z 平面。注意在定义中，对n 求和是在－∞、 ＋∞之间求和，可以称为双边Z 变换。还有一种称为单边Z 变换的定义，如下式： 这种单边z 变换的求和限是从零到无限大，因此对于因果序列，用两种Z 变换定义计算的结果是一样的。z 变换是傅里叶变换的推广，傅里叶变换是z 变换的特例—单位圆上的变换。 二、简述z 变换的收敛域。 对于任意给定的序列x(n)，使其z 变换X(z)收敛的所有z 值的集合称为X(z)的收敛域（region of convergence ，Roc ）。使级数收敛的充要条件是满足绝对可和,如下式所示： 一般收敛域用环状域表示如下： 常用的Z 变换是一个有理函数，用两个多项式之比表示为： ()()n n X z x n z ∞ -=-∞ = ∑0 ()()n n X z x n z ∞ -==∑() ()() P z X z Q z = 分子多项式P(z)的根是X(z)的零点，分母多项式Q(z)的根是X(z)的极点。在极点处Z 变换不存在，因此收敛域中没有极点，收敛域总是用极点限定其边界。 三、x(n)=u(n), 求其Z 变换。 解： ()()n n n n X z u n z z ∞ ∞ --=-∞ == =∑∑ X(z)存在的条件是1||1<-z ，因此收敛域为|z|>1, 因此: X(z)表达式表明，极点是z=1，单位圆上的z 变换不存在，或者说收敛域不包含单位圆， 该例同时说明一个序列的傅里叶变换不存在，但在一定的收敛域内z 变换是可以存在的。 四、求)()(n u a n x n =的z 变换及其收敛域。 解： 10 1()( )1n n n n n n X z a u n z a z az ∞ ∞ ---=-∞ == == -∑ ∑ 在收敛域中必须满足1||1<-az ，因此收敛域为|z|>|a|。 五、已知 求z 反变换。 解： 1 1 ()||1 1X z z z -=>- 国贸专业生产实习报告 随着国际贸易的日益完善，以及中国在国际贸易的地位的不断上升，我们作为未来社会的国贸人员，为了加强社会竞争力，应培养较强的国贸工作的操作能力。于是，在结束了大三的课程后,学校给了我们一个很好的实习锻炼机会，让我们模拟国际贸易实务操作，从而从中掌握国际贸易流程。 一、实习目的 ①熟悉外贸实务的具体操作流程; ②了解、巩固与深化已经学过的理论和方法; ③增强对外贸实务的感性认识; ④提高发现问题、分析问题以及解决问题的能力。 二、实习方法： 通过进入SimTrade模拟平台，进行上机模拟操作 Simtrade外贸实习平台是一个十分成功的国际贸易模拟软件，它在很大程度上解决了学生实习难的问题。学生在网上进行国际货物买卖实务的具体操作，能很快掌握进出口的成本核算、询盘、发盘与还盘等各种基本技巧；熟悉国际贸易的物流、资金流与业务流的运作方式；切身体会到国际贸易中不同当事人面临的具体工作与他们之间的互动关系；学会外贸公司利用各种方式控制成本以达到利润最大化的思路；认识供求平衡、竞争等宏观经济现象，并且能够合理地加以利用。老师通过在网站发布新闻等行为对国际贸易环境实施宏观调控，使学生在实习中充分发挥主观能动性，真正理解并吸收课堂中所学到的知识，为将来走上工作岗位打下良好基础。 三、实习遇到的问题 1、预算错误 这是开始接触Simtrade时所最容易忽略的问题。虽然老师曾多次提醒，做贸易前一定要计算好了一切费用，选好贸易术语，最后再签定合同。但我们经常做出口商的在还没有调查进口商所在地市场的情况下就先去工厂进货了。如果工厂角色也没有做好预算，草草就签订了合同，那么可能出口商和工厂都赚不到钱。在最后交易完成后，我们经常大叫“啊，这个运费怎么比我的货物数量还多啊？”“这个保险费怎么这么贵，我要赔钱了！” 实验一 MATLAB 仿真软件的基本操作命令和使用方法 实验容 1、帮助命令 使用 help 命令，查找 sqrt （开方）函数的使用方法； 2、MATLAB 命令窗口 （1）在MATLAB 命令窗口直接输入命令行计算3 1)5.0sin(21+=πy 的值； （2）求多项式 p(x) = x3 + 2x+ 4的根； 3、矩阵运算 （1）矩阵的乘法 已知 A=[1 2;3 4]， B=[5 5;7 8]，求 A^2*B （2）矩阵的行列式 已知A=[1 2 3;4 5 6;7 8 9]，求A （3）矩阵的转置及共轭转置 已知A=[1 2 3;4 5 6;7 8 9]，求A' 已知B=[5+i,2-i,1;6*i,4,9-i], 求B.' , B' （4）特征值、特征向量、特征多项式 已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征值、特征向量、特征多项式； （5）使用冒号选出指定元素 已知：A=[1 2 3;4 5 6;7 8 9]；求A 中第3 列前2 个元素；A 中所有列第2，3 行的元素； 4、Matlab 基本编程方法 （1）编写命令文件：计算1+2+…+n<2000 时的最大n 值； （2）编写函数文件：分别用for 和while 循环结构编写程序，求 2 的0 到15 次幂的和。 5、MATLAB基本绘图命令 （1）绘制余弦曲线 y=cos(t)，t∈[0，2π] （2）在同一坐标系中绘制余弦曲线 y=cos(t-0.25)和正弦曲线 y=sin(t-0.5)， t∈[0，2π] （3）绘制[0，4π]区间上的 x1=10sint 曲线，并要求： （a）线形为点划线、颜色为红色、数据点标记为加号； （b）坐标轴控制：显示围、刻度线、比例、网络线 （c）标注控制：坐标轴名称、标题、相应文本； >> clear; 实验一 程序的控制与转移 一、实验目的 1、掌握条件算符的使用。 2、掌握循环操作指令（BNAZ ）和比较操作指令（CMPR ） 二、实验设备 计算机、ZY13DSP12BD 实验箱、5402EVM 板。 三、实验原理 程序控制指令主要包括分支转移、子程序调用、子程序返回、条件操作及循环操作等。通过传送控制到程序存储器的其他位置，转移会中断连续的指令流。转移会影响在PC 中产生和保护的程序地址。其中转移可以分为两种形式的，一种是有条件的，另一种是无条件的。 四、实验内容 编写程序，实现计算y= ∑=5 1 i i x 的值。 五、实验步骤 1、用仿真机将计算机与ZY13DSP12BD 实验箱连接好，并依次打开实验箱电源、仿真机电源，然后运行CCS 软件。 2、新建一个项目：点击Project －New ，将项目命名为example2，并将项目保存在自己定义的文件夹下。 3、新建一个源文件example2.asm 。将该文件添加到工程example2.pjt 中。 4、在工程管理器中双击example2.asm ，编写源程序： .tiltle ”example2.asm ” .mmregs STACK .usect ”STACK ”,10H ;堆栈的设置 .bss x,5 ;为变量分配6个字的存储空间 .bss y,1 .def start .data table: .word 10,20,3,4,5 ;x1,x2,x3,x4,x5 .text Start: STM #0,SWWWSR ;插入0个等待状态 STM #STACK+10H,sp ;设置堆栈指针 STM #x,AR1 ;AR1指向x RPT #4 ;下一条被重复执行5遍 MVPD table,*AR1+ ;把程序存储器中的数据传送到数据存储器 LD #0,A ;A 清零 CALL SUM ;调用求和函数 end: B end SUM: STM #x,AR3 ;AR3指向x STM #4,AR2 ;AR2=4 loop: ADD *AR3+,A ;*AR3+A-->A,然后AR3+ BANZ loop,*AR2- ;如果AR2的值不为0,则跳到loop 处;否则执行下一条指令 STL A,*(y) ;把A 的低16位赋给变量y 信息科学与技术学院本科三年级 数字信号处理实验报告 2011 年12 月21日 实验一 序列的傅立叶变换 实验目的 进一步加深理解DFS,DFT 算法的原理;研究补零问题;快速傅立叶变换 （FFT ）的应用。 实验步骤 1. 复习DFS 和DFT 的定义，性质和应用； 2. 熟悉MATLAB 语言的命令窗口、编程窗口和图形窗口的使用；利用提供的 程序例子编写实验用程序；按实验内容上机实验，并进行实验结果分析；写出完整的实验报告，并将程序附在后面。 实验内容 1. 周期方波序列的频谱试画出下面四种情况下的的幅度频谱,并分析补零后，对信号频谱的影响。 实验结果： 60 ,7)4(;60,5)3(; 40,5)2(;20,5)1()] (~[)(~,2,1,01 )1(,01,1)(~=========±±=???-+≤≤+-+≤≤=N L N L N L N L n x DFS k X m N m n L m N L m N n m N n x ) 52.0cos()48.0cos()(n n n x ππ+= 2. 有限长序列x(n)的DFT （1） 取x(n)(n=0:10)时，画出x(n)的频谱X(k) 的幅度； （2） 将（1）中的x(n)以补零的方式，使x(n)加长到(n:0~100)时，画出 x(n)的频谱X(k) 的幅度； （3） 取x(n)(n:0~100)时，画出x(n)的频谱X(k) 的幅度。利用FFT 进行谱分析 已知：模拟信号 以t=0.01n(n=0:N-1)进行采样，求N 点DFT 的幅值谱。 请分别画出N=45; N=50;N=55;N=60时的幅值曲线。 实验结果： ) 8cos(5)4sin(2)(t t t x ππ+= 国际经济与贸易专业 外贸交易模拟实验 实验报告 学号__ _ 姓名__ _____ 班级_______ 指导老师___ _______ 实验地点 一、试验目的 在Simtrade这个虚拟贸易平台中，通过扮演不同的角色，熟练掌握各种业务技巧，了解到国际贸易的物流、资金流和业务流的运作方式，增强感性认识和社会适应能力，进一步巩固、深化已学过的理论知识，提高综合运用所学知识发现问题、解决问题的能力。 二、实习时间 本次实习时间为2011-03-01至2011-03-28，共计四周 三、试验小结 本次试验运用交易方式:①L/C+CIF②L/C+CFR③T/T+CIF④T/T+CFR⑤ D/P+FOB⑥T/T+FOB 涉及到的交易产品：①洋菇罐头( 柄及碎片)CANNED STEMS AND PIECES MUSHROOMS②甜玉米罐头CANNED SWEET CORN③荔枝罐头CANNED LITCHIS④芒果罐头CANNED MANGOES⑤名牌手提包FAMOUS-BRAND HANDBAG⑥香味蜡烛SCENTED CANDLE 四、试验遇到的问题及其解决方法 1.预算与实际支出相差较大 ?没有换算成本币 在CONTRACT111中，我的进口预算表和是实际发生额都有２+汇率的差异，原因就在于我公司注册资金为欧元，而样本中为美元，而我当初为了省时省力省脑，就全部依样画弧，没有转换成欧元导致。 ?集装箱计算的问题 在CONTRACT4442中，我出口商出口预算表海运费在计算过程中集装箱数为4，而实际确实12个集装箱，结果海运费预算和实际相差甚多，之后导致FOB 价格、利润等一系列数字有出入。 ?粗心大意 东南大学自动化学院 实验报告 课程名称： D SP 原理及C 程序开发 第二次实验 实验名称：基于DSP 系统的实验——指示灯、拨码开关和定时器院（系）：自动化专业：自动化 姓名：学号： 实验室：实验组别： 同组人员：实验时间：2012 年 4 月 18日 评定成绩：审阅教师： 第一部分实验：基于DSP 系统的实验——指示灯和拨码开关 一．实验目的 1. 了解ICETEK –F28335-A 评估板在TMS320F28335DSP 外部扩展存储空间上的扩展。 2. 了解ICETEK –F28335-A 评估板上指示灯和拨码开关扩展原理。 3. 学习在C 语言中使用扩展的控制寄存器的方法。 二．实验设备 计算机，ICETEK –F28335-A 实验箱（或ICETEK 仿真器+ICETEK–F28335-A 评估板+相关连线及电源）。 三．实验原理 1．TMS320F28335DSP 的存储器扩展接口 存储器扩展接口是DSP 扩展片外资源的主要接口，它提供了一组控制信号和地址、数据线，可以扩展各类存储器和存储器、寄存器映射的外设。 -ICETEK –F28335-A 评估板在扩展接口上除了扩展了片外SRAM 外，还扩展了指示灯、DIP 开关和D/A 设备。具体扩展地址如下： 0x180004- 0x180005：D/A 转换控制寄存器 0x180001：板上DIP 开关控制寄存器 0x180000：板上指示灯控制寄存器 -与ICETEK –F28335-A 评估板连接的ICETEK-CTR 显示控制模块也使用扩展空间控制主要设备： 208000-208004h ：读-键盘扫描值，写-液晶控制寄存器 208002-208002h ：液晶辅助控制寄存器 208003-208004h ：液晶显示数据寄存器 2．指示灯与拨码开关扩展原理 页码，1/3 吉林大学牡丹园站 -- Construction精华区文章阅读 发信人: arwang (旺旺), 信区: Construction 标 题: 流体力学核心期刊 发信站: 牡丹园新站 (Sun Dec 21 09:18:46 2003) 流体力学核心期刊 Journal of Fluid Mechanics = 流体力学杂志 . 英国.527C0001 International Journal of Heat and mass Transfer = 国际传热与传质杂志 \ 英国 .525C0006 AIAA Journal = 美国航空与航天学会志 . 美国.877B0001 The Physics of Fluids, A = 流体物理学，A辑 . 美国.527B0002 Fluids Dynamics = 流体动力学 （ 英译苏刊）. 美国.527B0054 Journal of Engineering Physics = 工程物理杂志（英译苏刊）. 美国.534B0053 Journal of Heat transfer,Transactions of the ASME = 传热杂志，ASME汇刊 . 美 国.725B0001 The Physics of Fluids, B = 流体物理学，B辑 . 美国.527B0002 International Journal for Numerical Methods in Fluids = 国际流体力学数值方法 杂志 . 英国.527C0004 Fluid MechanicsSoviet Research = 苏联流体力学研究（英译苏刊） . 美国.527B005 2 International Journal of Multiphase flow = 国际多相流杂志 . 英国.527C0003 Zeitschrift fur Angewandte Mathematik und Mechanik = 应用数学与力学杂志 . 德 国.519A0001 Magnetohydrodynamics = 磁流体动力学（英译苏刊）. 美国.527B0053 Journal of Applied Mechnaics and Technical physics = 应用力学与技术物理杂志（ 英 译苏刊）. 美国.529B0052 Journal of Fluids Engineering, Transactions of the ASME = 流体工程杂志，ASME 汇刊 . 美国.780B0001 Physical Review , A = 物理评论，A辑 . 美国.530B0002 Soviet PhysicsDOKLADY = 苏联物理学报告（英译苏刊）. 美国.530B0070 International Journal of Heat and Fluid Flow = 国际热与流体流杂志 . 英国.527 C0053 Journal of Non-Newtonian Fluid Mechanics = 非牛顿流体动力学杂志 . 荷兰.527LB 0053 International Communications in Heat and Mass Transfer = 国际传热与传质通讯 . 英 国.725C0056 Heat Transfer Soviet Research = 苏联传热研究 . 美国.725B0054 Physical Review Letters = 物理评论快报 . 美国.530B0003 International Journal of Engineering Science = 国际工程科学杂志 . 英国.710C0 009 Journal of Computational Physics = 计算物理杂志 . 美国.539B0002 Waerme-und Stoffuebertragung = 热力学与流体力学 . 德国.710E0008 Physica,D = 物理，D辑 . 荷兰.530LB001 High Temperature = 高温（英译苏刊）. 美国.534B0052 JSME International Journal, II = 日本机械工程师学会国际杂志，II辑 . 日本.780 D0063 Fluid Dynamics Research = 流体动力学研究 . 荷兰.527LB001 Journal of the Physical Society of Japan = 日本物理学会志 . 日本.530D0002 Computers and Fluids = 计算机与流体 . 英国.--─738C0074 Heat Transfer-Japanese Research = 日本传热研究 . 美国.525B0055 Chemical Engineering Science = 化学工程科学 . 英国.810C0004 Physics Letters, A = 物理快报，A辑 . 荷兰.530LB004 Thermal Engineering = 热力工程（英译苏刊）. 英国 .721C0058 AIChE Journal美国化学工程师协会会志 . 美国.810B0001 Applied Mathematics and Mechanics = 应用数学与力学（英译苏刊）. 美国.ISSN 00 66-5479 Applied Scientific Research = 应用科学研究 . 荷兰.500LB002 Comptes Rendus de l Acadecie des Sciences , Serie II = 法国科学院报告，II辑 . FRA.500F0003 Numerical Heat Transfer = 数值传热 . 美国.725B0059 Rheologica Acta = 流变学学报 . 德国.526E0051吉林大学本科“十二五”规划教材修订教材选题立项名单.doc
simtrade实训总结
数字信号处理实验报告
dsp课程设计实验报告
数字信号处理实验报告
dsp课程设计实验报告总结
吉林大学培养方案
simtrade外贸实务实训报告
数字信号处理实验报告一
DSP实验报告
吉林大学《数字信号处理》期末考试学习资料(二)
Simtrade实验报告
数字信号处理实验报告(实验1_4)
DSP实验报告
西南交大数字信号处理报告
simtrade实验报告
DSP技术与课程设计实验报告二(精)
流体力学核心期刊
相关主题
文本预览