clear all clc 定义自变量范围 eta = 10e-51010e3; 计算函数值 for i=15 alphai=0103050709; ri=2-alphaialphai^12; Ki = sinalphaipipieta^alphaieta^2+ri^2; ploteta K; hold on; end 添加图例和标题 legendalpha=01 alpha=03 al

该代码错误在于计算 r(i) 的公式有误，应该加括号，即 r(i) = (2 - alpha(i)/alpha(i))^(1/2)。

修改后的代码如下：

clear all clc

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

% 计算函数值 for i = 1:5 alpha(i) = [0.1, 0.3, 0.5, 0.7, 0.9]; r(i) = (2 - alpha(i)/alpha(i))^(1/2); K(i) = sin(alpha(i)*pi)/pi * eta^alpha(i) / (eta^2 + r(i)^2); plot(eta, K); hold on; end

% 添加图例和标题 legend('alpha=0.1', 'alpha=0.3', 'alpha=0.5', 'alpha=0.7', 'alpha=0.9'); title('K(alpha,eta) for different values of alpha∈(0,1)'); xlabel('eta'); ylabel('K(alpha,eta)');