发布网友 发布时间:2024-01-10 15:34
共1个回答
热心网友 时间:2024-01-11 00:51
clearall%清楚工作空间的所有变量a=load('imaging2.txt');%加载txt文档N1=a(1)/a(2)+1;N2=a(3)/a(4)+1;N3=a(5)/a(6)+1;%根据文档里面的数计算N1~N3t=N3;b=[];%空矩阵ft(512)=0;%定义一个1行512列的0矩阵forn1=1:N1%从1到N1循环forn2=1:N2%第二层从1到N2循环m1=7+(n1-1)*N2*t+(n2-1)*t;%计算m1m2=(n1-1)*N2*t+n2*t+6;%计算m2ft(1,1:N3)=a(m1:m2);%把a的第m1到m2个元素赋值给ft的第1到第N3个变量;这里估计会出错,不能保证维度一致aa(1,:)=fft(ft);%把fft的ft个下标赋值给aa第1行所有列b(n1,n2,:)=abs(aa(1,:));%计算绝对值endend%结束循环surf(b(:,:,40));%画三维图view(0,90)%调整视角gridoff%不显示网格线