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