声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1616|回复: 1

[综合] 请教Peakpicking 程序问题

[复制链接]
发表于 2010-6-2 21:13 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
这是我写的一个峰值拾取法的模态识别程序,算出来的参数和理论相差很大,而且得到的振型不对,不知道哪个地方编错了,请高手指点。或者哪位已经编好了可否提供参考。谢谢!!

clear;clc
load('strainsig.mat');                                               %信号是15通道的长度为1024的动态应变信号
fs=2000;                                                                  %采样频率
n=length(strainsig);                                                  %信号的长度
r=fft(strainsig);r=fftshift(r);                                       %对strainsig进行傅立叶变换
amp=abs(r(513:1024,:))*2/n;                                    %幅值谱
f=linspace(0,fs/2,n/2);                                               %生成一个n点的向量(从-fs/2到fs)
t=(0:1023)/fs;
figure(1)
plot(f,amp);      
set(gca,'FontSize',13)
xlabel('Frequence(Hz)');ylabel('Amplitude');
title('Spectrogram')
for k=1:15
        P(k,:)=cpsd(strainsig(:,k),strainsig(:,k),1024);               %求自功率谱密度
end
figure(2)
PSD=log(P);
ff=linspace(0,fs/2,n/2+1);
plot(ff,PSD);
title('PSD')
xlabel('frequence(Hz)');ylabel('PSD(dB)')
fig=0;
f=1;
rr=r(513:1024,:);
while f~=0
            f=input('input a frequency for modal shape or 0 to stop : ');
                if f==0,
                    break
                end
                if fig==3,
                   close(3)
                end
            Ph=180*(abs(angle(rr(round(f),:)))-pi/2)/pi                            %计算相位角
            disp(Ph');
            for i=1:15
               H(i)=(P(i,f)/P(5,f))*sign(Ph(i));                                    %选取5节点作为参考点求振型
            end
            figure(3)
            fig=3;
            hold on
            whitebg(3,'w')                                                               %Change axes background color
            set(3,'position',[200,200,600,400])
            plot([0:16],[0,H,0],'k','linewidth',2)
            plot([0,16],[0,0],'k')
            axis([0,16,-1.2,1.2])                                               %轴线显示大小
end
回复
分享到:

使用道具 举报

发表于 2010-6-28 13:39 | 显示全部楼层
感觉没原始数据(strainsig.mat), 别人不好测试!?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-29 16:14 , Processed in 0.057967 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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