MATLAB新手实验(4)--函数的编写及调试

MATLAB新手实验(4)–函数的编写及调试

一、实验目的：

掌握MATLAB函数的编写及调试方法。

二、实验内容：

函数1

1.在 MATLAB 里，建议用 .m 文件来写函数，编写函数（保存为myfun1.m）。

%这是matlab标准函数定义 x可以是标量或矢量
function [y] = myfun1(x)

y = zeros(size(x)); % 预分配结果，保持和x同大小

for i = 1 : length(x)
    if x(i) <= 0
        y(i) = sin(x(i));
    elseif x(i) <= 3
        y(i) = x(i);
    else
        y(i) = -x(i) + 6;
    end
end
end

2. 测试你的函数
在命令窗口或新建脚本里运行：

x = [-2, 0, 1, 3, 4];   % 随便选几个测试点
y = myfun1(x)

测试成功：

函数2

给出myfun2的代码段：

function [m_x, max_x, min_x, rms_x] = myfun2(x)
%计算均值、最大值、最小值、均方根值
m_x = mean(x);            % 平均值
max_x = max(x);           % 最大值
min_x = min(x);           % 最小值
rms_x = sqrt(mean(x.^2)); % 均方根值（RMS）
end

给出test2的测试代码：

x = sin(0:0.01:6*pi);
[m, mx, mn, rmsv] = myfun2(x);
fprintf('均值: %.4f, 最大值: %.4f, 最小值: %.4f, 均方根值: %.4f\n', m, mx, mn, rmsv);
%%
x = rand(1,200);
[m, mx, mn, rmsv] = myfun2(x);
fprintf('均值: %.4f, 最大值: %.4f, 最小值: %.4f, 均方根值: %.4f\n', m, mx, mn, rmsv);

测试结果：

函数3

这边是要写一个函数，输入一个行向量，输出一个矩阵，每一列是xi的0到n-1次幂，给出myvander.m的代码：

function v = myvander(x)
% 生成范德蒙德矩阵
n = length(x);      % x的元素个数
v = zeros(n);       % 生成n×n的零矩阵

for i = 1:n
    v(:,i) = x(i).^(0:n-1)';   % 每一列是x(i)的0~n-1次幂（竖着）
end

% 或等价写法：每行是幂次
% for k = 1:n
%     v(k,:) = x.^(k-1);
% end
end

再给出test测试的代码：

v = myvander([2 3 4 5])

测试结果如下：

函数4

根据公式迭代编写代码sqrt_iter_custom如下给出：

%迭代结果 迭代次数 误差
function [root, iter, err] = sqrt_iter_custom(a)
% 用迭代法计算sqrt(a)
x_old = 1.0;%迭代初值
tol = 1e-5;%终止条件
%迭代次数不允许超过100次
max_iter = 100;

for iter = 1:max_iter
    x_new = x_old/2 + a/(2*x_old);
    if abs(x_new - x_old) <= tol
        break;
    end
    x_old = x_new;
end
root = x_new;
err = abs(root - sqrt(a));
end

测试代码如下，写个循环直接全测试了：

a_list = [3, 17, 113];
for k = 1:length(a_list)
    a = a_list(k);
    [root, iter, err] = sqrt_iter_custom(a);
    fprintf('a = %d:\n', a);
    fprintf('  迭代结果 = %.8f\n', root);
    fprintf('  精确值   = %.8f\n', sqrt(a));
    fprintf('  误差     = %.2e\n', err);
    fprintf('  迭代次数 = %d\n\n', iter);
end

测试结果：