
syms a b c d D L k d_L
format shortE
L=[0.0005 0.001 0.01];
d_L=0:0.1:10
for i=1:length(L)
for j=1:length(d_L)
D = 2*10^(-5) ; c = 1*10^(17) ; a = 5*(10^3) ;
b=L(i)^2/D ;
P(i,j)=-(a*b*c*(L(i)*a + (L(i)*a*(exp(-a*d) - exp(j)))/sinh(j) + 1))/(L(i)^2*a^2 - 1)
end
end
plot(d_L,P/1e20,'.-'),grid on
legend('L=0.0005','L=0.001','L=0.01')
xlabel('d/L');ylabel('P(×10^{20})')
每次输入plot那行就会显示错误使用plot
刚刚才学matlab,所以不太懂,希望大佬可以帮忙解答一下
