clear all;clc;
global hx0 hy0 hz0 % 外加磁场参数
hz0=0; hy0=0;
i=0;
mz3=0;
mz4=0;
for hx0=0:0.01:5;
y0=[0.11;0.01;pi-0.11;pi+0.01];
t0=0:0.1:500;
[t,y]=ode45('rhssys',t0,y0);
mx1=sin(y(:,1)).*cos(y(:,2));
my1=sin(y(:,1)).*sin(y(:,2));
mz1=cos(y(:,1));
mx2=sin(y(:,3)).*cos(y(:,4));
my2=sin(y(:,3)).*sin(y(:,4));
mz2=cos(y(:,3));
i=i+1;
mz3(i)=cos(y(5000,1));
mz4(i)=cos(y(5000,3));
end
hx0=0:0.01:5;
subplot(2,1,1);
plot(hx0,mz3,'LineWidth',2);
xlabel('h_x0','fontsize',20);ylabel('m_z1','fontsize',20);
title('eta=0.33','fontsize',20);
set(gca,'FontSize',20);
subplot(2,1,2);
plot(hx0,mz4,'LineWidth',2);
xlabel('h_x0','fontsize',20);ylabel('m_z2','fontsize',20);
set(gca,'FontSize',20);
此循环程序现在实现了hx0和t0的双循环吗?即每一个hx0,此程序算了t0从0到500所对应的mx,my,mz吗?如果没有,如何实现hx0和t0的双循环?最后怎样改变语句,可以读出所有t0对应的y0值,而不是只是最后一组y0.