How to have the data for each iteration for a for loop?

1 次查看（过去 30 天）
Tharindu
Tharindu 2014-11-19
评论： Tharindu 2014-11-21
Hi I have two sets of acceleration excel sheets data (8192x50 each) and I need to perform Fast Fourier Transformation for the total set of columns taking the same column from each file. As I have 50 columns I have to do this for 50 times. I used the for loop but it will only give the results for the last iteration. please help me to get the whole set of data.
here is the script I used for the calculation;
a=xlsread('V02_All_column5.xlsx');
b=xlsread('V02_All_column6.xlsx')
for i =1:50
z=a(:,i);
x=b(:,i);
Fs=3200;
L=length(z);
NFFT=2^nextpow2(L);
window=hann(L/2);
noverlap=8;
[Pff(:,:),f]=pwelch(z(:,:),window,noverlap,NFFT,Fs);
[Pxx(:,:),f]=pwelch(x(:,:),window,noverlap,NFFT,Fs);
[Pfx(:,:),f]=cpsd(z(:,:),x(:,:),window,noverlap,NFFT,Fs);
[Pxf(:,:),f]=cpsd(x(:,:),z(:,:),window,noverlap,NFFT,Fs);
H1(:,:)=Pfx(:,:)./Pff(:,:);
H2(:,:)=Pxx(:,:)./Pxf(:,:);
end
I want to get the H1 values for the whole iteration. Please help me.
Thank You .........!!!

请先登录，再回答此问题。

采纳的回答

per isakson
per isakson 2014-11-19
编辑：per isakson 2014-11-20
Replace
H1(:,:)
by
H1(:,:,i)
and initialize H1 before the loop.
&nbsp
Addendum
"initialize" is short for Preallocating Memory
Add this line before the loop
H1 = zeros( n1, n2, 50 );
choose n1 and n2 so that H1 gets its final size
  3 个评论
显示 1更早的评论
Tharindu
Tharindu 2014-11-21
Hi many thanks for your kind assistance.
I think it is working now. But don't know how to see the data. Because in the command window it only shows some values of the H1. When I try to export it MS excel using xlswrite('1.xlswrite',real(H1));
it returns the error message 'Dimension of input array cannot be higher than two.'
Can you explain me how to write it to an excel sheet. Please note that for each iteration H1 should have 4097 measurements. I would love if I can have the results in a 4097x50 matrix rather than having all the iteration results in the same column.
Many thanks again.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by