用MATLAB编写函数求y=1/1²+1/2²+……+1/n²

以下是MATLAB编写函数求y=1/1²+1/2²+……+1/n²的代码：

function y = sum_of_squares(n)
    y = 0;
    for i = 1:n
        y = y + 1/i^2;
    end
end

在MATLAB中，可以使用for循环来实现累加求和。在这个函数中，我们首先将y初始化为0，然后使用for循环从1到n遍历每个数，将1/i²加到y中。最后，返回y作为函数的输出。

相关问题

1.Matlab计算根据 π²/6=1/1²+1/2²+1/3²+……+1/n².求π的近似值。当n分别取100、1000、10000时,结果是多少要求:分别用循环结构和向量运算(使用sum函数)来实现。 ​2.Matlab计算a=1+2+3+……+100。

1. 使用循环结构实现：

n = 100;
sum1 = 0;
for i = 1:n
    sum1 = sum1 + 1/i^2;
end
approx_pi1 = sqrt(6 * sum1);

n = 1000;
sum2 = 0;
for i = 1:n
    sum2 = sum2 + 1/i^2;
end
approx_pi2 = sqrt(6 * sum2);

n = 10000;
sum3 = 0;
for i = 1:n
    sum3 = sum3 + 1/i^2;
end
approx_pi3 = sqrt(6 * sum3);

使用向量运算（使用sum函数）实现：

n = 100;
sum1 = sum(1./(1:n).^2);
approx_pi1 = sqrt(6 * sum1);

n = 1000;
sum2 = sum(1./(1:n).^2);
approx_pi2 = sqrt(6 * sum2);

n = 10000;
sum3 = sum(1./(1:n).^2);
approx_pi3 = sqrt(6 * sum3);

2. Matlab计算a=1+2+3+……+100的结果可以使用等差数列求和公式来计算，也可以使用循环结构实现：

% 使用等差数列求和公式
n = 100;
a = n * (n + 1) / 2;

% 使用循环结构
n = 100;
sum_a = 0;
for i = 1:n
    sum_a = sum_a + i;
end

帮我写一个用matlab设计模拟滤波器系统函数为Ha(s)=3/(s²+4s+3)的AF的幅度特性的程序

好的，以下是用 MATLAB 设计模拟滤波器系统函数为 Ha(s) = 3/(s²+4s+3) 的 AF 幅度特性的程序：

% 定义系统函数 Ha(s)
num = [3];
den = [1 4 3];
Ha = tf(num, den);

% 绘制幅度特性曲线
w = logspace(-2, 2, 1000); % 频率范围
[mag, phase] = bode(Ha, w); % 计算幅度和相位
magdB = 20*log10(mag); % 将幅度转换为 dB
semilogx(w, magdB); % 绘制幅度特性曲线
grid on; % 添加网格线
title('AF幅度特性曲线'); % 添加标题
xlabel('频率 (rad/s)'); % 添加横坐标标签
ylabel('幅度 (dB)'); % 添加纵坐标标签

执行以上程序，即可得到 AF 幅度特性曲线。