Matlab实验三答案

function y=ywjz(A)
x=0;
n=numel(A);
for i=1:n
x=x+A(i);
end
v=x/n;
运行结果：
>>y=ywjz([1 2 3 4 5])
y=
3
编程2：
function y=ewjz(A)
x=0;
n=numel(A);
for i=1:n
x=x+A(i);
end
y=x/n;
if x==0
A(j)=y(i);j=j+1;
end
if x==1
B(k)=y(i);k=k+1;
end
if x==2
C(m)=y(i);m=m+1;
end
end
disp(A)
disp(B)
disp(C)
运行结果：
>> yush
3 9
7 16
5 11
7、已知一维数组A=[2,4,5,8,10]、B=[4,9,6,7,4]，用for循环语句实现 。求和函数可用sum().
2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位置加break语句、以便程序能正常执行。（执行循环体的次数不确定；每一次执行循环体后，一定会改变while后面所跟关系式的值。）
3、while循环也可以嵌套、其结构如下：
while逻辑表达式1
循环体语句1
while逻辑表达式2
循环体语句2
编程：
A=[2 4 5 8 10;];
B=[4 9 6 7 4;];
C=[];
y=0;
for i=1:5
C(i)=(A(i)*B(5-i+1));
y=y+C(i);
end
y
运行结果：
>> timu7
y =
178
8、（1）编写求解一维数组平均值的函数文件。
（2）编写求解二维数组平均值的函数文件。
编程1：
程序流程控制指令
break、return、pause
2、函数文件的结构和编写方法
函数只能在函数体内对变量进行操作，也就是只能访问函数本身工作空间中的变量。
M函数文件的结构
（1）函数定义行（function）
（2）H1行（函数帮助文本的第一行）
（3）函数帮助文本
（4）函数体
（5）注释
函数文件编写后，保存时，其文件名必须与函数名相同。注意：函数名不要与MATLAB自身的函数命令相同。
end
循环体语句3
end
if-else-end分支结构
if表达式1
语句1
else if表达式2(可选)
语句2
else(可选)
语句3
end
end
说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。
2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。
y =
1.9826e+087
y =
1.2689e+089
y =
8.2477e+090
y =
5.4434e+092
y =
3.6471e+094
y =
2.4800e+096
y =
1.7112e+098
y =
1.1979e+100
i =
70
5、编写一个乘法表（9×9）
编程：
for i=1:9
A=[];
for j=1:i;
y =
3.3453e+049
y =
1.4050e+051
y =
6.0415e+052
y =
2.6583e+054
y =
1.1962e+056
y =
5.5026e+057
y =
2.5862e+059
y =
1.2414e+061
y =
6.0828e+062
y =
3.0414e+064
y =
1.5511e+066
三、实验内容：
2、用for循环语句实现求1～100的和
编程：
sum=0;
for i=1:100
sum=sum+i;
end
sum
运行结果：
>> ljia
sum =
5050
3、用for循环语句实现编写一个求阶乘的函数文件。
编程：
function y=jcheng(n)
y=1;
for i=1:n
y=y*i;
A=[A,sprintf('%dx%d=%d ',j,i,i*j)];
end
disp(A)
end
运行结果：
>> cfb
1x1=1
1x2=2 2x2=4
1x3=3 2x3=6 3x3=9
1x4=4 2x4=8 3x4=12 4x4=16
1x5=5 2x5=10 3x5=15 4x5=20 5x5=25
运行结果：
>> y=ewjz([2 3;4 5;6 7;])
y =
4.5000
调用验证：
A=[1 3 3;5 5 7;];
y=ewjz(A);
y
>> aa
y =
4
9、编写验证魔方矩阵的函数文件，输出要求如下：
（1）如果输入矩阵的维数小于3，输出显示’error’
（2）如果输入矩阵的不是方阵，输出显示’the size of matrix X must be N-by-N matrix’
end
y
运行结果：
假设n=5
>> jcheng(5)
y =
120
ans =
120
4、找到一个n!>10100的值（利用上题的阶乘函数文件）
编程：
i=1;
while(1)
y=jcheng(i);
if y>10^100
break;
end
i=i+1;
end
i
运行结果：
>> zhaon
y =
1
y =
2
y =
实验三MATLAB语言的程序设计
一、实验目的：
1、熟悉MATLAB程序编辑与设计环境
2、掌握各种编程语句语法规则及程序设计方法
3、函数文件的编写和设计
4、了解和熟悉跨空间变量传递和赋值
二、实验基本知识：
1、程序流程控制语句
for循环结构
语法：for i＝初值：增量：终值
语句1
……
语句n
end
说明：1．i＝初值：终值，则增量为1。
2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。
while循环结构
语法：while逻辑表达式
循环体语句
end
说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达式的值为假，则程序执行end之后的语句。
6、判断y向量的元素属于3n、3n+1或3n+2，设y= [3,5,7,9,11,16]。
编程：
假设：将3n放进数组A，将3n+1放进数组B，将3n+2放进数组C。
y=[3 5 7 9 11 16;];
A=[];
B=[];
C=[];
j=1;k=1;m=1;
for i=1:6
x=rem(y(i),3);
y =
2.6525e+032
y =
8.2228e+033
y =
2.6313e+035
y =
8.6833e+036
y =
2.9523e+038
y =
1.0333e+040
y =
3.7199e+041
y =
1.3764e+043
y =
5.2302e+044
y =
2.0398e+046
y =
8.1592e+047
disp('eroor')
end
end
S=[sum(x,1),sum(x,2)',sum(diag(x))',sum(diag(rot90(x)))'];
if S(1)==S(2)&S(2)==S(3)&S(3)==S(4)&S(1)==S(4)
disp('yes')
else
disp('no')
end
验证运行：
>> yanzh([8 1 6;3 5 7;4 9 2;])
yes
>> yanzh([8 1 6;3 5 7;])
the size of matrix X must be N-by-N matrix
no
>> yanzh([8 6;3 5;])
eroor
no
>> yanzh([8 2 6;3 5 7;4 9 2;])
（3）显示行、列和及其对角线求和后的值，并判断其和是否相同，若不同，显示‘No’,相同显示‘Yes’
编程：
function yanzh(x)