设为首页收藏本站
切换到宽版

电力仿真论坛

 找回密码
 立即注册

QQ登录

只需一步，快速开始

搜索
热搜: pscad atp VIP会员
电力仿真论坛»论坛 › 软件区 › MATLAB › MATLAB中prony分析法的问题

小喇叭+ 发布

03-01 22:05
admin：等待验证会员：1）邮箱注册需要激活邮件，找不到邮件请查找垃圾箱；2）先注册论坛账号再绑定QQ
04-04 15:39
03-12 12:40
新来朋友

dyy_zhangwm

han

Tupolev144

林柒

打工也要做研究

freedom

wzh

Sisyphus

wu_haifei

pop撒

cyj-ATP

继电保护

fy030509

yoonae

fengxian

luzhiyuan

节度使

szg0933

Iris

YANG

柔直小白

messiliu10

ywwy

18720669717

qq954646921
返回列表 发新帖
查看: 3.7K|回复: 1

MATLAB中prony分析法的问题

[复制链接]
发表于 2020-10-27 17:19:10 | 显示全部楼层 |阅读模式

注册可看大图、可下载

您需要 登录 才可以下载或查看，没有账号？立即注册

x
现在我根据论文搞了一个prony的程序，但是只有输出的四个量，我想把导入的数据画成横轴为频率，纵轴为幅值的频谱图，应该怎么改

  1. x=0:0.0001:1;   %原信号
  2. y=V22(3:end);   %提取
  3. % h=fopen('ws.txt','w');

  5. p=6     %p为模型阶数
  6. sf=1/0.0001    %sf为采样频率
  7. [F,D,A,theta]=exprony_ma(y,p,sf)
  8. %调用prony算法函数
  9. function [F,D,A,theta]=exprony_ma(y,p,sf)
  10. nm=2*p;                          %%% p为模型的阶数
  11. n=fix(length(y)/2);
  12. dt=1/sf;                        %%% 由采样频率计算时间间隔
  13. h=y(1:2*n);
  14. L=length(h);
  15. M=L/2;
  16. for k=1:nm
  17.     x1(:,k)=h(k:M-1+k);
  18. end
  19. for k=1:M
  20.     x2(k,:)=-h(nm+k);
  21. end
  22. B=x1\x2;           %% 最小二乘法求解prony 多项式系数
  23. B(nm+1)=1;
  24. B1=B(nm+1:-1:1);
  25. V=roots(B1);        %% 求特征方程的根
  26. F1=abs(log(V))/(2*pi*dt);  %% 初步计算频率
  27. D1=log(abs(V))/(2*pi*dt);  %% 初步计算衰减因子
  28. for k=0:(2*n-1)
  29.     Va(k+1,:)=conj(V').^k;  %% 范得蒙矩阵 Va
  30. end
  31. S1=2*(inv(conj(Va')*Va)*conj(Va')*h);
  32. [F2,I]=sort(F1);        %%%频率F1从小到达排序
  33. m=0;
  34. for k=1:nm-1
  35.     if F2(k)~=F2(k+1)   %%去掉其中的共轭项
  36.         continue;
  37.     end
  38.     m=m+1;
  39.     ii=I(k);           %% 提取有效的元素序数
  40.     F(m)=F1(ii);
  41.     D(m)=D1(ii);
  42.     A(m)=abs(S1(ii));
  43.     theta(m)=angle(S1(ii))*180/pi;
  44. end
  45. end
复制代码

谁来看过此贴
此贴被TA们浏览了3656次
回复

使用道具 举报

发表于 2020-12-15 20:10:14 | 显示全部楼层
请问楼主，是不是Plot(x,y)，横轴是表示频率的变量，纵轴是幅值的变量，就可以了哇。。plot(A(m),F1)。我觉得是这样修改的，你也是这样绘制的嘛
回复

使用道具 举报

返回列表 发新帖
发贴规则：　
1.严禁将帖子发至无关版面，请选择对应版块发贴，以维护论坛的系统性和整洁性。
2.提问题需要将问题描述清楚，涉及到仿真模型问题需要添加报错图片或描述，上传仿真模型效果更佳。
3.由于论坛的时效性不足，可以发帖后点击楼层下分享到：QQ好友和群快速得到关注。
4.保持和谐。
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Archiver|手机版|小黑屋|电力仿真论坛

GMT+8, 2025-5-3 11:21 , Processed in 0.369068 second(s), 40 queries .

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表