这是我的程序

function [wa]=Wa0(x0,y0)

syms t tao

[m,~]=size(x0);  %求xx的维数

xx=0:0.001:13;  %定义绘图坐标轴

t=0:0.001:13;

plot(x0,y0,'*');

hold on;

fs=2000;N=20000;   %采样频率和数据点数

n=0:N-1;

t=n/fs;   %时间序列

%x=y0; %信号

W1=detrend(y0); % 消除时间序列中的线性趋势项

W1=detrend(W1,'constant'); % 消除时间序列中的均值

subplot(2,2,1),plot(x0,W1,'*');

yy=fftshift(fft(W1,N));    %对信号进行快速Fourier变换

mag=abs(yy);     %求得Fourier变换后的振幅

mag=mag/(N/2);   %换算成实际的幅度

mag(1)=mag(1)/2;

f=(-N/2:N/2-1)*fs/N;%进行对应的频率转换

subplot(2,2,2),plot(f,mag);   %绘出随频率变化的振幅

xlabel('频率/Hz');

ylabel('振幅');grid on;

subplot(2,2,3),plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅

xlabel('频率/Hz');

ylabel('振幅');

title('振幅-频率曲线图');grid on;

Pyy=[1:N/2];

for i=1:N/2

Pyy(i)=phase(yy(i)); %计算相位

Pyy(i)=Pyy(i)*180/pi; %换算为角度

end;

subplot(2,2,4),plot(f(1:N/2),Pyy(1:N/2));   %显示相位图

title('相位-频率曲线图');

xlabel('相位/度');

ylabel('振幅');grid on;

1.jpg

(169.4 KB, 下载次数: 6)

2015-10-14 11:13 上传

2015-10-14 11:12 上传

点击文件名下载附件

2.81 MB, 下载次数: 40

matlab求出过冲的最大幅值,频谱分析为什么最大幅值处频率为0?相关推荐

  1. 已知 方程 用 matlab 求表达式,已知自变量,因变量和函数表达式,可以用matlab求出函数表达式中的未知参数吗...

    答:这不是解方程,这是曲线拟合.如果你的已知的dr,R只有一个点,那么你的未知数是有无数个解的. cftool,拟合函数有一个选项是自定义函数. 答:用polyfit函数,即多项式拟合函数 设x=X- ...

  2. Java - 用数组求出斐波那契数列的前20项值

    斐波那契数列 斐波那契数列(Fibonacci sequence),又称黄金分割数列.因数学家列昂纳多·斐波那契(Leonardoda Fibonacci)以兔子繁殖为例子而引入,故又称为" ...

  3. matlab矩阵保存到表格,将matlab求出的矩阵保存在Excel表格中

    将matlab求出的矩阵保存在Excel表格中 转贴Genial分享程序 function xlswrite(m,header,colnames,filename); % xlswrite Easil ...

  4. 如何用matlab求出覆冰图像像素点,输电线路覆冰图像的Canny算子改进算法研究

    输电线路覆冰图像的Canny算子改进算法研究 摘 要: 每到冬季 关键词: 形态学滤波; Canny算子; 边缘检测; 覆冰检测; 微分算子 中图分类号:TM755 文献标志码:A文章编号:1006- ...

  5. c编程:求出4#215;4矩阵中最大和最小元素值及其所在行下标和列下标,求出两条主对角线元素之和。...

    //求出4×4矩阵中最大和最小元素值及其所在行下标和列下标,求出两条主对角线元素之和 #include <stdio.h> int main() {int sum=0;int max,mi ...

  6. matlab求曲线每点的斜率_如何用matlab求出图中各条直线的斜率

    展开全部 1.瞬时斜率62616964757a686964616fe4b893e5b19e31333431373236:求导. 方法:d(k)=(u(k)-u(k-1))/T(k);d(0)=0; d ...

  7. matlab求桌子频率,怎样用matlab求出桌子的边缘曲线和开槽长度

    匿名用户 1级 2014-09-14 回答 L=120;  %长 D=50;   %宽,圆桌面直径 d=2.5;  %木板宽 hL=L/2; %半长 R=D/2;  %圆桌面半径 y=-R+d/2:d ...

  8. Matlab求出图像中每个细胞的细胞核与核仁面积之比代码

    clc;clear;close all; img = imread('1.bmp'); I = rgb2gray(img); figure(1),imshow(I) level =graythresh ...

  9. matlab怎么求两个数的和,matlab怎么求出两个函数的交点

    matlab中,两个自变量的函数怎么求最大值(急!1) 需求:利用matlab求解二元函数y=f(x1,x2)=(339-0.01*x1-0.003*x2)*x1+(399-0.004*x1-0.01 ...

最新文章

  1. Redis 在CentOS 6上的 安装和部署以及redis的主从复制sentinel实现HA
  2. 我的世界java版怎么打不开_JAVA版我的世界打不开,求助!
  3. jquery自动触发事件
  4. 计算机创造奇迹英语作文,高中英语作文模板 第274期:Creating Miracle 创造奇迹
  5. 登录iis7 网站服务器上网站网页可以访问 但远程访问网站网页不能访问的解决方法
  6. Java作业基础知识梳理--2020-12-03
  7. STL之accumulate
  8. 科学家之殁,竟是区块链之过?
  9. 服务器端 viewstate
  10. 2.nginx 配置
  11. 利用ajaxSubmit()提交表单后,子窗口关闭,父窗口显示查询页面
  12. 数据结构 队列的结构特点及基本操作
  13. Redundant Variable 多余变量
  14. Java开发报表——Grid++Report 报表设计器
  15. InputStreamReader乱码
  16. (以三星S8为例)安卓全面屏手势设置教程
  17. Linux磁盘术语描述
  18. 在Linux下驱动D-link DFE-530TX(最终稿)(转)
  19. 2022全网最全的爆款抓包工具,各有千秋
  20. [zz]64bitnbsp;ubuntunbsp;安装32bit软件

热门文章

  1. 担心写出的汽车软文没人看?
  2. 2022年全球市场落石防护网总体规模、主要生产商、主要地区、产品和应用细分研究报告
  3. 对公司绩效考核办法的一些个人意见
  4. 如何推动电信 计算机和信息服务贸易,信息通信技术推动国际贸易变革 我国数字贸易规模扩大...
  5. 基于维基百科的用户意图分类
  6. 一网打尽、详解面试最常被问到的几种 Python 设计模式
  7. 基于FRR全面解析BGP协议(八):FRR配置实例
  8. 玩转华为ENSP模拟器系列 | 配置SR-MPLS BE隧道的TI-LFA FRR功能示例
  9. 计算器算贝塞尔公式_怪物猎人:近战技能收益计算器Ver.4打死不改版
  10. 独家丨专访雅捷信息董事长、NVIDIA全球副总裁,看“非主流”的GPU数据库如何升级银行数据查询与加工...