matlab计算截断误差,Matlab相位截断误差仿真解说.doc
Matlab仿真
误差频谱仿真:
N=14 M=4 l=8 k=1
程序:
clear;
Fs =1;%采样频率1MHz
M=4;%截取累加器的高M位;
N = 14;%累加器的位数;
l=8;
k=1;
K=2^l*(2*k+1);
L=2^(N-M-l);
pe= 100000;
n=1:pe;
pp=pe+1;%总的采样点数
add_y=(n*K)/2^N;
error=(round(n*(2*k+1))/L)/(2^M);
NFFT = 2^nextpow2(pp);%设定频谱分析点数
f = Fs/2*linspace(0,1,NFFT/2+1);% 只显示 0?Fs/2 的频率
rom_y=sin(2*pi*error);%有相位截断误差
Y = fft(rom_y,NFFT)/pp;% NFFT 点的 FFT 变换
value_Y = abs(Y(1:NFFT/2+1));% 求 FFT 变换值的绝对值
DB_Y = 20*log10(value_Y);% 求功率的分贝
plot(f,value_Y); %绘制频谱图,只有相位截断误差
axis([0 1 0 0.5]);
N=14 M=4 l=7 k=4
程序:
clear;
Fs =1;%采样频率1MHz
M=4;%截取累加器的高M位;
N=14;%累加器的位数;
l=7;
k=4;
K=2^l*(2*k+1);%频率控制字
L=2^(N-M-l);
pe= 100000;
n=1:pe;
pp=pe+1;%总的采样点数
add_y=(n*K)/2^N;
error=(round(n*(2*k+1))/L)/(2^M);
NFFT = 2^nextpow2(pp);%设定频谱分析点数
f = Fs/2*linspace(0,1,NFFT/2+1);% 只显示 0?Fs/2 的频率
rom_y=sin(2*pi*error);%有相位截断误差
Y = fft(rom_y,NFFT)/pp;% NFFT 点的 FFT 变换
value_Y = abs(Y(1:NFFT/2+1));% 求 FFT 变换值的绝对值
DB_Y = 20*log10(value_Y);% 求功率的分贝
plot(f,value_Y); %绘制频谱图,只有相位截断误差
axis([0 1 0 0.5]);%包括坐标轴范围,axis([xmin xmax ymin ymax])
N=14 M=4 l=9 k=0
程序:
clear;
Fs =1;%采样频率1MHz
M=4;%截取累加器的高M位;
N = 14;%累加器的位数;
l=9;
k=0;%频率控制字,即累加的步长;
K=2^l*(2*k+1);
L=2^(N-M-l);
m=1:L-1
pe= 100000;
n=1:pe;
pp=pe+1;%总的采样点数
error=(2*m*pi)/2^(N-l);
NFFT = 2^nextpow2(pp);%设定频谱分析点数
f = Fs/2*linspace(0,1,NFFT/2+1);% 只显示 0?Fs/2 的频率
rom_y=sin(2*pi*error);%有相位截断误差
Y = fft(rom_y,NFFT)/pp;% NFFT 点的 FFT 变换
value_Y = abs(Y(1:NFFT/2+1));% 求 FFT 变换值的绝对值
DB_Y = 20*log10(value_Y);% 求功率的分贝
plot(f,value_Y); %绘制频谱图 没有幅度量化误差,只有相位截断误差
axis([0 1 0 0.5]);
N=14 M=5 l=10 k=0
程序:
clear;
Fs =1;%采样频率1MHz
M=5;%截取累加器的高M位;
N = 14;%累加器的位数;
l=10;
k=0;%频率控制字,即累加的步长;
K=2^l*(2*k+1);
L=2^(N-M-l);
pe= 100000;
n=1:pe;
pp=pe+1;%总的采样点数
add_y=(n*K)/2^N;
NFFT = 2^nextpow2(pp);%设定频谱分析点数
f = Fs/2*linspace(0,1,NFFT/2+1);% 只显示 0?Fs/2 的频率
rom_y=sin(2*pi*add_y);%有相位截断误差
Y = fft(rom_y,NFFT)/pp;% NFFT 点的 FFT 变换
value_Y = abs(Y(1:NFFT/2+1));% 求 FFT 变换值的绝对值
DB_Y = 20*log10(value_Y);% 求功率的分贝
plot(f,value_Y); %绘制频谱图
matlab计算截断误差,Matlab相位截断误差仿真解说.doc相关推荐
- 传输预编码matlab,基于MATLAB的MIMO系统预编码性能仿真教程.doc
基于MATLAB的MIMO系统预编码性能仿真教程 PAGE \* MERGEFORMAT - 33 - 摘要在现今的移动通信系统中,被极多的国际通信标准采纳为基础性关键技术的一种方法是多输入多输出的技 ...
- matlab 计算 r,MATLAB R2015b*化计算
目录 第1章MATLAB R2015b概述 1.1MATLAB简介 1.1.1MATLAB的发展史 1.1.2MATLAB的优点 1.1.3MATLAB系统组成 1.1.4MATLAB应用程序 1.1 ...
- matlab 计算 工程,matlab工程计算.doc
您所在位置:网站首页 > 海量文档  > 计算机 > matlab matlab工程计算.doc167页 本文档一共被 ...
- 求pi的数学模型matlab计算pi,MATLAB实验-pi的计算 圆周率的近似计算 数学软件与数学实验 教学课件.ppt...
MATLAB实验-pi的计算 圆周率的近似计算 数学软件与数学实验 教学课件 1.圆周率π的计算历程 所谓"圆周率"是指一个圆的周长与其直径的比值.古今中外,许多人致力于圆周率的研 ...
- matlab计算 var,matlab 在险价值 VaR 的计算
matlab 在险价值 VaR 的计算 matlab 在险价值 VaR 的计算 VaR 模型 数据获取 历史模拟法 蒙特卡罗模拟法 参数模型法 代码和数据下载 VaR 模型 Value at Risk ...
- 基于matlab信道均衡仿真,基于MATLAB的数字自适应信道均衡器仿真研究.doc
摘要:科技的迅速发展带动了我们生活中各个领域的进步,通信领域也不例外.在基带传输的通信系统中,实际波形不可避免的会发生失真,因为时变通信质量与传输速度的主要干扰因素是码间串扰ISI.信道均衡技术[2] ...
- matlab计算导纳矩阵,matlab实现导纳矩阵
<matlab实现导纳矩阵>由会员分享,可在线阅读,更多相关<matlab实现导纳矩阵(3页珍藏版)>请在人人文库网上搜索. 1.Matlab形成节点导纳矩阵学号: 姓名:侯成 ...
- 计算关联度矩阵matlab,计算关联度matlab
% 计算关联度 m=9; %子序列 k=11;%序列长度 X0=zeros(11,9);%差序列矩阵 R0=zeros(11,9);%关联系数矩阵 % 第一步:载入源数据 a0=[10353 1058 ...
- matlab 计算det,matlab数学运算与函数
① 基本代数运算操作 +,-,*,\,/,^, ② 矩阵运算函数:求行列式(det),矩阵求逆(inv),求秩(rank),求迹(trace), 求模(norm),d=eig(A)求矩阵A 的特征值, ...
最新文章
- 基于RFID的防伪系统设计
- 安装xml2js出现npm ERR! code E404 npm ERR! 404 Not Found: event-stream@3.3.6
- COCO数据集数据转换为XML格式
- mysql 4 中文模糊查询_解决MySQL中文模糊查询问题
- java 根据类名示例化类_Java LocalDateTime类| 带示例的getNano()方法
- 获取网关_阿里二面问了这道题:如何设计一个微服务网关系统
- 触发器中的 临时表 old 与 new
- BZOJ1283 序列(费用流)
- 腾讯AI开放平台使用尝试:通过文本翻译API进行汉译英
- SQL语句实现增删改查(1)
- N卡自带录屏软件geforce 双屏录制问题
- 收集的13个杀毒软件和安全防护软件(有图哦)
- 初来乍到:新用户冷启的算法技巧
- 瑞利分布的平方是什么分布
- java catch空指针异常处理_Java基础学习:java文件空指针异常处理
- python 多态app_多态简介 | Python从入门到精通:高阶篇之三十三-阿里云开发者社区...
- mac 无法打开“XXXXXXX.app”,因为无法确认开发者的身份。
- java IO读写文件详细介绍
- 列举常用字节输入流和输出流
- 520 情人节 七夕 表白动画
热门文章
- 太划算了! 1 元秒杀 1000 本爆款电子书!
- UVA 10242 || Fourth Point !!(求平行四边形第四点
- 格式化jsp/html文件
- MacOS Big Sur 系统 【您没有权限来打开应用程序“xf-adesk19”】
- Android 图片选择库美哭了
- Spring Boot: Bean definition overriding
- layui的html页面支持转jsp吗,layui jsp
- 五次多项式对2R平面机器人关节的轨迹控制
- 达梦数据库DCA培训笔记(基于dm8)
- 欧盟对谷歌开出第三张天价罚单 累计罚款总额达93.7亿美元