clear;
fid=fopen('C:\','r','ieee-be') %文件路徑名字 讀寫方式 字節序
for i=1:512
for j=1:512
I(i,j)=fread(fid,1,'int16'); %實數
Q(i,j)=fread(fid,1,'int16'); %虛數
A(i,j)=(Q(i,j)^2+I(i,j)^2)^0.5; %計算振幅
end
fclose(fid);
figure;imshow(A,(0,1000)); %顯示振幅強度圖像