matlab_题库(编辑修改稿)

matlab_题库(编辑修改稿)内容摘要：

yxz 。 ► x=linspace(1,1,100)。 %分割 [1,1]区间生成 x ► y=x。 %y 与 x相同 ► [X,Y]=meshgrid(x,y)。 %生成矩形域 [1,1] [1,1]网格节点坐标矩阵 ► Z=X.^2+Y.^2。 %生成 22 yxz  函数值矩阵 ► subplot(1,2,1) ► mesh(X,Y,Z)。 %在第一个子图中画 22 yxz  网格曲面 ► subplot(1,2,2) ► surf(X,Y,Z)。 %在第二个子图中画 22 yxz  光滑曲面 ► shading flat。 %对曲面 22 yxz  平滑并除去网格 例 28 在圆形域 122 yx 上绘制旋转抛物面： 22 yxz 。 ► x=linspace(1,1,300)。 %分割 [1,1]区间生成 x ► y=x。 %生成 y ► [X,Y]=meshgrid(x,y)。 %生成矩形域 [1,1]X[1,1]网格节点坐标矩阵 ► Z=X.^2+Y.^2。 %生成 22 yxz  函数值矩阵 ► i=find(Z1)。 %找出圆域 122 yx 之外的函数值（ z1）坐标点i ► Z(i)=NaN。 %对圆域 122 yx 之外的坐标点 i 处函数值进行“赋空” ► subplot(1,2,1) ► mesh(X,Y,Z)。 %在第一个子图中画 22 yxz  网格曲面 ► subplot(1,2,2) ► surf(X,Y,Z)。 %在第二个子图中画 22 yxz  光滑曲面 ► shading flat。 %对曲面 22 yxz  平滑并除去网格 例 29 画出 2222sinyxyxz在 ||,||  yx 上的图形。 ► x=::。 ► y=x。 ► [X,Y]=meshgrid(x,y)。 ► u=sqrt(X.^2+Y.^2)+eps。 %加 eps 使得 u 不等于 0，保证 z 有意义 ► Z=sin(u)./u。 ► surf(X,Y,Z) 例 210 有一组实验数据如下表所示，试绘图表示。 时 间 1 2 3 4 5 6 7 8 9 数据 1 数据 2 数据 3 ► t=1:9。 ► d1=[ ]。 ► d2=[ ]。 ► d3=[ ]。 ► plot(t,d1,39。 r+39。 ,t,d2,39。 kx:39。 ,t,d3,39。 b*39。 ,39。 linewidth39。 ,2,39。 markersize39。 ,8)。 ► title(39。 time amp。 data39。 )。 ► xlabel(39。 time39。 )。 ylabel(39。 data39。 )。 ► axis([0 10 0 100])。 ► text(,39。 \leftarrowdata139。 )。 % 39。 \leftarrow 39。 表示画一左箭头 ← ,且在标识前 ► text(3,39。 data2\rightarrow39。 )。 % 39。 \rightarrow 39。 表示画一右箭头 → ,且在标识后 ► text(,39。 \leftarrowdata339。 )。 ► grid 实验 3 MATLAB 编程介绍与循环结构 例 31： 求 n（ n=100）个奇数的和： s=1+3+5+…+(2 n1). clear。 clc。 %清除内存变量，清理命令窗口 n=100。 %赋值给定奇数的个数 s=0。 %设定存放和的变量 s并赋初值 0 for i=1:n %定义循环变量 i从 1到 n，以 1为步长，即为奇数序号 s=s+(2*i1)。 %先计算右端奇数并累加后再赋给左端的变量 s fprintf(39。 i=%.0f, s=%.0f\n39。 ,i,s) %逐行显示出累加求和的过程 end %循环结构结束 例 32： 求正整数 n 的阶乘： p=1 2 3 … n = n!，并求出 n=20 时的结果。 clear。 clc。 %清除内存变量，清理命令窗口 n=20。 %赋值给定正整数 p=1。 %设定存放阶乘的变量 p并赋初值 1 for i=1:n %定义循环变量 i从 1到 n，以 1为步长，即连续正整数 p=p*i。 %先计算右端乘积后再赋给左端的变量 p fprintf(39。 i=%.0f, p=%.0f\n39。 ,i,p) %逐行显示出 i! end %循环结构结束 例 33： 根据麦克劳林公式可以得到 e≈ 1+1+1/2!+1/3!+„ +1/n!，试求 e的近似值。 clear。 clc。 %清除内存变量，清理命令窗口 n=10。 %赋值给定正整数 p=1。 %设定存放阶乘的变量 p并赋初值 1 s=1。 %设定存放累加和的变量 s并赋初值 1 for i=1:n %定义循环变量 i从 1到 n，以 1为步长 p=p*i。 %先计算右端乘积后再赋给左 端的变量 p，此时 p为 i的阶乘 s=s+1/p。 %先计算右端阶乘倒数的累加后再赋给左端的变量 s fprintf(39。 i=%.0f, s=%.8f\n39。 ,i,s) %逐行显示出第 i次 e的近似值 end %循环结构结束 例 34： 对于 数列   ,2,1, nn ，求其前 n项和不超过 1000时的 n的值及和 . clear。 clc。 %清除内存变量，清理命令窗口 n=0。 %设定正整数并赋初值 0 s=0。 %设定存放累加和的变量 s并赋初值 0 while s=1000 %用累加和 s与 1000进行比较作为循环条件 n=n+1。 %改变 n为连续正整数 s=s+sqrt(n)。 %先计算右端开方数的累加后再赋给左端的变量 s fprintf(39。 n=%.0f, s=%.4f\n39。 ,n,s) %逐行显 示正整数及部分和 end %循环结构结束 例 35： 根据 e≈ 1+1+1/2!+1/3!+„ +1/n! 求 e的近似值，要求精确到 810。 clear。 clc。 %清除内存变量，清理命令窗口 p=1。 %设定存放阶乘的变量 p并赋初值 1 s=1。 %设定存放累加和的变量 s并赋初值 1 r=1。 %设定前后两次近似 值的误差 r并赋初值 1 k=0。 %设定构造连续正整数的变量 k赋初值 0又为循环次数 while r= %当近似值的精度 r没达到 910 时继续循环 k=k+1。 %累计循环次数并作为下一个正整数 k p=p*k。 %计算 k的阶乘 p r=1/p。 %计算前后两次近似值的误差 r s=s+r。 %计算 e的近似值 s fprintf(39。 k=%.0f, s=%.10f\n39。 ,k,s) %逐行显示出第 k次 e的近似值 s end %循环结构结束 实验 4 MATLAB 选择结构与应用实验 例 41： 求任意有限数组 a=[a(1),a(2),„ ,a(n)] 中数值最大的元素 M 以及所在位置 k. function [M,k]=findM(a) %定义函数 findM，输入数组 a，返回最大元素 M及 位置 k n=length(a)。 %获取数组的长度即元素的个数 n M=a(1)。 k=1。 %将第一个元素作为最大值赋值给 M，位置为 1。 for i=2:n %从第二个元素到最后一个元素依次进行 if a(i)M %比较后续元素与目前最大值 M的大小 M=a(i)。 k=i。 %将数值较大的元素赋值给 M，同时保留位置 i end %选择结构结束 end %循环结构结束 ► a=[1,pi,0]。 %任意给定一数组 ► [M,k]=findM(a) %调用函数 findM ◄ M = ◄ k =5 例 42： 编写一个函数将百分制成绩转换为优 (A)，良 (B)，中 (C)，差 (D)四等级 . function jb=dengji(fs) %定义函数 dengji，输入分数 fs，返回等级 A,B,C,D if fs=90 %判断分数 fs是否处在优秀级别上 jb=39。 A 39。 %定义为 A级 elseif fs=78 %判断分数 fs是否处在良好级别上 jb=39。 B 39。 %定义为 B级 elseif fs=60 %判断分数 fs是否处在合格级别上 jb=39。 C 39。 %定义为 C级 else %分数 fs不处于以上任何级别上 jb=39。 D 39。 %定义为 D级 end %选择结构结 束 ► jb=dengji(81) %调用函数 dengji ◄ jb =B 例 43： Fibonacci 数组的元素满足 Fibonacci 规则： }{na ： 121 aa ， 12   kkk aaa , ,3,2,1k . 求出该数组中第一个大于 10000 的元素。 n=100。 %给定一个较大的 n作为数列的位置 a=[1,1]。 %设定数列的初始值 for i=3:n %从第 3个元素开始循环递推生成后续元素 p=a(i1)+a(i2)。 %前两个元素之和生成后续元素 p a=[a,p]。 %将刚产生的元素 p放置到数组 a的最后，拼接成新的数组 if p10000 %判断将刚产生的元素 p是否超过 10000 break。 %跳出所在的 for循环 end %选择结构结束 end %循环结构结束 disp([a]) %显示所生成的数列，最后一个元素。