matlab利用理想采样对信号sin(t)进行取样和恢复,同时画出频谱图,所有图在一张figure中显示

时间: 2023-07-16 11:16:32 浏览: 38
好的,下面是实现代码: ```matlab % 设置采样率和采样时间 fs = 200; % 采样率 Ts = 1/fs; % 采样时间 % 设置信号参数 f = 10; % 信号频率 t = 0:1e-3:1; % 时间向量 x = sin(2*pi*f*t); % 信号 % 理想采样 n = 0:1/fs:1; % 采样时间向量 xn = sin(2*pi*f*n); % 采样信号 % 理想重构 t_interp = 0:1e-5:1; % 重构时间向量 x_interp = interp1(n,xn,t_interp,'spline'); % 重构信号 % 计算频谱 X = fft(x_interp)/length(x_interp); % 信号频谱 f_axis = linspace(-fs/2,fs/2,length(X)); % 频率轴 % 绘图 figure; subplot(3,1,1); plot(t,x,'LineWidth',2); xlabel('时间(s)'); ylabel('幅度'); title('原始信号'); grid on; subplot(3,1,2); stem(n,xn,'LineWidth',2); xlabel('时间(s)'); ylabel('幅度'); title('采样信号'); grid on; subplot(3,1,3); plot(f_axis,fftshift(abs(X)),'LineWidth',2); xlabel('频率(Hz)'); ylabel('幅度'); title('恢复信号频谱'); grid on; ``` 这段代码会生成一个包含三个子图的 figure,第一个子图为原始信号,第二个子图为采样信号,第三个子图为恢复信号频谱。

相关推荐

优化以下代码 close all; clear all; f1=40000;f2=10000;f3=20000; %信号频率 F0=1e6; %采样频率 T0=1/F0; %采样间隔 t=0:T0:10; %设置时间区间和步长 xa=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t); %原信号 %信号曲线图 figure; plot(t,xa); axis([0 0.0002 -3 3]) title('原信号'); Fs=1e5; % 抽样率大于最大频率二倍 T=1/Fs; %采样间隔 N=1000; %采样点个数 n=(0:(N-1))*T; tn=0:T:10; xn=sin(2*pi*f1*n)+sin(2*pi*f2*n)+sin(2*pi*f3*n); figure; subplot(211); stem(n,xn,'filled'); %抽样信号曲线图 axis([0 0.0002 -3 3]); title('取样信号'); subplot(212); xn_f=fft(xn); %xn_f=fftshift(fft(xn)); %傅里叶变换 f_xn=(0:length(xn_f)-1)*Fs/length(xn_f); plot(f_xn,abs(xn_f)); title('取样信号频谱'); %内插恢复原信号 t1=0:1000-T; TN=ones(length(t1),1)*n-t1'*T*ones(1,length(n)); y=xn*sinc(2*pi*Fs*TN); figure; subplot(211); plot(t1,y); axis([0 20 -3 3]); subplot(212); y_f=fft(y); %傅里叶变换 f_y=(0:length(y_f)-1)*Fs/length(y_f); plot(f_y,abs(y_f)); low_filter=hanming_low; x2=filter(low_filter,y); figure; subplot(211); plot(x2); axis([0 100 -1 1]); subplot(212); x2_f=fft(x2); %傅里叶变换 f_x2=(0:length(x2_f)-1)*Fs/length(x2_f); plot(f_x2,abs(x2_f)); title('10KHz'); high_filter=hanming_high; x1=filter(high_filter,y); figure; subplot(211); plot(x1); axis([0 100 -1 1]); subplot(212); x1_f=fft(x1); %傅里叶变换 f_x1=(0:length(x1_f)-1)*Fs/length(x1_f); plot(f_x1,abs(x1_f)); title('40KHz'); band_filter=hanming_band; x3=filter(band_filter,y); figure; subplot(211); plot(x3); axis([0 100 -1 1]); subplot(212); x3_f=fft(x3); %傅里叶变换 f_x3=(0:length(x3_f)-1)*Fs/length(x3_f); plot(f_x3,abs(x3_f)); title('20KHz');

最新推荐

recommend-type

Matlab求信号响应与频谱分析.docx

求解问题为:利用MATLAB编程,自行定义一个连续系统(2阶),求解系统的冲激响应、阶跃响应。输入信号变化时,如为f(t)=exp(-t)*u(t)时系统的输出,并画出该系统的零极点图,频率响应特性。
recommend-type

基于Java的SaaS OA协同办公毕设(源码+使用文档)

系统概述 SaaS OA协同办公系统通常包括以下几个关键组件: 用户界面(UI):提供用户交互界面,用于任务管理、日程安排、文档共享等。 后端服务:处理业务逻辑,如用户认证、数据管理、服务集成等。 数据库:存储用户数据、任务数据、文档数据等。 服务层:提供业务逻辑服务,如权限管理、工作流程等。 集成API:与其他系统集成,如邮件服务、短信服务等。 主要功能 用户认证与管理:用户登录、权限分配、用户资料管理。 任务管理:创建、分配、跟踪和归档任务。 日程管理:安排会议、提醒事件、查看日历。 文档管理:上传、下载、共享和版本控制文档。 协同工作:实时编辑文档、团队讨论、任务协作。 技术架构 Java:作为主要的编程语言。 Spring Boot:用于快速开发基于Java的后端服务。 Apache Shiro或Spring Security:用于安全和认证。 Thymeleaf或JSF:用于构建Java Web应用的用户界面。 数据库:如MySQL、PostgreSQL或MongoDB。 开发优势 实用性:解决企业日常办公需求,提高工作效率。 技术先进:使用当前流行的Java技术栈和框架。
recommend-type

虎年春节送祝福微信小程序源码下载/新版UI/支持多种流量主

虎年春节送祝福微信小程序源码下载,新版UI支持多种流量主,这是一款网友用以前发过的一款端午送祝福改的一款小程序。 里面的背景图包括祝福语都已经修改成与虎年相关的内容了,总体来说找的背景图还是可以的,不过有些地方和细节小编也给完善了一下。 然后小编测试的时候发现还没有流量主,所以小编也给加了几个流量主进去,到时候大家直接替换流量主的ID就可以了。 另外支持更多小程序推荐,拥有独立的推荐界面 PS:进入送祝福的按钮,部分机型是在老虎的帽子那里,部分是在金元宝那里
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

MATLAB正态分布协方差分析:揭示正态分布变量之间的协方差

![MATLAB正态分布协方差分析:揭示正态分布变量之间的协方差](https://site.cdn.mengte.online/official/2021/11/20211128213137293.png) # 1. 正态分布概述 正态分布,又称高斯分布,是统计学中最重要的连续概率分布之一。它广泛应用于自然科学、社会科学和工程领域。 正态分布的概率密度函数为: ``` f(x) = (1 / (σ√(2π))) * exp(-(x - μ)² / (2σ²)) ``` 其中: - μ:正态分布的均值 - σ:正态分布的标准差 - π:圆周率 正态分布具有以下特性: - 对称性:
recommend-type

我正在开发一款个人碳足迹计算app,如何撰写其需求分析文档,请给我一个范例

为了更全面、清晰地定义个人碳足迹计算app的需求,需求分析文档应该包含以下内容: 1.项目简介:对该app项目的概述及目标进行说明。 2.用户分析:包括目标用户群、用户需求、行为等。 3.功能需求:对app的基本功能进行定义,如用户登录、数据录入、数据统计等。 4.非功能需求:对使用app的性能和质量等进行定义,如界面设计、数据安全、可扩展性等。 5.运行环境:包括app的开发环境和使用环境。 下面是一个范例: 需求分析文档 1. 项目简介 该app项目旨在为用户提供一款方便、易用、可定制的个人碳足迹计算平台,以促进环保和可持续性发展。 2. 用户分析 目标用户群:全球关
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依
recommend-type

MATLAB正态分布相关性分析:探索正态分布变量之间的关联

![MATLAB正态分布相关性分析:探索正态分布变量之间的关联](https://img-blog.csdnimg.cn/bd5a45b8a6e94357b7af2409fa3131ab.png) # 1. MATLAB中正态分布的理论基础 正态分布,又称高斯分布,是一种常见的概率分布,其概率密度函数为: ``` f(x) = (1 / (σ√(2π))) * e^(-(x-μ)² / (2σ²)) ``` 其中,μ表示正态分布的均值,σ表示标准差。正态分布具有以下特点: - **对称性:**正态分布的概率密度函数关于均值μ对称。 - **钟形曲线:**正态分布的概率密度函数呈钟形曲线