湖 北 汽 车 工 业 学 院 《MATLAB及应用》实验指导书
plot(scalar_x,y_linear),title('method=linear'); subplot(2,2,3),plot(x,y,'*'),hold on,
plot(scalar_x,y_spline),title('method=spline'); subplot(2,2,4),plot(x,y,'*'),hold on, plot(scalar_x,y_cubic),title('method=cubic');
2. 求下列函数的解，并绘制图形。
x5(1) y?e?x，初始点为x?8
(2) y?xsinx
clear,clc %第一小题
y_1=@(x)exp(x)-x^5; x0 =8;
subplot(1,2,1), hold on,fplot(y_1,[x0,x0+10]),title('exp(x)-x^5'); %第二小题
y_2=@(x)x*sin(x);
subplot(1,2,2), hold on,fplot(y_2,[-pi,pi]),title('x*sin(x)');
9
湖 北 汽 车 工 业 学 院 《MATLAB及应用》实验指导书
3. 求下列函数的极值。
(1) (2)
z?x2??y?1?z??x?y?1?22
clear,clc
z_1 = @(x)x(1)^2-(x(2)-1)^2;
[x,fvalue,flag,output]=fminsearch(z_1,[0,0]) disp('第二小题')
z_2 = @(x)(x(1)-x(2)+1)^2;
[x,fvalue,flag,output]=fminsearch(z_2,[0,0])
4. 计算下列积分。
(1) ? 1 ?1x?x3?x5dx
(2)
?? 110101sinyx?ydxdyx2?4
10
湖 北 汽 车 工 业 学 院 《MATLAB及应用》实验指导书
clear,clc
fun1=@(x)x+x.^3+x.^5; q = quad(fun1,-1,1)
fun2=@(x,y)sin(y).*((x+y)./(x.^2+4)); q = dblquad(fun2,1,10,1,10)
第八章
1. 编写程序，该程序在同一窗口中绘制函数在
?0,2??之间的正弦曲线和余弦曲线，
步长为 ?/10，
线宽为 4 个象素，正弦曲线设置为蓝色实线，余弦曲线颜色设置为红色虚线，两条曲线交点处，用红色星号标记。
clear,clc x=0:pi/10:2*pi; f=@(x)(cos(x)-sin(x)); x1=fzero(f,[0,pi]); x2=fzero(f,[pi,2*pi]);
plot(x,sin(x),'b-','LineWidth',4),hold on,plot(x,cos(x),'r:','LineWidth',4); plot(x1,sin(x1),'rh','markerfacecolor','y','markersize',10); plot(x2,sin(x2),'rh','markerfacecolor','y','markersize',10);
11
湖 北 汽 车 工 业 学 院 《MATLAB及应用》实验指导书
2． 绘制下列图像
（1）y?xsinx，0?x?10?
22z?x?6xy?y?6x?2y?1，?10?x?10，?10?y?10 （2）三维曲线：
x2y2z??164，?16?x?16，?4?y?4 （3）双曲抛物面：
clear,clc
x = 0:pi/100:10*pi; y=x.*sin(x);
subplot(1,3,1),plot(x,y,'b'),title('y=x*sinx') %μú￡¨2￡?D?ìa
[X,Y] = meshgrid(-10:0.2:10);
Z = X.^2 + 6*X*Y + Y.^2 + 6*X + 2*Y -1; subplot(1,3,2),mesh(X,Y,Z),title('三维曲面') %μú￡¨3￡?D?ìa
12
湖 北 汽 车 工 业 学 院 《MATLAB及应用》实验指导书
[X,Y] = meshgrid(-16:0.2:16,-4:0.1:4); Z = X.^2/16 - Y.^2/4;
subplot(1,3,3),mesh(X,Y,Z),title('双曲线抛物面')
3. 绘制下列图像
（1）绘制电脑磁盘使用情况的饼状图
（2）生成 100 个从 0 到 10 之间的随机整数，绘制其直方图 （3）生成 10个从 0 到 10 之间的随机整数，绘制其阶跃图
clear,clc x=[37 63];
subplot(1,3,1),pie(x,{'可用空间37%','已用空间63%'}),title('饼状图'); subplot(1,3,2),hist(round(rand(100,1)*10)),title('直方图'); subplot(1,3,3),stairs(round(rand(10,1)*10)),title('阶跃图');
13