MATLAB 第四章习题
1.绘制函数曲线y=2sin(3πt+π/4)，t的范围是0~2.
2.在同一图形窗口绘制曲线y1=sin(t),t的范围为0~4π；y2=2cos(2t),t的范围为π~3π。要求y1曲线为黑色点画线，y2曲线为红色虚线圆圈，并在图的右下角标注2条曲线的图例（用legend），横坐标以π为单位标注，运行界面如图所示。
3.在同一图形窗口中绘制函数y=2^1/2e^(-t)sin(2πt+π/4)及其包络线图形，t的范围为【0,2】。
其运行结果如图所示。
4.在同一图形窗口中分别绘制y1=sin(2πt)、y2=cos(2πt)、y3=e^(-4t)共三条函数曲线，t的范围为[0,2].给坐标轴加上标注，给整个图形加上标题，在图形窗口添加文本字符串用于对各曲线分别加以文字说明，其运行界面如图所示。
5.绘z=xe^(-x^2+y^2)的三维曲面图和曲线图，x的范围为[-2,2],y的范围为[-2,2]，并实现部分镂空。
6.绘制z=x^2+y^2的三维网线图和曲面图，x的范围为[-5,5],y的范围为[-5,5].将网线图用spring色图并用颜色标尺显示色图，将曲面图颜色用sshading命令做连续变化。
题目： ）
选项A：对
选项B：错
答案：对
题目：Matlab中注释行的开头是#。（ ）
选项A：对
选项B：错
答案：错
题目：已知程序如下，关于该程序的功能描述正确的是（ ）。
t=0:0.5:10;
y1=exp(-0.1*t);
y2=exp(-0.2*t);
y3=exp(-0.5*t);
plot(t,y1,'-ob',t,y2,':*r',t,y3,'-.^g')
选项A：a=0.1时为蓝色-o-o-o-o线
选项B：a=0.2时为绿色..* ..* ..* ..*线
选项C：a=0.5时为绿色-.^-.^-.^-.^线
选项D：以上程序是用不同的线型和标记点画出三种情况下的曲线
选项E：有一组测量数据满足 以上程序是用不同的线型和标记点画出三种情况下的曲线,
a=0.1时为蓝色-o-o-o-o线, a=0.5时为绿色-.^-.^-.^-.^线
第四章
2.求方程xln(+x)--0.5x=0的正根；
fun=inline('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x');
fzero(fun,2)
答：所求正根为2.1155。
8.作出下列函数图形，观察所有局部极大，局部极小和全局最大，全局最小值点的粗略位置，并用MATLB函数fminbnd和fminsearch求各极值点的确切位置。
（1）f(x)=x2sin(x2-x-2); [-2,2]
(2) f(x)=3x5-20x3+10; [-3,3]
(3) f(x)= ; [0,3]
解：（1）clear;
fun=inline('x.^2.*sin(x.^2-x-2)');
fplot(fun,[-2 2]);grid on;
x(1)=-2;
x(3)=fminbnd(fun,-1,-0.5);
x(5)=fminbnd(fun,1,2);
fun2=inline('-x.^2.*sin(x.^2-x-2)');
x(2)=fminbnd(fun2,-2,-1);
x(4)=fminbnd(fun2,-0.5,0.5);
x(6)=2
feval(fun,x)
（2）clear;
fun=inline('3*x.^5-20*x.^3+10');
fplot(fun,[-3 3]);grid on;
x(1)=-3;
x(3)=fminsearch(fun,2.5);
fun2=inline('-(3*x.^5-20*x.^3+10)');
x(2)=fminsearch(fun2,-2.5);
x(4)=3;
feval(fun,x)
（3）fun=inline('abs(x^3-x^2-x-2)');
fplot(fun,[0 3]);grid on;
fminbnd(fun,1.5,2.5)
fun2=inline('-abs(x^3-x^2-x-2)');
第三章
1．要求在闭区间]2,0[上产生具有10个等距采样点的一维数组。试用两种不同的指令实
现。
方法一t1=linspace(0,2*pi,10)
t1=00.69811.39632.09442.79253.49074.18884.88695.58516.2832
方法二t2=0:2*pi/9:2*pi
结果与t1相同
2．由指令rand('state',0),A=rand(3,5)生成二维数组A，试求该数组中所有大于0.5的元素
的位置，分别求出它们的“全下标”和“单下标”。rand('state',0),A=rand(3,5)
[ri,cj]=find(A>0.5);
id=sub2ind(size(A),ri,cj);
ri=ri';
cj=cj';
disp('')
disp('大于0.5的元素的全下标')
disp(['行号',int2str(ri)])
disp(['列号',int2str(cj)])
disp('')
disp('大于0.5的元素的单下标')
disp(id')
A=
0.95010.48600.45650.44470.9218
0.23110.89130.01850.61540.7382
0.60680.76210.82140.79190.1763
大于0.5的元素的全下标
行号132332312
列号112234455
大于0.5的元素的单下标13569111213143．已知矩阵4321A，运行指令B1=A.^(0.5),B2=A^(0.5),可以观察到不同运算方法
所得结果不同。（1）请分别写出根据B1,B2恢复原矩阵A的程序。（2）用指令检验
所得的两个恢复矩阵是否相等。A=[1,2;3,4];
B1=A.^0.5
B2=A^0.5
A1=B1.*B1;
A2=B2*B2;
norm(A1-A2,'fro')
B1=
1.00001.4142
1.73212.0000
B2=
0.5537+0.4644i0.8070-0.2124i
