陷波器程序如下,麻烦帮忙看看

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消除谐波,如何在含有整次谐波和非整次谐波的信号中去除整次谐波?...相关推荐

  1. 如何用matlab消除谐波,哈工大 matlab电机学仿真实验消除谐波

    一.实验内容及目的 1.实验内容 一个单相合成电动势中含有明显的5次谐波和7次谐波,其表达式如下所示.采用短距方法,验证削弱谐波电动势的效果. E=sin(Vs)+0.2sin(Vs)*sin(5Vs ...

  2. 为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换

    [纯技术帖]为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换来源: 陈诚--WECN的日志 写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,内容 ...

  3. 为什么要进行傅立叶变换?如何用Matlab实现快速傅立叶变换?

    http://blog.sina.com.cn/s/blog_9015f32301014jbp.html 写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,绝大部分内容非我所原创.在 ...

  4. 为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换?

    https://www.douban.com/note/164400821/ 写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,绝大部分内容非我所原创.在此向多位原创作者致敬!!! ...

  5. 【转】为什么要进行傅立叶变换?傅立叶变换究竟有何意义?如何用Matlab实现快速傅立叶变换?...

    写在最前面:本文是我阅读了多篇相关文章后对它们进行分析重组整合而得,绝大部分内容非我所原创.在此向多位原创作者致敬!!!一.傅立叶变换的由来关于傅立叶变换,无论是书本还是在网上可以很容易找到关于傅立叶 ...

  6. 用matlab处理表格,matlab删除excel表格数据-如何用matlab处理多个excel表格中的数据...

    如何用matlab处理多个excel表格中的数据 biao='D:\Program Files\matlab\bin\filename.xls'; A1=xlsread (biao,'Sheet1', ...

  7. 如何用matlab解异或方程,Matlab-6:解非线性方程组newton迭代法

    函数文件: function x=newton_Iterative_method(f,n,Initial) x0=Initial; tol=1e-11; x1=x0-Jacobian(f,n,x0)\ ...

  8. 傅立叶变换的原理、意义以及如何用Matlab实现快速傅立叶变换

    本帖最后由 xiaoliu 于 2011-7-28 21:00 编辑 一.傅立叶变换的由来 关于傅立叶变换,无论是书本还是在网上可以很容易找到关于傅立叶变换的描述,但是大都是些故弄玄虚的文章,太过抽象 ...

  9. matlab基带信号速率,如何用MATLAB产生25Gbaud的基带QPSK调制信号

    如何用MATLAB产生25Gbaud的基带QPSK调制信号 包括生成随机符号序列,进行QPSK调制,过采样,脉冲成型滤波和降采样.有错误敬请指正. 参数设置:symbol rate=25Gbaud D ...

最新文章

  1. 7-7 宇宙无敌大招呼 (C语言)
  2. windows 7系统搭建PHP网站环境
  3. linux下使用expect+scp+shell实现分布式集群系统安装,升级,部署
  4. [python调试笔记] 编译运行h5py
  5. SQL注入之时间盲注 和 报错注入(sql-lab第一关为例)
  6. 计算机课实验三,成都信息工程学院计算机网络课程实验三
  7. 【51单片机快速入门指南】1:基础知识和工程创建
  8. Sentence-BERT: 一种能快速计算句子相似度的孪生网络
  9. Linux关闭占用端口的进程
  10. MySql可视化工具MySQL Workbench使用笔记
  11. 小米3文件与电脑连接到服务器,小米3如何连接电脑_小米3连接电脑发送文件的步骤...
  12. 滴滴的2019:巨亏和裁员之后,群狼将至
  13. lighttpd http响应报文(Response)增加安全头Referrer-Policy和X-Permitted-Cross-Domain-Policies方法
  14. waterfall.js
  15. ID精度损失处理-Long类型转换成JSON的String类型
  16. 【每日一句】名人名言学英语(11.26)
  17. (附源码)计算机毕业设计ssm个人人际关系管理软件
  18. 深度学习部署:Windows安装pycocotools报错解决方法
  19. 【SAP消息号AW003】
  20. 软件设计师_计算机系统——进制加减法

热门文章

  1. 10分钟了解一致性hash算法
  2. ASP.NET Core on K8S深入学习(2)部署过程解析与部署Dashboard
  3. CAP 2.4版本发布,支持版本隔离特性
  4. Apache SkyWalking的架构设计【译文】
  5. C#使用Xamarin开发可移植移动应用进阶篇(8.打包生成安卓APK并精简大小),附源码
  6. ArcGIS实验教程——实验十七:缓冲区分析(Buffer Analysis)
  7. ArcGIS实验教程——实验三:矢量数据采集与编辑(矢量化)
  8. C和指针之数组编程练习5 (矩阵相乘)
  9. Tree前序反序列化
  10. Android之支付宝设计与开发