以下是修改后的代码,应该能够正常运行:

clear all clc

% 定义自变量范围 eta = 10e-5:10:10e3;

% 定义初始向量
alpha=[0.1,0.3,0.5,0.7,0.9];

% 计算函数值 for i=1:5 r(i)=((2-alpha(i))/alpha(i))^(1/2); K(i,:) = sin(alpha(i).pi)/pieta.^alpha(i)./(eta.^2+r(i).^2); end

% 绘图 figure; for i=1:5 plot(eta, K(i,:)); hold on; end

xlabel('eta'); ylabel('K'); legend('alpha=0.1', 'alpha=0.3', 'alpha=0.5', 'alpha=0.7', 'alpha=0.9'); title('K vs eta for different alpha values');


原文地址: https://www.cveoy.top/t/topic/bpiO 著作权归作者所有。请勿转载和采集!

免费AI点我,无需注册和登录