MATLAB程序设计基础（三）

目录

1、实验目的：

2、实验内容

---

1、实验目的：

1）熟悉MATLAB程序编辑与设计环境；

2）掌握各种编程语言语法规则及程序设计方法；

3）熟悉函数文件的编写与设计；

4）了解和熟悉跨空间变量传递和赋值。

---

实验要求：

1）实验前，要认真预习实验任务，了解实验目的和实验内容；

2）实验时，每题要编写程序代码形成独立的M文件，并调试程序使其能正确运行；

3）实验后，按要求编写实验报告，源程序要有适当的注释，以提高程序的可读性。

2、实验内容

2-1）、编写脚本M文件，实现对下图矩阵指定区域的访问。要求将原矩阵赋给变量A，绿色区域矩阵为B，粉色区域矩阵为C，黄色区域矩阵为D。

A = [4 10 1 6 2; 8 2 9 4 7;7 5 7 1 5; 0 3 4 5 8];
B = A(5)
C = A([2:4],[2 3])
D = A(:,end)

---

 新建一个M函数文件，输入函数表达式，代码如下：

function y = func(x)
% 求解表达式
% 注意：x是矩阵，所以用 点乘 点除 点乘方！！！！！！
y = 1 ./ ( (x-2).^2 + 0.1 ) + 1 ./ ( (x-3).^4 + 0.01 );

end

命令文件中调用如下代码

x = rand(2, 2)
a = func(x)

效果如下：

---

 2-3）、编写求向量均值和中值的程序，并求出向量A=[1,2,3,4,5,6,7,8,9,10]的均值和中值。

编写M函数文件

function [average, med] = junzhizhongzhi(u) %主函数
%junzhizhongzhi使用内部函数求均值和中值
n = length(u);   %长度
average = mean(u, n);   %平均值
med = median(u, n);     %中值

function a = mean(v,n)      %子函数，计算均值
    a = sum(v) / n;
function m = median(v, n)   %子函数，计算中值
    w = sort(v);            %sort函数：按升序或降序对数组元素排序
    if rem(n,2)==1          %rem函数 ：求除法的后余数
        m = w((n+1)/2);
    else
        m = (w(n/2)+w(n/2+1))/2;
end

 调用函数

A=[1,2,3,4,5,6,7,8,9,10];
[zhongzhi, junzhi__] = junzhizhongzhi(A)    %[均值, 中值]

运行结果如下：

---

2-4）、已知一作匀速直线运动的物体的初速度v0，加速度a，求任意时刻t的速度和位移。要求构建函数f，并在命令窗口调用该函数求在v0=2、a=4、t=5时，物体的速度和移动的距离。

注：速度公式 v = v0+at  ,   位移公式 x = v0*t+1/2*a*t*t

 代码如下：

function [v, x] = V_S(v0, a, t)
%构建函数f,输入初速度V0，加速度a，任意时刻t，求解速度和位移
v = v0 + a*t;
x = v0*t + 1/2*a*t*t;
end

[v,x] = V_S(2, 4, 5)    %调用函数