资源描述：
《mr原理、技术与设备（冯衍秋）模拟bloch方程的受激（excitation）过程》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、1.模拟Bloch方程的受激(Excitation)过程。答：模拟的MATLAB脚本内容如下，共1个dT=0.5;%间隔ImsT=500;%总时间N=ceil(T/dT)+1;time=[0:0.1:N-l]*dT;df=10;wl=2*pi*df/10000;w2=2*pi*df/500;z=cos(wl•*time);x=sin(wl•壮ime)•*sin(w2.壮ime);y=sin(wl•壮ime)•*cos(w2.*time);plot3(x,y,z);grid
2、onwhile1comet3(x,z);pause(0.2)end结果如下三图（动画截图）:图一：图二图三:QFigure11•模拟Bloch方程的弛豫(Relaxation)过程。答：模拟MATLAB脚本如下，共3个：dT=0.05;%间隔1msT=1000;%总吋间N=ceil(T/dT)+l;%时间间隔数df=10;%0.1s后关闭共振场T1=600;T2=100;%生成适当的种子矩阵[AB]=freeprecess(dTzTl,T2Zdf);t
3、ime=[0:0.1:N-l]*dT;%=====模拟衰减=====M=zeros(3ZN);%实时跟踪磁化矢量末端位置M(:,l)=[l;0;0];%加磁场fork=2:NM(:,k)=A*M(:zk-1)+B;end;plot3(M(l,:),M(2,:),M(3,:))gridonwhile1comet3(M(lz:),M(2,:),M(3,:)z0.01);pause(0.5);end<1—文件结束—>function[AfpzBfp]=freeprecess(
4、TzTlzT2Zdf)%%模拟口rfl进动和消共振的函数%经过T时间，给定弛豫时间Tl、T2%和消共振频率df・时间单位：msz消共振频率：Hz.phi=2*pi*df*T/1000;%弛豫进动半径El=exp(-T/T1);E2=exp(-T/T2);Afp=[E200;0E20;00El]*zrot(phi);Bfp=[001-E1]1;<1—文件结束—>functionRz=zrot(phi)%%这是z轴旋转变换函数Rz=[cos(phi)-sin(phi)0;sin(phi)
5、cos(phi)0;001]结果如下三图(动画截图)：图一:图二图三:1.给定频谱图像headraw,按下列要求进行傅里叶逆变换重建图像:A.将数据第1到64行数据置0B.将数据1到64行与193到256行置0C.选择数据奇数行D.选择数据偶数行E.将数据放到512*512大小的矩阵中心，其他位置处为0答：各条件下代码分列如下，附结果图像：A.fori=l:64forj=l:256newheadraw(izj)=0;endfori=65:256forj=l:256newh
6、eadraw(i,j)=headraw(i,j);endendI=ifft2(newheadraw);Is=fftshift(I);Ir=real(Is);figure,imshow(Ir)fori=l:64forj=l:256newheadraw(izj)=0;endfori=65:192forj=l:256newheadraw(irj)=headraw(iAj);endendfori=193:256forj=l:256newheadraw(i,j)=0;en
7、dendI=ifft2(newheadraw);Is=fftshift(I);Ir=real(Is);figure,imshow(Ir)newheadraw=heaclraw;fori=l:2:256%奇数彳亍forj=l:256newheadraw(i,j)=0;endendI=ifft2(newheadraw);Is=fftshift(I);Ir=real(Is);figure,imshow(Ir)
8、!—文件：D.m—>newheadraw=headraw;fori=2:2:256%偶数行forj=l:256newheadraw(i,j)=0;endI=ifft2(newheadraw);Is=fftshift(I);Ir=real(Is);figure,imshow(Ir)
