数字信号处理课程设计离散时间信号处理
- 格式:doc
- 大小:670.96 KB
- 文档页数:14
目录1、软件介绍2、设计原理3、设计过程3.1设计内容3.2设计步骤4、程序运行结果与总结5、心得体会6、参考文献1、软件介绍MATLAB是一种科学计算软件,适用于工程应用个领域的分析设计与复杂计算,它使用方便,输入简捷,运算高效且内容丰富,很容易由用户自行扩展。
因此,当前已成为美国和其他发达国家大学教学和科学研究中最常用且必不可少的工具。
MATLAB以矩阵作为基本编程单元,它提供了各种矩阵的运算与操作,并有较强的绘图功能。
MATLAB集科学计算、图像处理、声音处理于一身,是一个高度的集成系统,有良好的用户界面,并有良好的帮助功能。
MATLAB不仅流行于控制界,在机械工程、生物工程、语音处理、图像处理、信号分析、计算机技术等各行各业中都有极广泛的应用。
MATLAB语言的特点:1.编程效率高;2.用户使用方便;3.扩充能力强;4.语句简单,内涵丰富;5.高效方便的矩阵和数组运算;6.方便的绘图功能。
数字信号处理的主要对象是数字信号,且是采用运算的方法达到处理目的,因此,其实现方法不同于模拟信号的实现方法,基本上可以分为两种实现方法,即软件实现方法和硬件实现的方法。
而硬件实现指的是选用合适的DSP芯片,配有适合芯片语言及任务要求的软件,实现某种信号处理功能的一种方法。
数字信号处理的特点:1.灵活性;2.高精度和高稳定性;3.便于大规模集成;4.对数字信号可以存储、运算;5系统可以获得高性能指标。
2、设计原理1、常用典型序列(1)单位阶跃序列单位阶跃序列)(n u 定义为在21n n n ≤≤区间的值,MATLAB 中可用下列表达式描述: n = [ n 1 : n 2 ]; x = [ (n - n 0) > =0 ];(2)实数指数序列na n x n∀= ,)(例如要实现序列100 ,)9.0()(≤≤=n n x n,MATLAB 中可用下列表达式描述: n = [ 0 : 10 ]; x =(0.9 ).^ n ;(3)复数指数序列nen x nj ∀=+ ,)()(0ωσ例如要实现序列100 ,)()32(≤≤=+n en x nj ,MATLAB 中可用下列表达式描述:n = [ 0 : 10 ]; x = exp( (2+3j)*n ) ;(4)矩形序列矩形序列)(n R N 定义为),0()10(01)(N n n N n n R N≥<-≤≤⎩⎨⎧= (2-3)矩形序列有一个重要的参数,就是序列宽度N 。
离散时间信号处理课程设计某雷达系统接收机框架如图1所示。
接收机输入信号,其中,中心频率=40MHz ,为相位调制函数(调制带宽B=2MHz),信号谱如图2所示。
输出基带信号I[n]、Q[n]数据率(采样率)为2.5MHz 。
图1 接收机原理框架图2 接收信号参数1. 请设计ADC 的采样率fs ,画出x[n]的频谱X(e jw )2. 请推导正交解调器的输出S_I[n]和S_Q[n]的表达式,画出频谱示意图;若要将信号谱搬移到零中频,请确定NCO 的频率W 03. 请设计一个FIR 线性相位数字LPF 对正交解调的输出进行处理,要求杂散抑制比超过50dBc ,确定滤波器设计指标、给出设计过程和结果。
4. 请确定采样率变换模块的参数(抽取或内插系数),画出输出信号谱。
fXc(f)1.设计ADC的采样率fs,画出x[n]的频谱X(e jw)Adc采样率fs>=2*f0,fs取100MHZ。
φ(t)为相位调制函数(调制带宽B=2MHz),取φ(t)=cos(2π*1000000*t)作为单音调制信号,则,,,程序如下:clear;clc;clf;n0=20000;%点数sprate=100000000;%采样频率100Mt0=n0/sprate;%总时长fca=40000000;%载波频率fsi=1000000;%信号频率t=linspace(0,t0,n0);Dsb=cos(2*pi*fca*t+cos(fsi*2*pi*t));subplot(3,1,1);plot(t(1:100),Dsb(1:100));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);w=linspace(-1,1,n0);plot(w,fftshift(abs(fft(Dsb))));title('幅频');xlabel('w');ylabel('|X(ejw)|');grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Dsb))));title('相频');xlabel('w');ylabel('∠X(ejw)');grid on;legend('modulated signal');上图中频域图0.8代表0.8*100000000/2 = 40M2.推导正交解调器的输出S_I[n]和S_Q[n]的表达式,画出频谱示意图;若要将信号谱搬移到零中频,请确定NCO的频率WW0==4π/5S_I[n] = x[n]* cos[wn]=cos[4πn/5 +cos[πn/50]]*cos[4πn/5]Sl=Dsb.*cos(2*pi*fca*t);subplot(3,1,1);plot(t(1:100),Sl(1:100));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);plot(w,fftshift(abs(fft(Sl))));title('幅频');xlabel('w');ylabel('|X(ejw)|');grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sl))));title('相频');xlabel('w');ylabel('∠X(ejw)');grid on;legend('modulated signal');n]=cos[4πn/5 +cos[πn/50]]*sin[4πn/5] S_Q[n] =x[n]* sin[wSq=Dsb.*sin(2*pi*fca*t);subplot(3,1,1);plot(t(1:100),Sq(1:100));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);plot(w,fftshift(abs(fft(Sq))));title('幅频');xlabel('w');ylabel('|X(ejw)|');grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sq))));title('相频');xlabel('w');ylabel('∠X(ejw)');grid on;legend('modulated signal');3.设计一个FIR线性相位数字LPF对正交解调的输出进行处理,要求杂散抑制比超过50dBc,确定滤波器设计指标、给出设计过程和结果。
离散时间信号处理教学设计一、教学目标本课程的主要目标是让学生掌握离散时间信号的基本概念、特性、分类、表示以及常见离散时间信号处理方法。
具体目标如下:1.理解离散时间信号的基本概念和特性;2.掌握离散时间信号的分类及其表示方法;3.熟练掌握离散时间信号处理的常见方法及其应用。
二、教学内容本课程主要包括以下内容:1. 离散时间信号的基本概念和特性1.离散时间序列的概念;2.离散时间序列的特性:周期性、奇偶性、对称性、能量、功率等;3.离散时间信号的采样和量化。
2. 离散时间信号的分类及其表示方法1.常见离散时间信号的分类;2.离散时间信号的表示方法:采样序列、单位取样序列、单位脉冲序列、单位斜坡序列、正弦序列等。
3. 离散时间信号处理的常见方法及其应用1.离散时间信号处理的模型和抽象表示;2.离散时间信号的线性时不变系统;3.离散时间信号的卷积与相关;4.快速傅里叶变换在离散时间信号处理中的应用。
三、教学方法本课程主要采用课堂讲授与实例教学相结合的方式进行。
1. 课堂讲授本课程将采用讲授法,由教师系统讲解离散时间信号处理的相关概念、理论和方法。
2. 实例教学本课程将通过实例讲解,将理论知识和应用实践相结合,以便学生更好地掌握离散时间信号处理的相关方法、技巧和应用。
四、教学评价为了确保学生掌握了离散时间信号处理的相关知识和方法,本课程将采用以下方式进行评价:1. 课堂参与度学生在课堂上的表现和参与度将作为评价因素之一。
2. 作业和考试学生将进行课堂作业和期末考试,以检验他们对离散时间信号处理相关知识和方法的掌握程度。
3. 综合成绩学生的综合成绩将由其在课堂参与度、作业和考试中所得分数计算而来。
五、教学资源1. 教材1.信号与系统分析;2.数字信号处理;3.数字信号处理的快速算法。
2. 工具软件1.Matlab;2.Python;六、教学时长本课程设置为总时长 48 学时,授课时间为每周 2 学时,共计 24 周。
一、实验目的1. 理解数字信号处理的基本概念和原理。
2. 掌握离散时间信号的基本运算和变换方法。
3. 熟悉数字滤波器的设计和实现。
4. 培养实验操作能力和数据分析能力。
二、实验原理数字信号处理(Digital Signal Processing,DSP)是利用计算机对信号进行采样、量化、处理和分析的一种技术。
本实验主要涉及以下内容:1. 离散时间信号:离散时间信号是指时间上离散的信号,通常用序列表示。
2. 离散时间系统的时域分析:分析离散时间系统的时域特性,如稳定性、因果性、线性等。
3. 离散时间信号的变换:包括离散时间傅里叶变换(DTFT)、离散傅里叶变换(DFT)和快速傅里叶变换(FFT)等。
4. 数字滤波器:设计、实现和分析数字滤波器,如低通、高通、带通、带阻滤波器等。
三、实验内容1. 离散时间信号的时域运算(1)实验目的:掌握离散时间信号的时域运算方法。
(2)实验步骤:a. 使用MATLAB生成两个离散时间信号;b. 进行时域运算,如加、减、乘、除等;c. 绘制运算结果的时域波形图。
2. 离散时间信号的变换(1)实验目的:掌握离散时间信号的变换方法。
(2)实验步骤:a. 使用MATLAB生成一个离散时间信号;b. 进行DTFT、DFT和FFT变换;c. 绘制变换结果的频域波形图。
3. 数字滤波器的设计和实现(1)实验目的:掌握数字滤波器的设计和实现方法。
(2)实验步骤:a. 设计一个低通滤波器,如巴特沃斯滤波器、切比雪夫滤波器等;b. 使用MATLAB实现滤波器;c. 使用MATLAB对滤波器进行时域和频域分析。
4. 数字滤波器的应用(1)实验目的:掌握数字滤波器的应用。
(2)实验步骤:a. 采集一段语音信号;b. 使用数字滤波器对语音信号进行降噪处理;c. 比较降噪前后的语音信号,分析滤波器的效果。
四、实验结果与分析1. 离散时间信号的时域运算实验结果显示,通过MATLAB可以方便地进行离散时间信号的时域运算,并绘制出运算结果的时域波形图。
离散时间信号处理课程设计1. 概述离散时间信号处理是数字信号处理的重要学科,它研究的是离散时间信号的获取、处理和分析,是现代通信、控制工程等领域中必不可少的一门技术。
在本次课程设计中,我们将学习如何使用 MATLAB 软件进行离散时间信号的处理,包括信号的采样、重构、滤波、时域与频域分析等内容。
2. 实验原理2.1 信号采样采样是将连续时间信号转化为离散时间信号的过程。
在实验中,我们将学习如何进行周期采样、非周期采样、抽样保持等方式的信号采样。
采样频率决定着信号的重构质量,通常情况下,采样频率应该远高于信号的最高频率分量。
2.2 信号重构信号重构是将离散时间信号转化为连续时间信号的过程。
在实验中,我们将学习如何进行零阶保持、一阶保持、插值、最小二乘法等方式的信号重构。
信号重构的质量取决于采样频率以及信号重构方法。
2.3 信号滤波信号滤波是指对离散时间信号进行滤波处理,以达到去除噪声、平滑信号、高通或低通滤波等目的。
在实验中,我们将学习 IIR 滤波器、FIR 滤波器等常见的数字信号滤波器的设计与应用。
2.4 时域分析时域分析是指对离散时间信号进行时间域分析,如均值、方差、相关系数、自相关函数等分析方法,以便于更好地理解信号特征。
在实验中,我们将学习MATLAB 软件中常用的时域分析函数。
2.5 频域分析频域分析是指对离散时间信号进行频域分析,如傅里叶变换、功率谱密度、自谱密度等分析方法,以便于更好地理解信号频率特征。
在实验中,我们将学习MATLAB 软件中常用的频域分析函数。
3. 实验内容3.1 实验一:信号采样与重构在实验一中,我们将学习如何进行周期采样、非周期采样、抽样保持等方式的信号采样,并对采样后的信号进行零阶、一阶保持、最小二乘法等方式的信号重构,最终对重构信号进行比较和评价。
3.2 实验二:数字信号滤波在实验二中,我们将学习 IIR 滤波器、FIR 滤波器的设计与应用,以及通过MATLAB 软件进行数字信号滤波的实现。
离散时间信号处理第三版课程设计一、课程设计概述本文介绍离散时间信号处理第三版课程设计内容,旨在掌握离散时间信号处理的基本概念、信号分析和滤波器设计技术等。
本课程设计包括以下内容:•信号取样与重构•时域信号分析•频域分析•离散滤波器设计•快速傅里叶变换学生需要使用Matlab软件完成相关实验,包括程序设计、仿真、算法实现和实验报告等。
二、实验设计2.1 信号取样与重构任务描述:对连续时间信号进行采样,并用重构公式将重建的信号与原始信号进行比较,分析采样频率对重构误差的影响。
实验步骤:1.设计一个连续时间信号,如三角波,正弦波等;2.对信号进行采样,采样频率分别为2倍、1.5倍、1倍、0.8倍信号带宽;3.使用重构公式重建信号,并与原始信号进行比较,分析采样频率对重构误差的影响;4.编写程序实现上述过程,并撰写实验报告。
2.2 时域信号分析任务描述:分析时域信号的基本特征,如幅值、相位、周期性等,并对信号进行傅里叶变换,分析频谱特性。
实验步骤:1.设计一个周期性矩形脉冲信号;2.使用Matlab自带的fft函数进行傅里叶变换,得到幅频响应;3.分析信号周期、幅值、相位以及频谱特性;4.仿真不同信号,分析时域波形与频域特征,并撰写实验报告。
2.3 频域分析任务描述:对采样信号进行傅里叶变换,并进行频域分析。
实验步骤:1.设计一个带通滤波器,设计通带边界、阻带边界和折返带宽等参数;2.对信号进行采样,并进行傅里叶变换;3.计算频率响应曲线和幅度、相位响应值;4.分析设计参数对滤波器性能的影响;5.编写程序实现上述过程,并撰写实验报告。
2.4 离散滤波器设计任务描述:对数字滤波器进行设计并进行性能测试。
实验步骤:1.设计一个数字低通滤波器;2.使用Butterworth滤波器设计方法进行设计;3.对设计参数进行分析和调整,寻找最优方案;4.使用Matlab实现滤波器,测试其性能;5.分析设计参数、阶数、截止频率等对滤波器性能的影响;6.编写程序实现上述过程,并撰写实验报告。
数字信号处理教程第五版教学设计课程简介本课程是数字信号处理教程的第五版,旨在教授数字信号处理的基础理论、算法和应用。
通过本课程的学习,学生将了解如何在数字领域中进行信号处理,包括滤波、采样、功率谱估计和谱分析等。
课程目标本课程的目标是:1.掌握数字信号处理的基本理论和概念。
2.熟悉数字信号处理中的常用算法和技术。
3.能够在实际应用中运用所学知识进行数字信号处理。
4.培养学生的理论研究和实践能力。
教学安排本课程将分为以下几个模块:第一模块:信号与系统基础本模块将讲解信号与系统的基础知识,包括信号的分类、信号的时域和频域表示、系统的线性性和时不变性等内容。
第二模块:离散信号与系统本模块将介绍离散信号和离散系统的基础知识,包括离散时间信号和连续时间信号的转换、离散时间系统和连续时间系统的转换、离散时间卷积和相关等内容。
第三模块:数字滤波器本模块将讲解数字滤波的基本概念和分类、实现数字滤波的不同方法、数字滤波器的设计和优化等内容。
第四模块:数字信号的采样和重构本模块将讲解数字信号的采样和重构,包括采样定理、插值和抽样等内容。
第五模块:功率谱估计和谱分析本模块将讲解数字信号的功率谱估计和谱分析,包括周期图和谱密度函数、特征值分解和Prony方法等内容。
教学方法本课程的教学方法主要采用讲授和实践相结合的方式。
其中,讲授部分将使用教材和辅助课件进行,包括教授基本概念、算法和应用;实践部分将开展编程实验和课程设计,引导学生解决实际问题,提高学生的独立思考和实践能力。
评估方式本课程的评估方式包括以下几个方面:1.考试成绩:占总成绩的50%。
2.实验成绩:占总成绩的30%。
3.课程设计成绩:占总成绩的20%。
总结本课程是一门重要的基础课程,对于数字信号处理及其应用的学习具有重要的意义。
希望学生通过本课程的学习,能够掌握数字信号处理的基本概念和方法,并在实际应用中发挥出所学知识的价值。
《数字信号处理》课程教案数字信号处理课程教案第一部分:课程概述数字信号处理是现代通信和信号处理领域中的重要学科,本课程旨在介绍数字信号处理的基本概念和理论,并探讨其在实际应用中的应用和技术。
第二部分:教学目标1. 理解数字信号处理的基本原理和基础知识;2. 掌握数字信号的采样、量化和编码技术;3. 了解常见的数字滤波器设计方法;4. 学习数字信号处理中的快速傅里叶变换(FFT)算法;5. 探讨数字信号处理在音频、图像和视频信号处理中的应用。
第三部分:教学内容1. 数字信号处理基础知识1.1 数字信号与模拟信号的比较1.2 采样和量化1.3 数字信号编码1.4 常见信号的时域和频域表示2. 离散时间信号和系统2.1 离散时间信号的表示和性质2.2 线性时不变系统2.3 离散时间系统的性质和分类3. 离散时间系统的频域分析3.1 离散时间信号的傅里叶变换3.2 离散频域系统的频率响应3.3 滤波器的设计和实现4. 数字滤波器设计4.1 IIR滤波器的设计方法4.2 FIR滤波器的设计方法4.3 改进的滤波器设计方法5. 快速傅里叶变换(FFT)算法5.1 傅里叶变换的基本概念及性质5.2 离散傅里叶变换(DFT)及其性质5.3 快速傅里叶变换算法及其应用6. 数字信号处理在多媒体中的应用6.1 音频信号处理技术6.2 图像信号处理技术6.3 视频信号处理技术第四部分:教学方法1. 理论讲授与案例分析相结合,通过实际应用案例来深化理解;2. 课堂互动,鼓励学生提问和参与讨论;3. 实验操作,通过实际操作提升学生的实践能力;4. 小组合作,鼓励学生进行小组项目研究和报告。
第五部分:教学评估1. 平时表现:出勤、课堂参与和作业完成情况;2. 期中考试:对课程前半部分内容的回顾和检验;3. 实验报告:根据实验内容,撰写实验报告并提交;4. 期末考试:综合检验对整个课程的掌握情况。
第六部分:教材与参考书目主教材:《数字信号处理导论》(第四版),作者:约翰·G·普罗阿基斯;参考书目:1. 《数字信号处理》(第四版),作者:阿兰·V·奥泽;2. 《数字信号处理:实用方法与应用》(第三版),作者:埃密里奥·马其尔夏兰德。
离散时间信号处理课程设计某雷达系统接收机框架如图1所示。
接收机输入信号,其中,中心频率=40MHz ,为相位调制函数(调制带宽B=2MHz),信号谱如图2所示。
输出基带信号I[n]、Q[n]数据率(采样率)为2.5MHz 。
图1 接收机原理框架图2 接收信号参数1. 请设计ADC 的采样率fs ,画出x[n]的频谱X(e jw )2. 请推导正交解调器的输出S_I[n]和S_Q[n]的表达式,画出频谱示意图;若要将信号谱搬移到零中频,请确定NCO 的频率W 03. 请设计一个FIR 线性相位数字LPF 对正交解调的输出进行处理,要求杂散抑制比超过50dBc ,确定滤波器设计指标、给出设计过程和结果。
4. 请确定采样率变换模块的参数(抽取或内插系数),画出输出信号谱。
fXc(f)1.设计ADC的采样率fs,画出x[n]的频谱X(e jw)Adc采样率fs>=2*f0,fs取100MHZ。
φ(t)为相位调制函数(调制带宽B=2MHz),取φ(t)=cos(2π*1000000*t)作为单音调制信号,则,,,程序如下:clear;clc;clf;n0=20000;%点数sprate=100000000;%采样频率100Mt0=n0/sprate;%总时长fca=40000000;%载波频率fsi=1000000;%信号频率t=linspace(0,t0,n0);Dsb=cos(2*pi*fca*t+cos(fsi*2*pi*t));subplot(3,1,1);plot(t(1:100),Dsb(1:100));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);w=linspace(-1,1,n0);plot(w,fftshift(abs(fft(Dsb))));title('幅频');xlabel('w');ylabel('|X(ejw)|');grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Dsb))));title('相频');xlabel('w');ylabel('∠X(ejw)');grid on;legend('modulated signal');上图中频域图0.8代表0.8*100000000/2 = 40M2.推导正交解调器的输出S_I[n]和S_Q[n]的表达式,画出频谱示意图;若要将信号谱搬移到零中频,请确定NCO的频率WW0==4π/5S_I[n] = x[n]* cos[wn]=cos[4πn/5 +cos[πn/50]]*cos[4πn/5]Sl=Dsb.*cos(2*pi*fca*t);subplot(3,1,1);plot(t(1:100),Sl(1:100));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);plot(w,fftshift(abs(fft(Sl))));title('幅频');xlabel('w');ylabel('|X(ejw)|');grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sl))));title('相频');xlabel('w');ylabel('∠X(ejw)');grid on;legend('modulated signal');n]=cos[4πn/5 +cos[πn/50]]*sin[4πn/5] S_Q[n] =x[n]* sin[wSq=Dsb.*sin(2*pi*fca*t);subplot(3,1,1);plot(t(1:100),Sq(1:100));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);plot(w,fftshift(abs(fft(Sq))));title('幅频');xlabel('w');ylabel('|X(ejw)|');grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sq))));title('相频');xlabel('w');ylabel('∠X(ejw)');grid on;legend('modulated signal');3.设计一个FIR线性相位数字LPF对正交解调的输出进行处理,要求杂散抑制比超过50dBc,确定滤波器设计指标、给出设计过程和结果。
在20Mhz左右范围内混叠产生一些频率需滤除,所以设计fir的指标为:fpass=5MHZ,fstop=15MHZ,要求杂散抑制比超过50dBc,取Astop=50db时就已足够了。
窗函数采用Kaiser窗,wp=π/10,ws=3π/10,wc=π/5,△w=π/5,A=50β=0.5842*(A-21)^0.4+0.07886*(A-21)=4.5336M=(A-8)/(2.285*pi/5) = 29.254M取30pp=linspace(0,29,30);hlp=(sin((pp-15)*pi/5))./((pp-15)*pi);hlp(16)=0.2;h=hlp.*FIRLPF';scatter(pp,h,'o');subplot(3,1,2);ll=linspace(0,2*pi,20000);h3=zeros(1,20000);for k=1:20000h3(k)=0;for i=1:30h3(k) = h3(k) + h(i)*exp(-1j*2*pi*(k-1)*(i-1)/20000);endendplot(ll,10*log(abs(h3)));grid on;xlabel('w');ylabel('dB');4.确定采样率变换模块的参数(抽取或内插系数),画出输出信号谱。
S_I[n]通过低通滤波器后Sl=Dsb.*cos(2*pi*fca*t);Sl=conv(Sl,h);subplot(3,1,1);plot(t(1:1000),Sl(1:1000));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);w=linspace(-1,1,20029);plot(w,fftshift(abs(fft(Sl))));title('幅频');xlabel('w');ylabel('|X(ejw)|');grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sl))));title('相频');xlabel('w');ylabel('∠X(ejw)');grid on;legend('modulated signal');采样率降为2.5M,即抽取系数为40,程序如下Sll=zeros(1,500);for i=1:500Sll(i)=Sl(i*40);endsubplot(3,1,1);plot(t(1:200),Sll(1:200));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);w=linspace(-1,1,500);plot(w,fftshift(abs(fft(Sll))));title('幅频');xlabel('w');ylabel('|X(ejw)|');grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sll)))); title('相频');xlabel('w');ylabel('∠X(ejw)'); grid on;legend('modulated signal');对S_Q做同样的操作Sq=Dsb.*sin(2*pi*fca*t);Sq=conv(Sq,h);subplot(3,1,1);plot(t(1:1000),Sq(1:1000));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);w=linspace(-1,1,20029);plot(w,fftshift(abs(fft(Sq)))); title('幅频');xlabel('w');ylabel('|X(ejw)|'); grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sq))));title('相频');xlabel('w');ylabel('∠X(ejw)'); grid on;legend('modulated signal');降采样后Sqq=zeros(1,500);for i=1:500Sqq(i)=Sq(i*40);endsubplot(3,1,1);plot(t(1:200),Sqq(1:200));title('时域');xlabel('t');ylabel('voltage');grid on;legend('modulated signal');subplot(3,1,2);w=linspace(-1,1,500);plot(w,fftshift(abs(fft(Sqq)))); title('幅频');xlabel('w');ylabel('|X(ejw)|'); grid on;legend('modulated signal');subplot(3,1,3);plot(w,fftshift(angle(fft(Sqq))));title('相频');xlabel('w');ylabel('∠X(ejw)');grid on;legend('modulated signal');目录第一章总论 ......................... 错误!未定义书签。