case num2cell(7:10) y=log(abs(b+c/x)); end
disp(['y=',num2str(y)]); 3. 答:
A=fix(10+89*rand(1,20)); sum=0; for i=1:20
sum=sum+A(i);
end
B=A(find(A<(sum/20))); C=B(find(rem(B,2)==0)); disp(C);
4. 答:
(1)
用循环结构实现:
输入20个数,求其中最大数和最小数。要求分别用循环结构和调用MATLAB的max函数、min函数来实现。 产生20个两位随机整数,输出其中小于平均值的偶数。
v_max=0; v_min=0; for i=1:20
x=input(['请输入第', num2str(i), '数:']); if x> v_max
end; if x< v_min
end;
end
disp(['最大数为:', num2str(v_max)]); disp(['最小数为:', num2str(v_min)]); (2)
用max函数、min函数实现:
v_min=x; v_max=x;
for i=1:5 end
disp(['最大数为:', num2str(max(A))]); disp(['最小数为:', num2str(min(A))]);
A(i)=input(['请输入第', num2str(i), '数:']);
5. 答:
已知:
s?1?2?22?23???263,分别用循环结构和调用MATLAB的sum函数求s的值。
(1) 用循环结构实现:
s=0; for i=0:63 s=s+2^i; end s
(2) 调用sum函数实现: s=0:63; s=2.^s; sum(s)
6.
当n分别取100、1000、10000时,求下列各式的值。
11111??????(?1)n?1??(?ln2)
234n111?????(?) (2) 1?357411111?????n??(?) (3)
4166443(1) (4)
?2?2??4?4??6?6??(2n)?(2n)????????? ????????1?33?55?7(2n?1)?(2n?1)?????????2? 答:
要求分别用循环结构和向量运算(使用sum或prod函数)来实现。
(1) 用循环结构实现: sum=0; for k=1:100 end sum
使用sum函数:
sum=sum+(-1)^(k+1)/k;
x=[]; for k=1:10000 end
x=[x, (-1)^(k+1)/k];
sum(x)
(2) 用循环结构实现: sum=0; for k=1:100 end sum
sum=sum+(-1)^(k+1)/(2*k-1);
使用sum函数: x=[];
for k=1:100 end
(3) 用循环结构实现: sum=0;
for k=1:100 end sum
使用sum函数实现:
x=[]; for k=1:100 end sum(x)
(4) 用循环结构实现: t=1; for k=1:100 end t
使用prod函数实现: x=[]; for k=1:100 end prod(x)
7. 编写一个函数文件,求小于任意自然数n的斐波那契(Fibnacci)数列各项。斐波那契数列定义如下:
x=[x, ((2*k)*(2*k))/((2*k-1)*(2*k+1))]; t=t*(((2*k)*(2*k))/((2*k-1)*(2*k+1))); x=[x, 1/(4^k)]; sum=sum+1/(4^k);
sum(x)
x=[x, (-1)^(k+1)/(2*k-1)];
?f1?1, n?1??f2?1, n?2 ?f?f?f, n?2n?1n?2?n答:
function x=fibnacci(n) for i=1:n
if i<=2 else end
x(i)=x(i-1)+x(i-2); x(i)=1;
end
8. 答:
函数文件myfnc.m: function [x, y]= myfnc(A, B) try
x=A*B;
catch end
y=A.*B;
命令文件myexe.m: A=input('请输入矩阵A:'); B=input('请输入矩阵B:'); [x, y]=myfnc(A, B); if length(x)==0
display('两矩阵的维数不匹配,无法进行乘积运算!');
else
disp('矩阵A和矩阵B的乘积为:'); x
end
disp('矩阵A和矩阵B的点乘为:'); y
x=[];
编写一个函数文件,用于求两个矩阵的乘积和点乘,然后在命令文件中调用该函数。
9. 先用函数的递归调用定义一个函数文件求
?ii?1nm,然后调用该函数文件求
?k??k??2k?1k?1100501。 kk?110答:
函数文件myfnc.m: function sum=myfnc(n, m) if n<=1 else end
sum= myfnc (n-1, m)+n^m; sum=1;
在命令窗口中调用myfnc.m文件,计算
?k??k??2k?1k?1100501:
k?1k10
sum=myfnc(100, 1)+ myfnc(50, 2)+myfnc(10,-1)
10. 写出下列程序的输出结果。
①
s=0;
a=[12,13,14;15,16,17;18,19,20;21,22,23];
for k=a
for j=1:4 end
end s
答:执行结果为 ②
s=108
命令文件exe.m执行后的结果为: x =
4
y=
2
4
6
12
20
if rem(k(j),2)~=0 end
s=s+k(j);