MR原理、技术与设备（冯衍秋）模拟Bloch方程的受激（Excitation）过程.docxVIP

模拟Bloch方程的受激(Excitation)过程。答：模拟的MATLAB脚本内容如下，共1个!-- 文件：excitation.m --dT=0.5; %间隔1msT=500; %总时间N=ceil(T/dT)+1;time=[0:0.1:N-1]*dT;df = 10;w1=2*pi*df/10000;w2=2*pi*df/500;z=cos(w1.*time);x=sin(w1.*time).*sin(w2.*time);y=sin(w1.*time).*cos(w2.*time);plot3(x,y,z);grid onwhile 1comet3(x,y,z);pause(0.2)end!—- 文件结束--结果如下三图（动画截图）：图一：图二：图三：模拟Bloch方程的弛豫(Relaxation)过程。答：模拟MATLAB脚本如下，共3个：!-- 文件：relaxation.m -- dT = 0.05; %间隔1ms T= 1000; %总时间 N = ceil(T/dT)+1; %时间间隔数 df = 10; %0.1s后关闭共振场 T1 = 600; T2 = 100; %生成适当的种子矩阵 [A B] = freeprecess(dT,T1,T2,df); time = [0:0.1:N-1]*dT;%===== 模拟衰减===== M = zeros(3,N); %实时跟踪磁化矢量末端位置 M(:,1)=[1;0;0];%加磁场 for k=2:N M(:,k) = A*M(:,k-1)+B; end;% ===== 绘制结果 ====== plot3(M(1,:),M(2,:),M(3,:)) grid onwhile 1comet3(M(1,:),M(2,:),M(3,:),0.01);pause(0.5);end!—- 文件结束 --!-- 文件：freeprecess.m --function [Afp,Bfp]=freeprecess(T,T1,T2,df)%%模拟自由进动和消共振的函数%经过 T时间, 给定弛豫时间T1、T2%和消共振频率df. 时间单位：ms, 消共振频率：Hz.phi = 2*pi*df*T/1000; %弛豫进动半径E1 = exp(-T/T1);E2 = exp(-T/T2);Afp = [E2 0 0;0 E2 0;0 0 E1]*zrot(phi);Bfp = [0 0 1-E1];!—- 文件结束 --!-- 文件：zrot.m --function Rz=zrot(phi)% %这是Z轴旋转变换函数Rz = [cos(phi) -sin(phi) 0;sin(phi) cos(phi) 0; 0 0 1]!—- 文件结束 --结果如下三图（动画截图）：图一：图二：图三：给定频谱图像headraw，按下列要求进行傅里叶逆变换重建图像：将数据第1到64行数据置0将数据1到64行与193到256行置0选择数据奇数行选择数据偶数行将数据放到512*512大小的矩阵中心，其他位置处为0答：各条件下代码分列如下，附结果图像：!-- 文件：A.m --for i=1:64 for j=1:256 newheadraw(i,j)=0; endendfor i=65:256 for j=1:256 newheadraw(i,j)=headraw(i,j); endendI=ifft2(newheadraw);Is=fftshift(I);Ir=real(Is);figure,imshow(Ir)!—- 文件结束 --结果图：!-- 文件：B.m --for i=1:64 for j=1:256 newheadraw(i,j)=0; endendfor i=65:192 for j=1:256 newheadraw(i,j)=headraw(i,j); endendfor i=193:256 for j=1:256 newheadraw(i,j)=0; endendI=ifft2(newheadraw);Is=fftshift(I);Ir=real(Is);figure,imshow(Ir)!—- 文件结束 --结果图：!-- 文件：C.m --newheadraw=headraw;for i=1:2:256 %奇数行 for j=1:256 newheadraw(i,j)=0; endendI=ifft2(newheadraw);Is=fftshift(I);Ir=real(Is);figure,imshow