我审查完公式,又从新编辑,希望你能再提点一下
我的想法是 解四个方程,a,b是参数,t(1)_t(4)为变量,t(1)_t(4)分别随a,b变化,生成三维图,
t0=[0;7;0];
a=0:5;
b=0:5;
[A,B]=meshgrid(a,b)
p1=zeros(6,6);
p2=zeros(6,6);
p3=zeros(6,6);
p4=zeros(6,6);
for i=1:length(a)
for j=1:length(b)
p(:,i)=fsolve(@(t)SH_2(t,a(i),b(j)),t0);
p1(i,j)=real(sol(1));
p2(i,j)=real(sol(2));
p3(i,j)=real(sol(3));
p4(i,j)=real(sol(4));
end
end
figure
subplot(2,2,1)
surf(A,B,p1); shading interp; colorbar
subplot(2,2,2)
surf(A,B,p2); shading interp; colorbar
subplot(2,2,3)
surf(A,B,p3); shading interp; colorbar
subplot(2,2,4)
surf(A,B,p4); shading interp; colorbar
func