声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1531|回复: 3

[编程技巧] 书中自带程序中的'bearingout.dat文件不存在,程序运行不了,此问题该如何解决?

[复制链接]
发表于 2008-12-5 12:35 | 显示全部楼层 |阅读模式

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

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

x
《MATLAB时频分析技术及其应用》自带光盘有中如下一段程序,可是其中'bearingout.dat文件不存在啊,程序运行不了,我急需此程序实验图形,那位高手有此程序的运行结果啊,或是此问题该如何解决?谢谢各位达人了。
fid=fopen('bearingout.dat','r');
N=1024;
xdata=fread(fid,N,'int16');
fclose(fid);
xdata=(xdata-mean(xdata))/std(xdata,1);
figure(1);
plot(1:N,xdata);
xlabel('时间 t/n');
ylabel('电压 V/v');
[c,l]=wavedec(xdata,4,'db10');
d4=wrcoef('d',c,1,'db10',4);
d3=wrcoef('d',c,1,'db10',3);
d2=wrcoef('d',c,1,'db10',2);
d1=wrcoef('d',c,1,'db10',1);
figure(2)
subplot(4,1,1);
plot(d4,'LineWidth',2);
ylabel('d4');
subpot(4,1,2);
plot(d3,'LineWidth',2);
ylabel('d3');
subpot(4,1,3);
plot(d2,'LineWidth',2);
ylabel('d2');
subpot(4,1,4);
plot(d1,'LineWidth',2);
ylabel('d1');
xlabel('时间 t/s');
y=hilbert(d1);
ydata=abs(y);
y=y-mean(y);
nfft=1024;
p=abs(fft(ydata,nfft));
figure(3);
plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));
xlabel('频率 f/Hz');
ylabel('功率谱 P/W');
回复
分享到:

使用道具 举报

发表于 2011-10-2 19:04 | 显示全部楼层
《MATLAB时频分析技术及其应用》自带光盘有中如下一段程序,可是其中'bearingout.dat文件不存在啊,程序运行不了,我急需此程序实验图形,那位高手有此程序的运行结果啊,或是此问题该如何解决?谢谢各位达人了。
fid=fopen('bearingout.dat','r');
N=1024;
xdata=fread(fid,N,'int16');
fclose(fid);
xdata=(xdata-mean(xdata))/std(xdata,1);
figure(1);
plot(1:N,xdata);
xlabel('时间 t/n');
ylabel('电压 V/v');
[c,l]=wavedec(xdata,4,'db10');
d4=wrcoef('d',c,1,'db10',4);
d3=wrcoef('d',c,1,'db10',3);
d2=wrcoef('d',c,1,'db10',2);
d1=wrcoef('d',c,1,'db10',1);
figure(2)
subplot(4,1,1);
plot(d4,'LineWidth',2);
ylabel('d4');
subpot(4,1,2);
plot(d3,'LineWidth',2);
ylabel('d3');
subpot(4,1,3);
plot(d2,'LineWidth',2);
ylabel('d2');
subpot(4,1,4);
plot(d1,'LineWidth',2);
ylabel('d1');
xlabel('时间 t/s');
y=hilbert(d1);
ydata=abs(y);
y=y-mean(y);
nfft=1024;
p=abs(fft(ydata,nfft));
figure(3);
plot((0:nfft/2-1)/nfft*fs,p(1:nfft/2));
xlabel('频率 f/Hz');
ylabel('功率谱 P/W');
发表于 2012-1-5 09:36 | 显示全部楼层
你有光盘?
发表于 2012-3-19 19:30 | 显示全部楼层
你好,你找到实验用的数据了吗?可以给我发一份吗?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-9-21 07:59 , Processed in 0.054924 second(s), 19 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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