信号与系统上机实验报告

信号与系统上机实验报告学院：电气与信息工学院姓名：班级：学号：

实验目的验证采样定理会进行傅里叶分析，分析信号S（t）的频谱组成信号处理的基本过程（调制、采样、滤波）数字滤波器的设计与应用二.实验内容给定如下信号:QUOTE，QUOTEccQUOTEntnt为均值为0，均方差为5的高斯白噪声1)、编写MATLAB程序分析如下过程：使用不同的采样率对QUOTEs(t)s(t)进行采样，选用两组数字滤波器恢复QUOTEb(t)b(t)和QUOTEc2(t)c2(t)，画出QUOTEs(t)s(t)产生前后，采样后、滤波后的各阶段波形，并对现象进行分析；(2)、应用SIMULINK对上述信号调制、采样、滤波过程进行仿真。三．仿真分析1.代码如下clc,clear;T=10;%假定时宽为10fs=600;%采样频率t=0:1/fs:T;%时间;tN=length(t);%tN为50001，即t向量的max{行，列}f=0:(2*pi/tN):(2*pi);%频率f=f(1:tN);%f这个向量的第1个元素到tN个元素c1=sin(200*pi*t);figure(1);%c1tsubplot(2,1,1);plot(t(1:100),c1(1:100));subplot(2,1,2);plot(f,abs(fft(c1)));c2=2*sin(400*pi*t);figure(2);%c2tsubplot(2,1,1);plot(t(1:100),c2(1:100));subplot(2,1,2);plot(f,abs(fft(c2)));bt=zeros(1,tN);%表示一个1行tN列的0矩阵tM=0.1*fs;fori=0:(T/1-1);%bt((i*fs+1):(i*fs+tM))=ones(1,tM);endfigure(3);%btsubplot(2,1,1);plot(t(1:100),bt(1:100));subplot(2,1,2);plot(f,abs(fft(bt)));nt=wgn(1,tN,0.5,'linear');figure(4);%whitenoisesubplot(2,1,1);plot(t(1:100),nt(1:100));subplot(2,1,2);plot(f,abs(fft(nt)));st=bt.*c1+c2+nt;figure(5);%stsubplot(2,1,1);plot(t(1:100),st(1:100));subplot(2,1,2);plot(f,abs(fft(st)));bandfilter1=designfilt('bandpassfir','FilterOrder',20,...'CutoffFrequency1',190,'CutoffFrequency2',210,...'SampleRate',10000);c1after=filter(bandfilter1,st);figure(6);%c2_afterthroughbandpassfiltersubplot(2,1,1);plot(t(1:100),c1after(1:100));subplot(2,1,2);plot(f,abs(fft(c1after)));st1=st.*c1;figure(7);%st_aftersubplot(2,1,1);plot(t(1:100),st1(1:100));subplot(2,1,2);plot(f,abs(fft(st1)));bandfilter2=designfilt('lowpassfir','PassbandFrequency',40/10000,...'StopbandFrequency',50/10000,'PassbandRipple',0.2,...'StopbandAttenuation',65,'DesignMethod','kaiserwin');btafter=filter(bandfilter2,st1);figure(8);%bt_afterthroughlowpassfiltersubplot(2,1,1);plot(t(1:100),btafter(1:100));subplot(2,1,2);plot(f,abs(fft(btafter)));具体分析分析过程：为了方便采样，我们选取了时宽t=10的波形，接着我们用频率fs为500Hz对波形进行了采样。得到下图所示波形：通过波形图我们发现这波形与我们所想得到的正弦波相差甚远，经过我们分析发现这由于我们的取样频率较低（满足采样定律），使得采样之后的波形不够明显，为此我们决定采用更大的fs来进行采样，所以选择了fs=10000Hz来进行采样采样之后的波形如下：由于采样频率的提高，采样的波形增多，密度更大，在时域上所呈的取样波形十分接近正弦波，同时由于冲激信号的抽样特性，在频域上有两个十分明显的冲激信号，为此我们决定用频率fs为10000Hz的进行接下来的实验。为此c2t在时域上和s域上的波形如下所示：时域上有一个多周期的波形，而在频域上，根据sinwt的傅里叶变换可知在400pi和-400pi的位置上各有一个冲激信号，如图中所示。而bt在时域上和s域上的波形图就如下图所示：由于fs过大导致bt这个门函数在时域上的小区间内显得很大且只有一个门，根据矩形脉冲f(t)=E[u(t+τ/2)-u(t-τ/2)]的傅里叶变换为：F（w）=EτSa(wτ/2)可知在wτ/2处有一个幅度为0.1的冲激信号，由于该函数为连续函数，为此有无数的冲激，而根据此时所取的频率，图像上有两个冲激。根据题目要求n(t)s是均值为0，均方差为5的高斯白噪声，为此在时域上和频域上的波形如下图所示：现在我们已经得到了b(t),c1(t),c2(t)和n(t)波形图，为此我们得到s(t)在时域上和频域上的波形图:杨浦两个杨浦两个从频域上我们发现有两个干扰信号产生，为了得到c2t的波形，我们决定采用带通滤波器进行滤波消除干扰信号，滤波之后的图形如下：为此我们得到了c2t原来的图形同理我们也得到了c1t的波形图：为了得到bt的图形我们采用低通滤波器将其过滤出来，得到了如下：至此我们基本完成了采样，滤波，恢复等过程，得到了个阶段的波形图。

2Simulink仿真图源为两个正弦信号c1，c2，方波信号b(t)和噪声信号。实现对s(t)的采样分析和b(t),c2de恢复。低通滤波器设置：带通滤波器设置S(t)的频谱分析图b（t）与c2（t）通过滤波器后示波器上显示的图形。b(t)为上图，c2(t)为下图。

公司印章管理制度一、目的公司印章是公司对内对外行使权力的标志，也是公司名称的法律体现，因此，必须对印章进行规范化、合理化的严格管理，以保证公司各项业务的正常运作，由公司指定专人负责管理。二、印章的种类公章，是按照政府规定，由主管部门批准刻制的代表公司权力的印章。专用章，为方便工作专门刻制的用于某种特定用途的印章，如：合同专用章、财务专用章、业务专用章、仓库签收章等。3、手章(签名章)，是以公司法人代表名字刻制的用于公务的印章。三、印章的管理规定印章指定专人负责保管和使用，保管印章的地方(桌、柜等)要牢固加锁，印章使用后要及时收存。财务专用章由财务部负责保管，向银行备案的印章，应由财务部会计、总经办分别保管。3、印章要注意保养，防止碰撞，还要及时清洗，以保持印迹清晰。4、一般情况下不得将印章携出公司外使用，如确实因工作所需，则应由印章管理员携带印章到场盖章或监印。5、印章管理人员离职或调任时，须履行印章交接手续。四、公章刻制印章需本公司法人代表批准，并由印章管理专责人负责办理刻制并启用并交由专人进行保管。五、印章的使用使用任何的印章，需由相应负责人审核签字。为方便工作，总经理可授权印章管理专责人审核一般性事务用印。用印前印章管理人员须认真审核，明确了解用印的内容和目的，确认符合用印的手续后，在用印登记簿上逐项登记，方可盖章。3、对需要留存的材料，盖印后应留存一份立卷归档。4、不得在空白凭证、便笺上盖章。5、上报有关部门的文件资料，未经部门经理、总经理审签，不得盖章。6、以公司名义行文，未经总经理签发，不得盖章。7、按照合同会签制度的规定，所有合同和