马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
请问:我想读取数据全长求瞬时参量。下面程序中的前4句读取全长好像不对,如何改啊,读取的数据为*.dat
请指教
[FILENAME, PATHNAME] = uigetfile('*.dat', 'open new format VLR data file');
%打开文件并按double类型读取全长
fid = fopen([PATHNAME FILENAME],'r');
while ~(feof(fid))
[xr,n] = fread(fid,'double');
%xr = [1 4 9 3];
x = hilbert(xr');
xi=imag(x);
xr=xr';%转置成和xi一致的行列
a=sqrt(xr.^2+xi.^2);%瞬时幅度求法。
hold on;
subplot(4,1,1);
plot(xr);
title('原数据xr');
subplot(4,1,2);
plot(xi);
title('虚部xi');
subplot(4,1,3);
plot(x);
title('希尔伯特变换x');
subplot(4,1,4);
plot(a);
title('瞬时幅度a'); |