第二章
2.25 已知线性时不变系统的差分方程为
若系统的输入序列x(x)={1,2,3,4,2,1}编写利用递推法计算系统零状态响应的MATLAB程序，并计算出结果。
代码及运行结果：
>> A=[1,-0.5];
>> B=[1,0,2];
>> n=0:5;
>> xn=[1,2,3,4,2,1];
>> zx=[0,0,0];zy=0;
>> zi=filtic(B,A,zy,zx);
>> yn=filter(B,A,xn,zi);
>> figure(1)
>> stem(n,yn,'.');
>> grid on;
2.28图所示系统是由四个子系统T1、T2、T3和T4组成的，分别用单位脉冲响应或差分方程描述为
T1：
其他
T2:
其他
T3：
T4:
编写计算整个系统的单位脉冲响应h(n)，0≤n≤99的MATLAB程序，并计算结果。
代码及结果如下：
>> a=0.25;b=0.5;c=0.25;
>> ys=0;
>> xn=[1,zeros(1,99)];
>> B=[a,b,c];
>> A=1;
>> xi=filtic(B,A,ys);
>> yn1=filter(B,A,xn,xi);
>> h1=[1,1/2,1/4,1/8,1/16,1/32]; >> h2=[1,1,1,1,1,1];
>> h3=conv(h1,h2);
>> h31=[h3,zeros(1,89)]; >> yn2=yn1+h31;
>> D=[1,1];C=[1,-0.9,0.81]; >> xi2=filtic(D,C,yn2,xi); >> xi2=filtic(D,C,ys);
>> yn=filter(D,C,yn2,xi); >> n=0:99;
>> figure(1)
>> stem(n,yn,'.');
>> title('单位脉冲响应'); >> xlabel('n');ylabel('yn');
2.30 利用MATLAB画出受高斯噪声干扰的正弦信号的波形，表示为
其中v(n)是均值为零、方差为1的高斯噪声。
代码及结果如下：
>> N=100;
>> n=0:N;
>> xn=10*sin(0.02*pi*n);
>> R=randn(1,N+1);
>> x=xn+R;
>> figure(2);
>>plot(n,x,'.'),title('受高斯噪声干扰的正弦信号'),xlabel('n'),ylabel('x');
第三章
3.47 利用Matlab 工具箱函数zplane(b,a),画出下列Z 变换的零极点分布图，并给出所有可能的收敛域及对应序列的特性（左边序列，右边序列，双边序列）。
（1）12
181533325644162)(2
342341-+-+++++=z z z z z z z z z X
>> b=[2,16,44,56,32]; >> a=[3,3,-15,18,-12]; >> zplane(b,a)
（2）65
610204
.874.2698.1768.84)(2342342+++--+--=z z z z z z z z z X
>> b=[4,-8.68,-17.98,26.74,-8.04]; >> a=[1,-2,10,6,65]; >> zplane(b,a)
3.53 利用Matlab 语言，画出下列无限长脉冲响应系统的幅频响应特性曲线和相频响应特性曲线，并指出系统的类型。
（1）)
61.088.01)(4.01()1()(2
112
111-----+---=z z z z z z H >> syms z;
>> ps=z^-1*(1-z^-1)^2; >> ps1=expand(ps) ps1 =
1/z - 2/z^2 + 1/z^3 >> syms z;
>> ps = (1 - 0.4*z^ - 1)*(1 - 0.88*z^ - 1 + 0.61*z^ - 2); >> ps1=expand(ps) ps1 =
481/(500*z^2) - 32/(25*z) - 61/(250*z^3) + 1 >> a=[1,-32/25,481/500,-61/250]; >> b=[0,1,-2,1];
>> [H,w]=freqz(b,a,'whole');
>> subplot(2,1,1),plot(w/pi,abs(H));
>> xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|') >> subplot(2,1,2),plot(w/pi,angle(H)/pi);
>> xlabel('\omega/\pi');ylabel('\phi(\omega)/\pi')
(2) )
7957.04461.1.1)(683.01()0166.11)(1(0534.0)(2112
2112------+-+-+=z z z z z z z H
>> syms z;
>> ps=0.0534*(1+z^-1)*(1-1.0166*z^-1+z^-2)^2; >> ps1=expand(ps) ps1 =
6676839363/(125000000000*z^2) - 689661/(12500000*z) + 6676839363/(125000000000*z^3) - 689661/(12500000*z^4) + 267/(5000*z^5) + 267/5000
>> syms z;
>> ps=(1-0.683*z^-1)*(1-1.4461*z^-1+0.7957*z^-2);
>> ps1=expand(ps)
ps1 =
17833863/(10000000*z^2) - 21291/(10000*z) - 5434631/(10000000*z^3) + 1
>> a=[1,- 21291/10000,17833863/10000000, - 5434631/10000000,0,0];
>> b=[267/5000,- 689661/12500000,6676839363/125000000000,6676839363/125000000000,- 689661/12500000,267/5000];
>> [H,w]=freqz(b,a,'whole');
>> subplot(2,1,1),plot(w/pi,abs(H));
>> xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|')
>> subplot(2,1,2),plot(w/pi,angle(H)/pi);
>> xlabel('\omega/\pi');ylabel('\phi(\omega)/\pi')
（3）)
6493.05548.11)(8482.0499.11()1()(21214
13-----+-+--=z z z z z z H
>> syms z;
>>ps=(1-z^-1)^4; >>ps1=expand(ps) ps1 =
6/z^2 - 4/z - 4/z^3 + 1/z^4 + 1 >> syms z;
>> ps=(1-1.499*z^-1+0.8482*z^-2)*(1-1.5548*z^-1+0.6493*z^-2); >> ps1=expand(ps) ps1 =
9570363/(2500000*z^2) - 15269/(5000*z) - 114604103/(50000000*z^3) + 27536813/(50000000*z^4) + 1
>> a=[1,- 15269/5000,9570363/2500000,- 114604103/50000000,27536813/50000000]; >> b=[1,-4,6,-4,1];
>> [H,w]=freqz(b,a,'whole');
>> subplot(2,1,1),plot(w/pi,abs(H));
>> xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|') >> subplot(2,1,2),plot(w/pi,angle(H)/pi);
>> xlabel('\omega/\pi');ylabel('\phi(\omega)/\pi')
第四章
4.32 已知两个序列分别为
其他
其他
利用MATLAB，采用离散傅里叶变换计算出与的4点循环卷积。验证傅里叶变换的时域循环卷积定理。
代码及运行结果如下：
>> h=[2,1,1,2];
>> x=[1,-1,-1,1];
>> y=ifft(fft(h,4).*fft(x,4))
y =
0 -2 0 2
4.34 选择合适的变换区间长度N，编写matlab程序，用离散傅里叶变换对下列信号进行谱分析，画出幅频特性曲线和相频特性曲线。
（1）
（2）
代码及运行结果如下：
（1）
>> n=0:9;
>> x1=2*sin(0.2*pi*n);
>> X1=fft(x1);
>> subplot(2,1,1),stem(n,abs(X1)); >> title('幅频特性曲线');
>> ylabel('幅频'),xlabel('k');
>> subplot(2,1,2),stem(n,angle(X1)); >> ylabel('相频'),xlabel('k');
（2））
>> n=0:50;
>> x1=sin(0.45*pi*n).*sin(0.55*pi*n);
>> X1=fft(x1);
>> subplot(2,1,1),stem(n,abs(X1),'.');
>> ylabel('幅频'),xlabel('k');
>> title('幅频特性曲线');
>> subplot(2,1,2),stem(n,angle(X1),'.');
>> title('相频特性曲线');
>> ylabel('相频'),xlabel('k');
第五章
5.13 设连续时间信号,其中，,
,。
（1）如果FFT对进行频谱分析，问采样频率和采样点数N应如何选择，才能准确地求出，，的中心频率，为什么？
（2）按照所选择的和N，对进行采样，得序列x(n)。试用MATLAB工具箱函数fft，对x(n)做FFT，得到X(k)，并画出X(k)的幅频特性| X(k) |的曲线，标出| X(k) |的三个峰值所对应的k值分别是多少？
代码及运行结果如下：
（1）由奈奎斯特采样定理，所以≥20，取=40，为了分辨相邻谱峰f=2，所以N==20
代码及运行结果如下：
>> fs=40;%采样频率
>> f=1;
>> N=fs;
>> t=0:1/(fs-1):1;
>> x=cos(8*pi*t)+cos(16*pi*t)+cos(20*pi*t);
>> y=fft(x);
>> a=abs(y);
>> n=0:length(y)-1;
>> subplot(2,1,1);stem(t,x);
>> title('x(n)离散');
>> subplot(2,1,2);stem(n,a);
>> title('幅频特性');
>> gtext('k1');
>> gtext('k2');
>> gtext('k3');
5.16 设连续时间信号，试用MATLAB语言分析该信号的频率特性，要求信号的频率分辨率f≤0.5Hz，画出其幅频特性曲线。请分析误差来源，指出减小误差的方法。
代码及运行结果如下：
>> fs=100;%采样频率
>> f=0.4;
>> N=250;
>> t=0:1/(fs-1):2.5;
>> x=1+cos(100*pi*t);
>> y=fft(x);
>> a=abs(y);
>> n=0:length(y)-1;
>> subplot(2,1,1);stem(t,x);
>> title('x(n)');
>> subplot(2,1,2);stem(n,a);
>> title('幅频响应');
第六章
6.32 利用MATLAB工具箱函数bessel，设计三阶模拟低通贝塞尔滤波器，要求通带截止频率。画出滤波器的幅频响应特性曲线和相频响应特性曲线。
代码及运行结果如下：
>> [Bz,Az]=besself(3,2*pi*2000);
>> [H,W]=freqs(Bz,Az);%计算频率向量w上的滤波器频率响应
>> a=20*log10(abs(H));
>> figure(1)
>> plot(W,a)%幅频响应特性
>> title('幅频特性响应');
>> b=angle(H);
>> figure(2);
>> plot(W,b);
>> title('相频响应曲线');
6.37 利用MATLAB工具箱函数，采用模拟低通切比雪夫1型滤波器及脉冲响应不变法，设计满足下列指标的数字低通滤波器，并画出幅频响应曲线。
通带截止频率，阻带截止频率；
通带最大衰减，阻带最小衰减。
代码及运行结果如下：
Wp=0.2;Ws=0.4;Rp=1;Rs=15;
[N,Wpo]=cheb1ord(Wp,Ws,Rp,Rs,'s');
[Bs,As]=cheby1(N,Rp,Wpo,'s');
[Bz,Az]=impinvar(Bs,As);
w=0:0.1:pi;
[H,w1]=freqz(Bz,Az,w);
H=20*log10(abs(H));
plot(w/pi,H),grid on;
xlabel('\omega/\pi');ylabel('|H(e^j^\omega)|/dB');
第七章
7.11 用矩形窗设计线性相位FIR数字低通滤波器，逼近理想低通滤波器的为
（1）求出所设计的单位脉冲响应h（n）的表示式，确定与h（n）的长度N的关系式。（2）取N=31，，用MATLAB画出该FIR数字低通滤波器的幅频响应和相频响应特性曲线。
运行结果如下：
7.13 希望对输入低频模拟信号做数字低通滤波处理，以提取所需要的信号，设系统的采样频率为50hz，要求通带截止频率为10kHZ，阻带截止频率为25kHZ，阻带最小衰减为60dB。（1）用窗函数设计FIR数字低通滤波器，选择合适的窗函数及滤波器长度N，求出单位脉冲响应h（n）。
（2）用matlab画出该FIR数字低通滤波器的扶贫响应特性曲线和相频特性曲线。
运行结果如下：
7.21 利用matlab工具箱函数，分别画出长度为15的矩形窗，汉宁窗，汉明窗和布莱克窗的幅频特性曲线，观察他们的主要参数（主瓣宽度，过渡带宽度，旁瓣峰值幅度）的差别。运行结果如下：
（1）矩形窗
（2）汉宁窗
（3）汉明窗
（4）布莱克曼窗
由图可以看出旁瓣峰值幅度，矩形窗，汉宁窗，汉明窗，布莱克窗是依次减小的而主瓣宽度，和过渡带都是依次加宽的。
7.22 利用matlab工具箱函数firl，设计线性相位FIR数字低通滤波器，要求希望逼近的理想低通滤波器的通带截止频率wp=pi/4rad，滤波器长度N=21，试分别选用矩形窗，汉宁窗，汉明窗和布莱克窗进行设计，画出用每种窗函数所设计的滤波器的单位脉冲响应h（n）波形及幅频响应特性曲线，并进行比较，观察各种窗函数的设计性能。
运行结果如下：
（1）矩形窗
随机信号大作业 班级：02xxxx 姓名：xx
学号：02xxxxx 第一章 1.23上机题：设有随机初相信号X(t)=5cos（t+φ），其中相位φ是在区间（0,2π）上均匀分布的随机变量。试用Matlab编程产生其三个样本函数。 解：程序： clc clear m=unifrnd(0,2*pi,1,10); for k=1:3 t=1:0.1:10; X=5*cos(t+m(k)); plot(t,X); hold on
end title('其三个样本函数'); xlabel('t');ylabel('X(t)'); grid on ;axis tight ; 由 Matlab 产生的三个样本函数如下图所示： 第二章 2.22 上机题：利用Matlab 程序设计一正弦型信号加高斯白噪声的复合信号。 (3)分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。 1 2 3 4 5 6 7 8 9 10 -4-3-2-101 23 4其三个样本函数 t X (t )
解：取数据如下： 正弦信号的频率为：fc=10HZ，抽样频率为：fs=100HZ； 信号：x=sin(2*pi*fc*t)； 高斯白噪声产生复合信号y： y=awgn(x,10)； 复合信号y通过理想滤波器电路后得到信号y3 ，通过卷积计算可以得到y3 即：y3=conv2(y,sin(10*t)/(pi*t))； y3的幅度分布特性可以通过傅里叶变换得到Y3(jw)=fft(y3)，y3的功率谱密度：G3(w)=Y3(jw).*conj(Y3(jw)/length(Y3(jw)))。 程序： clear all; fs=100; fc=10; n=201; t=0:1/fs:2; x=sin(2*pi*fc*t); y=awgn(x,10); m=50; i=-0.49:1/fs:0.49; for j=1:m R(j)=sum(y(1:n-j-1).*y(j:199),2)/(n-j); Ry(49+j)=R(j);
射频大作业 基于PSpice仿真的振幅调制电路设计数字调制与解调的集成器件学习
目录 题目一：基于PSpice仿真的振幅调制电路设计与性能分析 一、实验设计要求 (3) 二、理论分析 1、问题的分析 (3) 2、差动放大器调幅的设计理论 (4) 2.1、单端输出差动放大器电路 2.2、双端输出差动放大器电路 2.3、单二极管振幅调制电路 2.4、平衡对消二极管调幅电路 三、PSpice仿真的振幅调制电路性能分析 (10) 1、单端输出差动放大器调幅电路设计图及仿真波形 2、双端输出差动放大器调幅电路设计图及仿真波形 3、单二极管振幅调制电路设计图及仿真波形 4、平衡对消二极管调幅电路设计图及仿真波形 四、实验总结 (16) 五、参考文献 题目二数字调制与解调的集成器件学习 一、实验设计要求 (17) 二、概述 (17) 三、引脚功能及组成原理 (18) 四、基本连接电路 (20) 五、参考文献 (21) 六、英文附录 (21)
题目一基于PSpice仿真的振幅调制电路设计 摘要 随着大规模集成电路的广泛发展,电子电路CAD及电子设计自动化（EDA）已成为电路分析和设计中不可缺少的工具。此次振幅调制电路仿真设计基于PSpice，利用其丰富的仿真元器件库和强大的行为建模工具，分别设计了差分对放大器和二极管振幅调制电路，由此对线性时变电路调幅有了更进一步的认识；同时，通过平衡对消技术分别衍生出双端输出的差分对放大器和双回路二极管振幅调制电路，消除了没用的频率分量，从而得到了更好的调幅效果。本文对比研究了单端输出和双端输出的差分对放大器调幅电路及单二极管和双回路二极管调幅电路，通过对比观察时域和频域波形图，可知平衡对消技术可以很好地减小失真。 关键词：PSpice 振幅调制差分对放大器二极管振幅调制电路平衡对消技术 一、实验设计要求 1.1 基本要求 参考教材《射频电路基础》第五章振幅调制与解调中有关差分对放大器调幅和二极管调幅的原理，选择元器件、调制信号和载波参数，完成PSpice电路设计、建模和仿真，实现振幅调制信号的输出和分析。 1.2 实践任务 (1) 选择合适的调制信号和载波的振幅、频率，通过理论计算分析，正确选择晶体管和其它元件；搭建单端输出的差分对放大器，实现载波作为差模输入电压，调制信号控制电流源情况下的振幅调制；调整二者振幅，实现基本无失真的线性时变电路调幅；观察记录电路参数、调制信号、载波和已调波的波形和频谱。 (2) 参考例5.3.1，修改电路为双端输出，对比研究平衡对消技术在该电路中的应用效果。 (3) 选择合适的调制信号和载波的振幅、频率，通过理论计算分析，正确选择二极管和其它元件；搭建单二极管振幅调制电路，实现载波作为大信号，调制信号为小信号情况下的振幅调制；调整二者振幅，实现基本无失真的线性时变电路调幅；观察记录电路参数、调制信号、载波和已调波的波形和频谱。 (4) 参考例5.3.2，修改电路为双回路，对比研究平衡对消技术在该电路中的应用效果。 1.3 写作报告 (1) 按论文形式撰写，包括摘要、正文和参考文献，等等。 (2) 正文包括振幅调制电路的设计原理、理论分析结果、实践任务中各阶段设计的电路、参数、波形和频谱，对观察记录的数据配以图像和表格，同时要有充分的文字做分析和对比，有规律性认识。 (3) 论文结构系统、完备、条理清晰、理论正确、数据翔实、分析完整。 1.4 相关提示 (1) 所有电路和信号参数需要各人自行决定，各人有不同的研究结果，锻炼学生的独立研究和实验分析能力。 (2) 为了提高仿真精度和减小调试难度，可以将调制信号和载波的频率设置得较低。 二、理论分析 1、问题的分析 根据题目的要求，差分对放大器和二极管振幅调制电路目的都是实现基本无
题 目： 数据挖掘 学 院： 电子工程学院 专 业： 智能科学和技术 学生姓名： ** 学 号： 02115*** k -means 实验报告 一、 waveform 数据 1、 算法描述 1. 从数据集{X n }n?1N 中任意选取k 个赋给初始的聚类中心c 1, c 2, …,
c k； 2.对数据集中的每个样本点x i，计算其和各个聚类中心c j的欧氏 距离并获取其类别标号： label(i)=arg min ||x i?c j||2,i=1,…,N,j=1,…,k 3.按下式重新计算k个聚类中心； c j=∑x j s:label(s)=j j ,j=1,2,…k 重复步骤2和步骤3，直到达到最大迭代次数为止2、实验结果 二、图像处理 1、算法描述 同上； 2、实验结果
代码： k_means: %%%%%%%%%K_means%%%%%%%% %%%%%%%%%函数说明%%%%%%%% %输入: % sample——样本集； % k ——聚类数目； %输出： % y ——类标(从0开始) % cnew ——聚类中心 % n ——迭代次数 function [y cnew n]=k_means(sample,k) [N V]=size(sample); %N为样本的个数 K为样本的维数 y=zeros(N,1); %记录样本类标 dist=zeros(1,k); rand_num=randperm(N); cnew=(sample(rand_num(1,1:k),:));%随机初始化聚类中心cold=zeros(k,V); n=0;
数字信号处理上机第一次实验 实验一： 设给定模拟信号()1000t a x t e -=，的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分 量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的进行采样。 ○1 。 ○2 。 比较两种采样率下的信号频谱，并解释。 实验一MATLAB 程序： （1） ○ 1 clc; fs=5000; ts=1/fs; N=1000; t=(-N:N)*ts; s=exp(-abs(t)); plot(t,s,'linewidth',1.5) xlabel('时间') ylabel('幅度') set(gca,'fontweight','b','fontsize',12) SPL=N*100; figure sp=fftshift(fft(s,SPL)); sp=sp/max(sp)*100; freqb=-fs/2:fs/SPL:fs/2-fs/SPL; plot(freqb,abs(sp)) xlabel('频率') ylabel('频谱幅度') set(gca,'fontweight','b','fontsize',12) yy=abs(abs(sp)-3); [aa,freqind]=min(yy); (freqind-SPL/2)*fs/SPL t ()a x t ()()15000s a f x t x n =以样本秒采样得到。()() 11j x n X e ω画出及其频谱()()11000s a f x t x n =以样本得到。()()11j x n X e ω画出及其频谱
随机信号分析大作业 学院：电子工程学院 班级：021151 学号：02115037 姓名：隋伟哲
第一题：设有随机信号X（t）=5cos(t+a),其中相位a是在区间（0,2π）上均匀分布的随机变量，使用Matlab编程产生其三个样本函数。 解： 源程序如下： clc;clear; C=2*pi*rand(1,3);%在[0,2π]产生均匀分布的相位角 t=1:.1:80; y1=5*cos(t+C(1)); %将产生的随机相位角逐一代入随机过程中 y2=5*cos(t+C(2)); %将产生的随机相位角逐一代入随机过程中 y3=5*cos(t+C(3)); %将产生的随机相位角逐一代入随机过程中 plot(t,y1,'r-'); hold on; plot(t,y2,'g--'); hold on; plot(t,y3,'k-'); xlabel('t');ylabel('X(t)'); grid on;axis([0 30 -8 8]); title('随机相位的三条样本曲线'); 产生的三条样本曲线：
第二题：利用Matlab程序设计一正弦型信号加高斯白噪声的复合信号。（1）分析复合信号的功率谱密度、幅度分布特性； （2）分析复合信号通过RC积分电路后的功率谱密度和相应的幅度分布特性； （3）分析复合信号通过理想低通系统后的功率谱密度和相应的幅度分布特性。 解：设定正选信号的频率为10HZ，抽样频率为100HZ x=sin(2*pi*fc*t)
（1）正弦函数加上高斯白噪声： y=awgn(x,10) y 的幅度分布特性可以通过傅里叶变换得到： Y(jw)=fft(y) y 的功率谱密度： G(w)=Y(jw).*conj(Y(jw)/length(Y(jw))) 随机序列自相关函数的无偏估计公式为： 1 01()()()N m xx n R m x n x n m N m --==+-∑ 01m N ≤≤- （2）复合信号 y 通过RC 积分电路后得到信号y2 通过卷积计算可以得到y2 即：y2= conv2(y,b*pi^-b*t) y2的幅度分布特性可以通过傅里叶变换得到： Y2(jw)=fft(y2) y2的功率谱密度： G2(w)=Y2(jw).*conj(Y2(jw)/length(Y2(jw))) （3）复合信号 y 通过理想滤波器电路后得到信号y3 通过卷积计算可以得到y3 即：y3=conv2(y,sin(10*t)/(pi*t)) y3的幅度分布特性可以通过傅里叶变换得到： Y3(jw)=fft(y3) y3的功率谱密度： G3(w)=Y3(jw).*conj(Y3(jw)/length(Y3(jw)))
数字水印技术 一、引言 随着互联网广泛普及的应用，各种各样的数据资源包括文本、图片、音频、视频等放在网络服务器上供用户访问。但是这种网络资源的幵放也带了许多弊端，比如一些用户非法下载、非法拷贝、恶意篡改等，因此数字媒体内容的安全和因特网上的侵权问题成为一个急需解决的问题。数字水印作为一项很有潜力的解决手段，正是在这种情况下应运而生。 数字水印（技术是将一些代表性的标识信息，一般需要经过某种适合的变换，变换后的秘密信息（即数字水印），通过某种方式嵌入数字载体（包括文档、音频、软件等）当中，但不影响原载体的使用价值，也不容易被人的知觉系统（如视觉或听觉系统）觉察或注意到。通过这些隐藏在载体中的信息，可以达到确认内容创建者、购买者、传送隐秘信息或者判断载体是否被篡改等目的。在发生产权和内容纠纷时，通过相应的算法可以提取该早已潜入的数字水印，从而验证版权的归属和内容的真伪。 二.算法原理 2.1、灰度图像水印 2.1.1基本原理 处理灰度图像数字水印，采用了LSB（最低有效位）、DCT变换域、DWT变换域三种算法来处理数字水印。在此过程中，处理水印首先将其预处理转化为二值图像，简化算法。 （1）LSB算法原理：最低有效位算法(Least Sig nificant Bit , LSB)是很常见的空间域信息隐藏算法, 该算法就是通过改变图像像素最不重要位来达到嵌入隐秘信息的效果, 该方法隐藏的信息在人的肉眼不能发现的情况下, 其嵌入方法简单、隐藏信息量大、提取方法简单等而获得广泛应用。LSB 信息嵌入过程如下: S′=S+f S ,M 其中,S 和S′分别代表载体信息和嵌入秘密信息后的载密信息;M为待嵌入的秘密信息, 而隐写分析则是从S′中检测出M以至提取M 。 （2）DCT算法原理：DCT 变换在图像压缩中有很多应用，它是JPEG，MPEG 等数据
说明 每个实验题目含有一个main函数和一些函数，与实验题目相关的基本运算的函数定义和main函数定义的代码在附录以及对应的文件夹中给出，供上机实验参考使用。对于每个题目，只需要根据题目要求设计算法，补充函数定义，然后对程序进行编译、调试。
实验一线性表 一、实验目的 1．熟悉线性表的顺序和链式存储结构 2．掌握线性表的基本运算 3．能够利用线性表的基本运算完成线性表应用的运算 二、实验内容 1．设有一个线性表E={e1, e2, … , e n-1, e n}，设计一个算法，将线性表逆置，即使元素排列次序颠倒过来，成为逆线性表E’={ e n , e n-1 , … , e2 , e1 }，要求逆线性表占用原线性表空间，并且用顺序表和单链表两种方法表示，分别用两个程序来完成。（文件夹：顺序表逆置、单链表逆置） 2．已知由不具有头结点的单链表表示的线性表中，含有三类字符的数据元素（字母、数字和其他字符），试编写算法构造三个以循环链表表示的线性表，使每个表中只含有同一类的字符，且利用原表中的结点空间，头结点可另辟空间。（文件夹：分解单链表） 实验二栈和队列 一、实验目的 1．熟悉栈和队列的顺序和链式存储结构 2．掌握栈和队列的基本运算 3．能够利用栈和队列的基本运算完成栈和队列应用的运算 二、实验内容 1．设单链表中存放有n个字符，试编写算法，判断该字符串是否有中心对称的关系，例如xyzzyx是中心对称的字符串。（提示：将单链表中的一半字符先依次进栈，然后依次出栈与单链表中的另一半字符进行比较。）（文件夹：判字符串中心对称） 2．假设以数组sequ[m]存放循环队列的元素，同时设变量rear和quelen 分别指示循环队列中队尾元素的位置和内含元素的个数。编写实现该循环队列的入队和出队操作的算法。 提示：队空的条件：sq->quelen==0；队满的条件：sq->quelen==m。（文件夹：循环队列）实验三串 一、实验目的 1．熟悉串的顺序存储结构 2．掌握串的基本运算及应用 二、实验内容 1．串采用顺序存储结构，编写朴素模式匹配算法，查找在串中是否存在给定的子串。（文件夹：模式匹配） 2．若S是一个采用顺序结构存储的串，利用C的库函数strlen和strcpy（或strncpy）编写
随机信号分析实验报告 信息25班 2120502123 赵梦然
作业题三： 利用Matlab 产生一个具有零均值、单位方差的的高斯白噪声随机序列X(n)，并通过一脉冲响应为 (0.8)(0)0 n n h n else =≥??? 的线性滤波器。 (1) 产生一个具有零均值、单位方差的的高斯白噪声随机序列X(n)，检验其一维概率密度函 数是否与理论相符。 (2) 绘出输入输出信号的均值、方差、自相关函数及功率谱密度的图形，讨论输出信号服从 何种分布。 (3) 试产生在[-1,+1]区间均匀分布的白噪声序列，并将其替换高斯白噪声通过上述系统。 画出此时的输出图形，并观察讨论输出信号服从何种分布。 作业要求 (1) 用MATLAB 编写程序。最终报告中附代码及实验结果截图。 (2) 实验报告中必须有对实验结果的分析讨论。 提示： (1) 可直接使用matlab 中已有函数产生高斯白噪声随机序列。可使用hist 函数画出序列的 直方图，并与标准高斯分布的概率密度函数做对比。 (2) 为便于卷积操作，当N 很大时，可近似认为h(N)=0。卷积使用matlab 自带的conv 函 数。 (3) 分析均值、方差等时，均可使用matlab 现有函数。功率谱密度和自相关函数可通过傅 里叶变换相互获得。傅里叶变换使用matlab 自带的fft 函数。 (4) 作图使用plot 函数。
一、作业分析： 本题主要考察的是加性高斯白噪声相关问题，因此构造一个高斯白噪声十分重要，故在本题中使用randn函数随机生成一个个符合高斯分布的数据，并由此构成高斯白噪声；而且由于白噪声是无法完全表示的，故此根据噪声长度远大于信号长度时可视为高斯白噪声，构造了一个长度为2000的高斯白噪声来进行试验。 二、作业解答： （1）matlab程序为： x-1000:1:1000; k=1*randn(1,length(x));% 生成零均值单位方差的高斯白噪声。 [f,xi]=ksdensity(x);%利用ksdensity函数估计样本的概率密度。 subplot(1,2,1); plot(x,k); subplot(1,2,2); plot(xi,f); 实验结果为：
西电人工智能大作业
八数码难题 一.实验目的 八数码难题：在3×3的方格棋盘上，摆放着1到8这八个数码，有1个方格是空的，其初始状态如图1所示，要求对空格执行空格左移、空格右移、空格上移和空格下移这四个操作使得棋盘从初始状态到目标状态。例如： (a) 初始状态 (b) 目标状态 图1 八数码问题示意图 请任选一种盲目搜索算法（深度优先搜索或宽度优先搜索）或任选一种启发式搜索方法（A 算法或 A* 算法）编程求解八数码问题（初始状态任选），并对实验结果进行分析，得出合理的结论。 本实验选择宽度优先搜索：选择一个起点，以接近起始点的程度依次扩展节点，逐层搜索，再对下一层节点搜索之前，必先搜索完本层节点。 二.实验设备及软件环境 Microsoft Visual C++，（简称Visual C++、MSVC、VC++或VC）微软公司的C++开发工具，具有集成开发环境，可提供编辑C语言，C++以及C++/CLI 等编程语言。 三.实验方法 算法描述： （1）将起始点放到OPEN表； （2）若OPEN空，无解，失败；否则继续； （3）把第一个点从OPEN移出，放到CLOSE表； （4）拓展节点，若无后继结点，转（2）； （5）把n的所有后继结点放到OPEN末端，提供从后继结点回到n的指针； （6）若n任意后继结点是目标节点，成功，输出；否则转（2）。
流程图：
代码： #include
学校:北京联合大学 系别:信息管理系 姓名：孙超 学号：2013110444006 《餐饮业信息管理系统的开发》 1、本项目的需求分析 随着今年来中国餐饮行业的日益火爆，在强烈的行业竞争中，一个高效的餐饮信息管理系统的应用，无疑是至关重要的。高效，便捷的管理系统，不仅仅极大的方便了食客的就餐，同时对于餐饮公司的各项信息管理有着很大的帮助，同时，我们的餐饮信息管理系统还能帮助餐厅降低错误率，扩大营业范围，增加知名度等。 为了使得系统在操作的过程中，更加便捷，具有针对性，本次系统设计主要分为：员工登陆操作信息系统，以及店主操作管理信息系统。不同的设计从而达到不同的功能，实现信息的有效传达与管理。 第一：在员工使用本餐饮信息管理系统应可以实现以下功能： 1.添加修改查询客户会员信息(修改客户信息需客户确认) 2.查询菜单 3.添加查询预定信息，为老顾客打折 4.客户可以在自己的会员账户里充值 5.顾客可以用现金买单也可以从会员账户里扣取 第二：管理员使用本餐饮信息管理系统应可以实现以下功能： 1.添加修改查询客户会员信息(修改客户信息需客户确认) 2.添加修改查询菜单信息，最好能看到菜品图片 3.添加查询预定信息，为老顾客打折 4.客户可以在自己的会员账户里充值 5.顾客可以用现金买单也可以从会员账户里扣取 6.设定具体的打折方法 7.添加职员信息，权限也可以定为管理员。 8.可以查询使用者的现金收款金额。 二、餐饮业管理数据库管理系统的E-R模型（概念结构设计） 1.用户（员工）的信息：
编号、密码、类型、姓名、电话、收款金额 2.客户信息： 用户编号、客户编号、姓名、电话、密码、开卡时间、卡内余额 3.食谱： 类型、名称、价格、配料、照片 4.预定： 用户编号、日期、预定时间、客户姓名、类型、预定食谱、桌号5桌台管理： 桌号、使用情况、 6.点餐管理： 用户编号、类型、菜品、数量、价格、照片 7.盈利管理： 日期、日支出金额、店内收入、外卖收入、盈利额度 各对象之间的联系图： 用户E-R图 主要存储一些用户信息，如用户的账号、密码和类型地点等等，主要用于用户登录，添加客户和添加预定时会使用到用户信息。
数字信号处理上机实验报告 14020710021 张吉凯 第一次上机 实验一： 设给定模拟信号()1000t a x t e -=，t 的单位是ms 。 (1) 利用MATLAB 绘制出其时域波形和频谱图（傅里叶变换），估计其等效带宽（忽略谱分量降低到峰值的3%以下的频谱）。 (2) 用两个不同的采样频率对给定的()a x t 进行采样。 ○1()()15000s a f x t x n =以样本秒采样得到。 ()()11j x n X e ω画出及其频谱。 ○2()()11000s a f x t x n =以样本秒采样得到。 ()() 11j x n X e ω画出及其频谱。 比较两种采样率下的信号频谱，并解释。 （1）MATLAB 程序： N=10; Fs=5; T s=1/Fs; n=[-N:T s:N]; xn=exp(-abs(n)); w=-4*pi:0.01:4*pi; X=xn*exp(-j*(n'*w)); subplot(211) plot(n,xn); title('x_a(t)时域波形'); xlabel('t/ms');ylabel('x_a(t)'); axis([-10, 10, 0, 1]); subplot(212); plot(w/pi,abs(X)); title('x_a(t)频谱图'); xlabel('\omega/\pi');ylabel('X_a(e^(j\omega))');
ind = find(X >=0.03*max(X))*0.01; eband = (max(ind) -min(ind)); fprintf('等效带宽为%fKHZ\n',eband); 运行结果： 等效带宽为12.110000KHZ
题目：数据结构上机报告学院：电子工程学院 专业：信息对抗技术 学生姓名：甘佳霖 学号：14020310092
西安电子科技大学 数据结构课程实验报告实验名称线性表 电子工程学院 1402031 班Array姓名甘佳霖学号 14020310092 同作者 实验日期 2017 年 3 月 18 日
实验一线性表 一、实验目的 1．熟悉线性表的顺序和链式存储结构 2．掌握线性表的基本运算 3．能够利用线性表的基本运算完成线性表应用的运算 二、实验要求 1.设有一个线性表E={e1, e2, … , e n-1, e n}，设计一个算法，将线性表逆置，即使元素排列次序颠倒过来，成为逆线性表E’={ e n, e n-1 , … , e2 , e1 }，要求逆线性表占用原线性表空间，并且用顺序表和单链表两种方法表示，分别用两个程序来完成。 2.已知由不具有头结点的单链表表示的线性表中，含有三类字符的数据元素（字母、数字和其他字符），试编写算法构造三个以循环链表表示的线性表，使每个表中只含有同一类的字符，且利用原表中的结点空间，头结点可另辟空间。 三、设计思路 1.顺序表做逆置操作时将对应的首尾元素位置交换，单链表的指针end指向链表的末尾，指针start指向链表头结点，指针s用来找到指向end节点的节点，将指向链表末尾和头结点的存储内容交换，然后头结点指针指向下一节点，s指针从start节点开始遍历寻找指向end 指针的节点，并将end指针赋值为s指针，就完成了单链表的逆置，可以看出单链表和顺序表都可以完成线性表的逆置。 2.分解单链表的实现思路是首先新建3个循环链表，然后顺序遍历单链表，ASCII码判断链表中的元素属于哪一类元素，然后将这个元素添加到对应的循环链表中，从而实现分解单链表的功能。 四、运行结果 1.单链表逆置：
《数字信号处理》课程教学大纲 课程代码：IB3123008 课程名称：数字信号处理英文名称：Digital Signal Processing 开课学期：第6学期 学分：3 学时：48 课程类別：必修课，专业基础课 适用专业：电子信息工程、信息对抗技术、遥感科学与技术、电磁场与无线技术、智能科学与技术 开课对象：三年级本科生 先修课程：信号与系统、MATLAB语言 后修课程：雷达原理、数字图像处理、数字音视频处理等 开课单位：电子工程学院 团队负责人：史林责任教授：史林 执笔人：史林核准院长：苏涛 一、课程性质、目的和任务 数字信号处理采用数字技术，研究信号和系统分析、处理、设计的基本原理和方法，是电子信息与电气工程类专业（电子信息工程专业、通信工程专业、信息工程专业等）的专业基础课，具有理论与实践紧密结合的特点。 通过本课程的学习，使学生建立数字信号处理的基本概念，掌握数字信号处理的基本原理、理论和方法，了解数字信号处理的新方法和新技术，熟练应用现代工具进行数字信号处理的仿真、分析和设计，达到能够对数字信号和系统进行分析、处理和设计的能力水平。为学习后续专业课程、进行创新性研究和解决复杂工程问题，奠定坚实的专业基础理论知识和工程实践能力。 本课程对学生达到如下毕业要求有贡献
二、教学内容、基本要求及学时分配 《数字信号处理》课程的教学内容、基本要求、学时分配和毕业要求指标点在教学中的具体体现如下。 （一）绪论 ( 2学时) 1.教学内容 介绍数字信号处理的基本概念、研究的内容及应用领域、发展概况和发展趋势，数字信号处理的基本特点，用数字方法处理信号的基本概念和一般方法。 2.基本要求 （1）了解数字信号处理研究的内容、应用领域、发展概况和发展趋势； （2）熟悉数字信号处理的基本概念和特点； （3）掌握用数字方法处理信号的基本概念和一般方法。 3.重点、难点 重点：数字信号处理的基本概念和特点。 难点：用数字方法处理信号的基本概念和一般方法 4.作业及课外学习要求 作业：分析数字信号处理的特点；熟悉用数字方法处理信号的一般方法，理解其中每个模块单元的作用。 课外学习：学习或复习MATLAB语言，掌握编程方法和技巧，做好后续的上机实验准备。 5.对毕业要求指标点的具体贡献 对指标点2-1的具体贡献：理解复杂工程问题中的数字系统； （二）离散时间信号和系统的时域分析( 4学时+4学时上机) 1.教学内容
(完整)2013西电射频大作业 编辑整理： 尊敬的读者朋友们： 这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（(完整)2013西电射频大作业）的内容能够给您的工作和学习带来便利。同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。 本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为(完整)2013西电射频大作业的全部内容。
班级 021171 射频电路基础期中大作业 学院电子工程学院 专业遥感科学与技术 学生姓名（02117017） (02117019） （02117024） 老师姓名朱天桥老师
题目一：基于 PSpice 仿真的振幅调制电路设计 1。1 基本要求 参考教材《射频电路基础》第五章振幅调制与解调中有关差分对放大器调幅和二极管调幅的原理，选择元器件、调制信号和载波参数，完成PSpice电路设计、建模和仿真,实现振幅调制信号的输出和分析。 1。2 实践任务 (1) 选择合适的调制信号和载波的振幅、频率，通过理论计算分析，正确选择晶体管和其它元件；搭建单端输出的差分对放大器，实现载波作为差模输入电压，调制信号控制电流源情况下的振幅调制;调整二者振幅,实现基本无失真的线性时变电路调幅；观察记录电路参数、调制信号、载波和已调波的波形和频谱。 （2）参考例5。3。1，修改电路为双端输出，对比研究平衡对消技术在该电路中的应用效果。 （3) 选择合适的调制信号和载波的振幅、频率，通过理论计算分析,正确选择二极管和其它元件；搭建单二极管振幅调制电路，实现载波作为大信号，调制信号为小信号情况下的振幅调制；调整二者振幅，实现基本无失真的线性时变电路调幅;观察记录电路参数、调制信号、载波和已调波的波形和频谱。 (4) 参考例5。3。2，修改电路为双回路,对比研究平衡对消技术在该电路中的应用效果. 摘要 此次振幅调制电路仿真设计基于PSpice,利用其丰富的仿真元器件库和强大的行为建模
西电排队论大作业 HEN system office room 【HEN16H-HENS2AHENS8Q8-HENH1688】
西安电子科技大学 (2016年度) 随机过程与排队论 班级： XXXXXXX 姓名： XXX XXX 学号： XXXXXXXXXX XXXXXXXXXXX 一步转移概率矩阵收敛快慢的影响因素 作者姓名：XXX XXX 指导老师姓名：XXX (西安电子科技大学计算机学院，陕西西安) 摘要：根据课程教材《排队现象的建模、解析与模拟【西安电子科技大学出版 社曾勇版】》，第[马尔可夫过程]中，马尔可夫过程链n时刻的k步转移概率结 果，当k=1时，得到一步转移概率。进而得到一步转移概率矩阵P（1）。为研究 此一步转移概率矩阵（下称一步矩阵）的收敛特性以及影响其收敛快慢的因素，使 用MATLAB实验工具进行仿真，先从特殊矩阵开始做起，发现规律，然后向普通矩 阵进行拓展猜想，并根据算术理论分析进行论证，最终得出一步矩阵收敛快慢的影 响因素。 关键词：一步转移概率矩阵 MATLAB 仿真猜想 一、问题概述 我们讨论时一步矩阵的特性应从以下两方面来分析： （1）矩阵P（n）在满足什么条件时具有收敛特性; 对于矩阵P（n）,当P（n）=P(n+1)时，我们说此矩阵 具有收敛特性，简称矩阵 P（n）收敛。 （2）若一个一步矩阵具有收敛特性，那么其收敛速度与什么有关
首先，我们需要明确什么是一步矩阵收敛： 对于一般的一步矩阵P 、矩阵An+1、矩阵An，若有： An+1=AnP=An 那么称该一步转移矩阵可收敛。 二、仿真实验 1、仿真环境 本次采用的是MATLAB仿真实验软件进行仿真实验 2、结果与分析 【1】、特殊矩阵：单位矩阵与类单位矩阵 从图（1）和图（2）可以看出，单位矩阵不具有收敛特性，类单位矩阵并非单位矩阵但是经过n次后也变为单位矩阵，所以此矩阵也不具有收敛特性。此类矩阵也易证明其不具有收敛性。 图（1）单位矩阵图（2）：类单位 矩阵 【2】、一般单位矩阵 图（3）：一般一步矩阵Ⅰ 图（4）：一般一步矩阵 从图（3）和（）可以看出他们分别在18次和4次后收敛到一个稳定的值 3、根据实验的猜想 根据在单位矩阵和一般单位矩阵和一般一步矩阵中得到的结果，可以对得出如下结论：类单位矩阵、单位矩阵是不具有收敛性的，而一般的一步矩阵是有收敛性的，而且收敛速率有快有慢。 对于上面结论中的状况，我们首先观察如上四个矩阵，不难发现，在矩阵收敛的最终结果矩阵中，其每行和均为1，而且每列上的值均为相同值。最终概率分布结果也是矩阵收敛后的一行。 所以根据上述的结果及分析做出如下猜想： 每一列比较均匀的矩阵收敛速度较快；与类单位矩阵类似的矩阵收敛速度较慢。 在极限情况下，有如下情况：
智能控制与应用实验报告神经网络控制器设计
一、 实验内容 考虑一个单连杆机器人控制系统，其可以描述为： 0.5sin()Mq mgl q y q τ+== 其中20.5M kgm =为杆的转动惯量，1m kg =为杆的质量，1l m =为杆长， 29.8/g m s =，q 为杆的角位置，q 为杆的角速度，q 为杆的角加速度， τ为系统的控制输入。具体要求： 1、设计神经网络控制器，对期望角度进行跟踪。 2、分析神经网络层数和神经元个数对控制性能的影响。 3、分析系统在神经网络控制和PID 控制作用下的抗干扰能力(加噪声干扰、加参数不确定)、抗非线性能力(加死区和饱和特性)、抗时滞的能力(对时滞大小加以改变)。 4、为系统设计神经网络PID 控制器(选作)。 二、 对象模型建立 根据公式(1)，令状态量121=,x q x x = 得到系统状态方程为： 12121 0.5**sin() x x mgl x x M y x τ=-= = (1) 由此建立单连杆机器人的模型如图1所示。
图1 单连杆机器人模型 三、系统结构搭建及神经网络训练 1.系统PID结构如图2所示: 图2 系统PID结构图 PID参数设置为Kp=16，Ki=10，Kd=8得到响应曲线如图3所示：
01234 5678910 0.2 0.4 0.6 0.8 1 1.2 1.4 t/s a n g l e /r a d 图3 PID 控制响应曲线 采样PID 控制器的输入和输出进行神经网络训练 p=[a1';a2';a3']; t=b'; net=newff([-1 1;-1 1;-1 1],[3 8 16 8 1],{'tansig' 'tansig' 'tansig' 'logsig' 'purelin'}); net.trainparam.epochs=2500; net.trainparam.goal=0.00001; net=train(net,p,t); gensim(net,-1) 产生的神经网络控制器如图4所示：
云南大学软件学院实验报告 课程：数据库原理与实用技术实验学期：任课教师： 专业：学号：姓名：成绩： 期末大作业：Electronic Ventor 数据库设计 一、实验目的 （1）掌握数据库设计的基本方法 （2）掌握各种数据库对象的设计方法 （3）熟练掌握DBA必须具备的技能 二、实验内容 1、根据项目的应用和项目的需求说明文档，进行详细的需求分析，给出需求分析的结果。 （1）客户可以在网站上注册，注册的客户要提供客户的姓名、电话、地址，以方便售后和联系，姓名即作为用户名，和密码一起用于注册和登录，客户编号可唯一识别用户，卡号可网上支付。其中地址、电话以方便联系和寄货； （2）网站管理员可以登记各种商品，供客户查询，订购。登记商品时要提供商品的名称、价格，商店中现有商品量，商品编号可唯一识别商品； （3）类别表示商品所属类别，类别编号可唯一识别类别，其中包含了，商品类别名称和制造厂商，可以对商品进行分类售卖； （4）客户可以在网上下订单，也可以到实体店购物，其在订单上所选择的支付方式不同（信用卡、借记卡、现金，现金代表实体店购物），网站管理员可以查看订单，并及时将订单的处理情况更新（比如货物已寄出的信息，订单状态：0：未处理，1：已处理，2：已发货）；订单编号可唯一识别订单，订单中包含订单产生时间，订单状态，支付方式和支付总额； （5）实体商店有自己的店名，卖多种商品，每个商店都有固定的地址，顾客可以到店中买商品，（注：在实体店中购买商品的顾客一律将顾客名默认为佚名），当商店中的库存量小于10时会有提醒到仓库中拿货； （6）配送单中包含查询号可唯一识别配送单，配送人，联系方式； （7）仓库中仓库编号可唯一识别仓库，其中每个仓库都有区号，代表其地址。 （8）各实体间关系 1）一个客户可以购买多种商品，一种商品可以被多个客户购买； 2）一个商品属于且仅属于一种类别，一种类别的商品可以包含多个商品或没有； 3）一种商品放在多个商店中销售，一个商店至少销售一种或销售多种商品； 4）一个订单对应一个客户，一个客户对应多个订单； 5）一个订单对应至少有一件商品或多件，一个商品对应多个订单； 6）一个订单可以有一个商品配送单 7）一个仓库可以存放多种商品，一种商品可以存放在一个仓库；
实验三、信号的频域与Z域分析 班级：学号：姓名：成绩： 1实验目的 （1）理解序列离散傅里叶变换（DTFT）的定义，熟悉序列DTFT的计算及其主要性质； （2）掌握Z变换的计算和主要性质，熟悉Z变换的收敛域及其与序列特性的关系，以及Z变换与DTFT的关系； （3）掌握时域离散线性时不变系统的频域分析方法，深刻理解系统的频率响应。了解系统的稳态响应和暂态响应、相位延迟和群延迟等概念； （4）掌握时域离散线性时不变系统的z域分析方法，深刻理解离散系统的系统函数及其零极点分布，熟悉零极点分布与系统的因果性和稳定性关系、零极点分布对系统频率特性的影响、差分方程的Z变换解法等； 2 实验内容 （1）设计计算机程序，产生序列并计算序列的DTFT，绘制其幅频特性和相频特性曲线； （2）根据系统的单位脉冲响应和差分方程，计算系统的频率响应，绘制系统频率响应的幅频特性和相频特性曲线； （3）根据系统的单位脉冲响应和差分方程，计算系统的系统函数、零极点分布；改变系统的零极点分布，观察系统频率响应的变化。 3实验步骤 （1）设计有限长序列Rn；计算序列的DTFT，绘制幅频特性和相频特性曲线 （2）改变系统的系统函数的零点分布，绘制系统改变前和改变后的频率响应的幅频特性和相频特性曲线 4 程序设计 x=[1,1,1,1];nx=[0:3];%x(n)=R(n) w=linspace(-2.8*pi,2.8*pi,100000);%取100000个点
X=x*exp(-j*nx'*w);%DTFT figure(1); subplot(3,2,1),plot(w/pi,abs(X));xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|') subplot(3,2,2),plot(w/pi,angle(X));xlabel('\omega/\pi');ylabel('\phi(\omega)/\pi') %差分方程求解 a=[1,-0.4];b=[1]; [H,w]=freqz(b,a,'whole'); subplot(3,2,3),plot(w/pi,abs(H));xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|') subplot(3,2,4),plot(w/pi,angle(H));xlabel('\omega/\pi');ylabel('\phi(\omega)/\pi') %零极点分布 a=[1,-1.6,0.9425];%分母 b1=[1,-0.3];b2=[1,-0.8];%分子 [F,w]=freqz(b1,a,'whole'); figure(2); subplot(2,2,1),zplane(b1,a);%零极点分布图 subplot(2,2,3),plot(w/pi,abs(F));xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|') subplot(2,2,4),plot(w/pi,angle(F));xlabel('\omega/\pi');ylabel('\phi(\omega)/\pi') figure(3);%改变零极点分布，观察频率响应变化 [F,w]=freqz(b2,a,'whole'); subplot(2,2,1),zplane(b2,a); subplot(2,2,3),plot(w/pi,abs(F));xlabel('\omega/\pi');ylabel('|X(e^j^\omega)|') subplot(2,2,4),plot(w/pi,angle(F));xlabel('\omega/\pi');ylabel('\phi(\omega)/\pi')
随机信号分析 实验报告 目录 随机信号分析 (1) 实验报告 (1) 理想白噪声和带限白噪声的产生与测试 (2) 一、摘要 (2) 二、实验的背景与目的 (2) 背景： (2) 实验目的: (2) 三、实验原理 (3) 四、实验的设计与结果 (4) 实验设计: (4) 实验结果: (5) 五、实验结论 (12) 六、参考文献 (13) 七、附件 (13) 1
理想白噪声和带限白噪声的产生与测试一、摘要 本文通过利用MATLAB软件仿真来对理想白噪声和带限白噪声进行研究。理想白噪声通过低通滤波器和带通滤波器分别得到低通带限白噪声和帯通带限白噪声。在仿真的过程中我们利用MATLAB工具箱中自带的一些函数来对理想白噪声和带限白噪声的均值、均方值、方差、功率谱密度、自相关函数、频谱以及概率密度进行研究，对对它们进行比较分析并讨论其物理意义。 关键词：理想白噪声带限白噪声均值均方值方差功率谱密度自相关函数、频谱以及概率密度 二、实验的背景与目的 背景： 在词典中噪声有两种定义：定义1：干扰人们休息、学习和工作的声音，引起人的心理和生理变化。定义2：不同频率、不同强度无规则地组合在一起的声音。如电噪声、机械噪声，可引伸为任何不希望有的干扰。第一种定义是人们在日常生活中可以感知的，从感性上很容易理解。而第二种定义则相对抽象一些，大部分应用于机械工程当中。在这一学期的好几门课程中我们都从不同的方面接触到噪声，如何的利用噪声，把噪声的危害减到最小是一个很热门的话题。为了加深对噪声的认识与了解，为后面的学习与工作做准备，我们对噪声进行了一些研究与测试。 实验目的: 了解理想白噪声和带限白噪声的基本概念并能够区分它们，掌握用MATLAB 或c/c++软件仿真和分析理想白噪声和带限白噪声的方法，掌握理想白噪声和带限白噪声的性质。