xdata=fread(fid,N,'int16');这句话有问题吗

swyzmde

fs=10000;
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');

happy

没看出什么问题来

[ 本帖最后由 happy 于 2006-10-24 11:12 编辑 ]

swyzmde

那为什么执行不了啊

swyzmde

哪位专家在MATLAB上试一下,帮帮小第吧,万分感谢

eight

原帖由 swyzmde 于 2006-10-24 17:03 发表
哪位专家在MATLAB上试一下,帮帮小第吧,万分感谢

1. 无数据，如何帮你运行？
2. 无错误代码，如何得知原因？

swyzmde

我是新手,什么是无数据呀,应该怎么办,请指点,虚心学习谢谢

jimin

原帖由 swyzmde 于 2006-10-24 21:56 发表
我是新手,什么是无数据呀,应该怎么办,请指点,虚心学习谢谢

fid=fopen('bearingout.dat','r');
这里面不是有个'bearingout.dat'数据呀，没这个数据别人运行不了程序，又不想自己编数据
怎么这程序不是你自己写的吗

swyzmde

谢谢啊

swyzmde

在问一个弱弱的问题,bearingout.dat数据怎么输入,着程序不是我写的,我只分析图象啊,谢谢各位多包涵

seagull.neu

把bearingout.dat换成你自己的数据文件就能运行