主程序
clc
clear
n=1;
Ars=zeros(1,4);
x=linspace(1,1.4,4);
for D=x
Ars(n)=hanshu(D);
n=n+1;
end
Ar=Ars;
figure(1)
D=linspace(1,1.4,4);
hold on
plot (D,Ar,'*-')
xlabel ('D')
ylabel ('Ar')
子程序
function [Ars]=hanshu(D)
syms Ar;
g=1*10^(-10);
Aa=0.0112
G=g/((Aa).^0.5)
k=1;
Q=0.01;
ac=G.^2./((k*Q/2).^(2./(D-1)));
g1=D.*((2-D)./D).^(D./2)./(3-2.*D);
g2=((2-D)./D).^((3-2.*D)./2);
g3=(D./(2-D)).^((2-D)./2)
P=8*10.^(-6);
equ=P==(1-3.*pi./16).*pi.*G.^(D-1).*g1.*Ar.^(D./2).*(g2.*Ar.^((3-2.*D)./2)-ac.^((3-2*D)./2))+k.*Q.*g3.*Ar.^(D./2).*ac.^((2-D)./2);
Ars=vpasolve(equ,Ar);