INSAR实验报告

实验一、SAR信号统计模型验证实验目的通过MATLAB编程，拟合有关SAR数据的实部、虚部、灰度、相位的概率密度曲线，并与实际曲线进行比较，掌握MATLAB编程的初步知识。实验数据txt文本格式的实部和虚部数据实验步骤一、打开MATLAB，根据老师给出的实部的MATLAB的程序和拟合曲线，试着编出虚部的程序，调试并运行。二、在实部实验的基础上，接着进行以后的改编，根据各个曲线的表达式，拟合各自的曲线。实验结果实部：虚部：灰度：相位：功率：实验心得这是首次用MATLAB编程实现SAR的统计特征图像，由于之前没有用MATLAB编程实现过，所以这次实验是在学习老师所给程序的基础上进行改编的。学会了MATLAB中画图函数的编写，函数曲线的颜色、线型、多幅图像同时呈现编程方法。但是功率power的图像画的并不理想，拟合的曲线并没有很好地显示出来。通过这次的练习，平时老师课件里的东西，我们也能编程写出来，确实让我们学会了理论联系实践，学以致用。这次的实验还比较顺利，希望在以后的学习和实验中可以更多的学习一些软件，多用编程的思想实现一些自己的想法。实验源代码实部：clear;clc;load'Real.txt';load'Imag.txt';[m,n]=size(Real);Realvector=reshape(Real,m*n,1);maxreal=max(Realvector);minreal=min(Realvector);x=minreal:10:maxreal;Realvar=var(Realvector);PDF_Re=exp(-x.^2/(2*Realvar))/sqrt(2*pi*Realvar);pixel_count1=PDF_Re*10*m*n;figure;plot(x,pixel_count1,'r');count2=hist(Realvector,x);holdon;plot(x,count2,'rs','MarkerEdgeColor','k',...'MarkerFaceColor','g',...'MarkerSize',2);ylabel('Pixelcount');xlabel('Realisticcomponent(Re)');虚部：clear;clc;load'Real.txt';load'Imag.txt';[m,n]=size(Imag);Imagvector=reshape(Imag,m*n,1);maximag=max(Imagvector);minimag=min(Imagvector);x=minimag:10:maximag;Imagvar=var(Imagvector);PDF_Im=exp(-x.^2/(2*Imagvar))/sqrt(2*pi*Imagvar);pixel_count2=PDF_Im*10*m*n;figure;plot(x,pixel_count2,'r');count=hist(Imagvector,x);holdon;plot(x,count,'rs','MarkerEdgeColor','k',...'MarkerFaceColor','g',...'MarkerSize',2);ylabel('Pixelcount');xlabel('Imaginarycomponent(Im)');灰度：load'Real.txt';load'Imag.txt';[m,n]=size(Imag);Amplitude=sqrt(Real.^2+Imag.^2);Amplitudevector=reshape(Amplitude,m*n,1);Ampvar=var(Amplitudevector)*2/(4-pi);maxAmp=max(Amplitudevector);minAmp=min(Amplitudevector);X=minAmp:10:maxAmp;PDF_A=X.*exp(-X.^2/(2*Ampvar))/Ampvar;pixel_count2=PDF_A*10*m*n;figure;plot(X,pixel_count2,'r');count=hist(Amplitudevector,X);holdon;plot(X,count,'rs','MarkerEdgeColor','k',...'MarkerFaceColor','g',...'MarkerSize',2);ylabel('Pixelcount');xlabel('Amplitudecomponent(A)');功率：clear;clc;load'Real.txt';load'Imag.txt';Powe=Real.^2+Imag.^2;[m,n]=size(Powe);Powevector=reshape(Powe,m*n,1);maxPowe=max(Powevector);minPowe=min(Powevector);x=minPowe:10000:maxPowe;Powevar=var(Powevector);PDF_P=(exp(-x/(2*Powevar)))/(2*Powevar);pixel_count4=PDF_P*10*m*n;figure;plot(x,pixel_count4,'y');count=hist(Powevector,x);holdon;plot(x,count,'rs','MarkerEdgeColor','k',...'MarkerFaceColor','g',...'MarkerSize',2);ylabel('Pixelcount');xlabel('Power(P)');相位：clear;clc;load'Real.txt';load'Imag.txt';[m,n]=size(Imag);Imagvector=reshape(Imag,m*n,1);Realvector=reshape(Real,m*n,1);Phase=atan2(Imagvector,Realvector);x=-pi:pi/100:pi;PDF_w=1/(2*pi);pixel_count2=PDF_w*pi/100*m*n;figure;plot(x,pixel_count2,'r');count=hist(Phase,x);holdon;plot(x,count,'rs','MarkerEdgeColor','k',...'MarkerFaceColor','g',...'MarkerSize',2);ylabel('Pixelcount');xlabel('Phasecomponent(Ph)');实验二INSAR信号统计模型验证实验目的通过MATLAB编程实现INSAR的数据统计特征的表达，掌握用MATLAB编写函数的方法，进一步掌握MATLAB在INSAR实验中的应用。实验环境及数据MATLAB7.0实验心得通过本次实验，我学会了使用MATLAB在同一窗口中呈现多个图像的方法。这次实验中要考虑L的取值不同给函数图像带来的差异。通过这次实验，我对MATLAB更熟悉了一些。用MATLAB编程主要要理清思路，掌握每个函数的用法，这样编程才会顺利。MATLAB这个软件的功能很强大，以后要好好学习，掌握基本的用法，为以后的学习打下基础。实验结果L=1L=10L=20相干图实验源代码L1：clc;clear;L=1;Gam=0.1;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.3;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.5;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.7;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.95;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);xlabel('Phai(rad)');ylabel('PDF(Phai)');L=10clear;L=10;Gam=0.1;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.3;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.5;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.7;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.95;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);xlabel('Phai(rad)');ylabel('PDF(Phai)');L=20clear;L=20;Gam=0.1;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.3;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.5;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.7;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);holdon;Gam=0.95;Phai=[-pi:(2*pi)/100:pi];PDF=pdf_phai(Phai,Gam,L);plot(Phai,PDF);xlabel('Phai(rad)');ylabel('PDF(Phai)');调用函数的源代码functionPDF=pdf_phai(Phai,Gam,L)Beita=Gam*cos(Phai);temp1=(1-Gam.^2).^L./2./pi;temp2=gamma(2*L-1)./(gamma(L)).^2./2.^(2*L-2);temp3=(2*L-1).*Beita.*(pi./2+asin(Beita))./(1-Beita.^2).^(L+0.5);temp4=1./(1-Beita.^2).^L;ifL==1PDF=temp1.*(temp2.*(temp3+temp4));elsetemp5=0;fori=0:(L-2)temp5=temp5+gamma(L-0.5).*gamma(L-1-i).*(1+(2*i+1).*Beita.^2)./gamma(L-0.5-i)./gamma(L-1)./(1-Beita.^2).^(i+2);endtemp5=temp5./2/(L-1);PDF=temp1.*(temp2.*(temp3+temp4)+temp5);endend相干图的源代码L=1;Gam=[0:0.01:1];forPhai=-pi:(2*pi)/100:piSD=sqrt(Phai^2*(pdf_phai(Phai,Gam,L)))endplot(Gam,SD,'y