24小时热门版块排行榜    

查看: 691  |  回复: 4
【悬赏金币】回答本帖问题,作者627218299将赠送您 5 个金币

627218299

新虫 (初入文坛)

[求助] 一段书上的故障诊断程序,xdada处报错求改(金币总额太少,抱歉) 已有1人参与

%采样频率
fs=10000;
nfft=10240;
%轴承外环信号
fid=fopen('D:12K\OR007_0.mat');
N=2048;
xdata=fread(fid,N,'int16');
fclose(fid);
xdata=(xdata-mean(xdata))/std(xdata,1);
%功率谱
figure(1);
Y=abs(fft(xdata,nfft));
plot((0:nfft/2-1)/nfft*fs,Y(1:nfft/2));
xlabel('频率 f/Hz');
ylabel('功率率 P/W');
%5层小波包分解
T=wpdec(xdata,5,'db2');
%重构低频信号
y1=wprcoef(T,[5,2]);
%y1的波形
figure(2);
supplot(2,2,1);
plot(1:N,y1);
xlabel('时间,t/s');
ylabel('电压,v/V');
%y1功率谱
Y1=abs(fft(y1,nfft));
subplot(2,2,2);
plot((0:nfft/2-1)/nfft*fs,Y1(1:nfft/2));
xlabel('频率,f/Hz');
ylabel('功率率,P/W');
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

soungch

银虫 (正式写手)

【答案】应助回帖

'D:12K\OR007_0.mat'
2楼2016-12-03 15:45:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

soungch

银虫 (正式写手)

【答案】应助回帖

'D:12K\OR007_0.mat'
你有这个位置的数据吗?你报错的地方,应该是没找到fid,也就是没这个数据。
3楼2016-12-03 15:46:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王家小胖

银虫 (小有名气)

人皇
4楼2017-03-04 21:36:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wkgo12

新虫 (初入文坛)

5楼2017-04-07 17:05:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 627218299 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
信息提示
请填处理意见