如何用matlab消除谐波,如何在含有整次谐波和非整次谐波的信号中去除整次谐波?...
陷波器程序如下,麻烦帮忙看看
clc;
clear all;
close all;
x = 0 : 0.05 : 10;
s= 210 * exp(-0.25 * x) .* cos(2 * pi * 1.5 * x + pi) + 110 * exp(0.25 * x) .* cos(2 * pi * 1.0 * x + pi/2 );
%load data;
%s=cc;
N=length(s); % 信号长度
fs=1000; % 采样频率
n=1:N;
n2=1:N/2;
tt=(n-1)/fs; % 时间刻度
ff=(n2-1)*fs/N; % 频率刻度
X=fft(s); % 谱分析
figure(1)
subplot 311;
plot(tt,s,'k');
title('原始数据'); xlabel('时间/s'); ylabel('幅值');
X=X/max(abs(X));
subplot 312;
plot(ff,abs(X(n2)),'k');
axis tight; title('谱分析');
xlabel('频率/Hz'); ylabel('幅值');
subplot 313;
plot(tt,abs(X(n)));
% for k=1:5 % 自适应陷波器
% j=(k-1)*2+1; % 设置50Hz和它的奇次谐波频率
% f0=116.73*j;
% j=(k-1)*2+1; % 设置50Hz和它的奇次谐波频率
f0=116.73;
x1=cos(2*pi*tt*f0); % 设置x1和x2
x2=sin(2*pi*tt*f0);
w1=0; % %初始化w1和w2
w2=100;
e=zeros(1,N); % %初始化e和y
y=zeros(1,N);
mu=0.0000001; % 设置迭代步长
for i=1:N % 自适应陷波器
y(i)=w1*x1(i)+w2*x2(i); % 计算y
e(i)=s(i)-y(i); % 计算e
w1=w1+mu*e(i)*x1(i); % 调整w
w2=w2+mu*e(i)*x2(i);
end
x=e;
% end
output1=e; % 陷波器输出
如何用matlab消除谐波,如何在含有整次谐波和非整次谐波的信号中去除整次谐波?...相关推荐
- 如何用matlab消除谐波,哈工大 matlab电机学仿真实验消除谐波
一.实验内容及目的 1.实验内容 一个单相合成电动势中含有明显的5次谐波和7次谐波,其表达式如下所示.采用短距方法,验证削弱谐波电动势的效果. E=sin(Vs)+0.2sin(Vs)*sin(5Vs ...
- 为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换
[纯技术帖]为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换来源: 陈诚--WECN的日志 写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,内容 ...
- 为什么要进行傅立叶变换?如何用Matlab实现快速傅立叶变换?
http://blog.sina.com.cn/s/blog_9015f32301014jbp.html 写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,绝大部分内容非我所原创.在 ...
- 为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换?
https://www.douban.com/note/164400821/ 写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,绝大部分内容非我所原创.在此向多位原创作者致敬!!! ...
- 【转】为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换?...
写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,绝大部分内容非我所原创.在此向多位原创作者致敬!!!一.傅立叶变换的由来关于傅立叶变换,无论是书本还是在网上可以很容易找到关于傅立叶 ...
- 用matlab处理表格,matlab删除excel表格数据-如何用matlab处理多个excel表格中的数据...
如何用matlab处理多个excel表格中的数据 biao='D:\Program Files\matlab\bin\filename.xls'; A1=xlsread (biao,'Sheet1', ...
- 如何用matlab解异或方程,Matlab-6:解非线性方程组newton迭代法
函数文件: function x=newton_Iterative_method(f,n,Initial) x0=Initial; tol=1e-11; x1=x0-Jacobian(f,n,x0)\ ...
- 傅立叶变换的原理、意义以及如何用Matlab实现快速傅立叶变换
本帖最后由 xiaoliu 于 2011-7-28 21:00 编辑 一.傅立叶变换的由来 关于傅立叶变换,无论是书本还是在网上可以很容易找到关于傅立叶变换的描述,但是大都是些故弄玄虚的文章,太过抽象 ...
- matlab基带信号速率,如何用MATLAB产生25Gbaud的基带QPSK调制信号
如何用MATLAB产生25Gbaud的基带QPSK调制信号 包括生成随机符号序列,进行QPSK调制,过采样,脉冲成型滤波和降采样.有错误敬请指正. 参数设置:symbol rate=25Gbaud D ...
最新文章
- 7-7 宇宙无敌大招呼 (C语言)
- windows 7系统搭建PHP网站环境
- linux下使用expect+scp+shell实现分布式集群系统安装,升级,部署
- [python调试笔记] 编译运行h5py
- SQL注入之时间盲注 和 报错注入(sql-lab第一关为例)
- 计算机课实验三,成都信息工程学院计算机网络课程实验三
- 【51单片机快速入门指南】1:基础知识和工程创建
- Sentence-BERT: 一种能快速计算句子相似度的孪生网络
- Linux关闭占用端口的进程
- MySql可视化工具MySQL Workbench使用笔记
- 小米3文件与电脑连接到服务器,小米3如何连接电脑_小米3连接电脑发送文件的步骤...
- 滴滴的2019:巨亏和裁员之后,群狼将至
- lighttpd http响应报文(Response)增加安全头Referrer-Policy和X-Permitted-Cross-Domain-Policies方法
- waterfall.js
- ID精度损失处理-Long类型转换成JSON的String类型
- 【每日一句】名人名言学英语(11.26)
- (附源码)计算机毕业设计ssm个人人际关系管理软件
- 深度学习部署:Windows安装pycocotools报错解决方法
- 【SAP消息号AW003】
- 软件设计师_计算机系统——进制加减法
热门文章
- 10分钟了解一致性hash算法
- ASP.NET Core on K8S深入学习(2)部署过程解析与部署Dashboard
- CAP 2.4版本发布,支持版本隔离特性
- Apache SkyWalking的架构设计【译文】
- C#使用Xamarin开发可移植移动应用进阶篇(8.打包生成安卓APK并精简大小),附源码
- ArcGIS实验教程——实验十七:缓冲区分析(Buffer Analysis)
- ArcGIS实验教程——实验三:矢量数据采集与编辑(矢量化)
- C和指针之数组编程练习5 (矩阵相乘)
- Tree前序反序列化
- Android之支付宝设计与开发