第5章MATLAB绘图习题答案.doc
MATLAB语言基础与应用(第二版)第5章 习题答案

第5章习题与答案5.1用矩阵三角分解方法解方程组123123123214453186920x x x x x x x x x +-=⎧⎪-+=⎨⎪+-=⎩ 解答:>>A=[2 1 -1;4 -1 3;6 9 -1] A =2 1 -1 4 -13 6 9 -1 >>b=[14 18 20]; b =14 18 20 >> [L, U, P]=lu(A) L =1.0000 0 0 0.6667 1.0000 0 0.3333 0.2857 1.0000 U =6.0000 9.0000 -1.0000 0 -7.0000 3.6667 0 0 -1.7143 P =0 0 1 0 1 0 1 0 0 >> y=backsub(L,P*b’) y =20.0000 4.6667 6.0000 >> x=backsub(U,y) x =6.5000 -2.5000 -3.5000 5.2 Cholesky 分解方法解方程组123121332352233127x x x x x x x ++=⎧⎪+=⎨⎪+=⎩ 解答:>> A=[3 2 3;2 2 0;3 0 12] A =3 2 32 2 03 0 12>> b=[5;3;7]b =537>> L=chol(A)L =1.7321 1.1547 1.73210 0.8165 -2.44950 0 1.7321>> y=backsub(L,b)y =-11.6871 15.7986 4.0415>> x=backsub(L',y)x =-6.7475 28.8917 49.93995.3解答:观察数据点图形>> x=0:0.5:2.5x =0 0.5000 1.0000 1.5000 2.0000 2.5000 >> y=[2.0 1.1 0.9 0.6 0.4 0.3]y =2.0000 1.1000 0.9000 0.6000 0.4000 0.3000 >> plot(x,y)图5.1 离散点分布示意图从图5.1观察数据点分布,用二次曲线拟合。
MATLAB课后习题集附标准答案

第2章MATLAB概论1、与其他计算机语言相比较,MA TLAB语言突出的特点是什么?答:起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富,可扩展性强.2、MA TLAB系统由那些部分组成?答:开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时,在选择组件窗口中哪些部分必须勾选,没有勾选的部分以后如何补安装?答:在安装MATLAB时,安装内容由选择组件窗口中各复选框是否被勾选来决定,可以根据自己的需要选择安装内容,但基本平台(即MATLAB选项)必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行,只是在选择组件时只勾选要补装的组件或工具箱即可.4、MATLAB操作桌面有几个窗口?如何使某个窗口脱离桌面成为独立窗口?又如何将脱离出去的窗口重新放置到桌面上?答:在MATLAB操作桌面上有五个窗口,在每个窗口的右下角有两个小按钮,一个是关闭窗口的Close 按钮,一个是可以使窗口称为独立的Undock按钮,点击Undock按钮就可以使该窗口脱离桌面称为独立窗口,在独立窗口的view菜单中选择Dock,菜单项就可以将独立的窗口重新防止的桌面上.5、如何启动M文件编辑/调试器?答:在操作桌面上选择“建立新文件”或“打开文件”操作时,M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.6、存储在工作空间中的数组能编辑吗?如何操作?答:存储在工作空间的数组可以通过数组编辑器进行编辑:在工作空间浏览器中双击要编辑的数组名打开数组编辑器,再选中要修改的数据单元,输入修改内容即可.7、命令历史窗口除了可以观察前面键入的命令外,还有什么用途?答:命令历史窗口除了用于查询以前键入的命令外,还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.8、如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件有什么区别?答:当前目录可以在当前目录浏览器窗口左上方的输入栏中设置,搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下,只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用,如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件,如果没有特别说明,数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径?答:(1)帮助浏览器:选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.(2)help命令:在命令窗口键入“help”命令可以列出帮助主题,键入“help 函数名”可以得到指定函数的在线帮助信息.(3)lookfor命令:在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.(4)模糊查询:输入命令的前几个字母,然后按Tab键,就可以列出所有以这几个字母开始的命令和函数.注意:lookfor和模糊查询查到的不是详细信息,通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375,并将其赋予变量a ? 答:在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法?各有什么优点? 答:(1)直接输入法,如a=[1 2 3 4],优点是输入方法方便简捷;(2)通过M 文件建立矩阵,该方法适用于建立尺寸较大的矩阵,并且易于修改; (3)由函数建立,如y=sin(x),可以由MATLAB 的内部函数建立一些特殊矩阵; (4)通过数据文件建立,该方法可以调用由其他软件产生数据.3.3在进行算术运算时,数组运算和矩阵运算各有什么要求?答:进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则,如矩阵a 与b 相乘(a*b )时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。
Get清风MATLAB教程a第5章习题解答张志涌

MATLAB教程2012a第5章习题解答-张志涌第5章 数据和函数的可视化习题5及解答1 椭圆的长、短轴2,4==b a ,用“小红点线〞画椭圆⎩⎨⎧==tb y ta x sin cos 。
〔参见图p5-1〕〖解答〗 clf a=4;b=2;t=0:pi/80:2*pi; x=a*cos(t); y=b*sin(t);plot(x,y,'r.','MarkerSize',15) axis equal xlabel('x') ylabel('y')shg-4-3-2-101234-3-2-1123xy2 根据表达式θρcos 1-=绘制如图p5-2的心脏线。
〔提示:采用极坐标绘线指令polar 〕〖解答〗 clftheta=0:pi/50:2*pi;rho=1-cos(theta);h=polar(theta,rho,'-r');%极坐标绘线指令。
h 是所画线的图柄。
set(h,'LineWidth',4) %利用set 设置h 图形对象的“线宽〞axis square %保证坐标的圆整性0.51 1.523021060240902701203001503301800ρ=1-cos θ3 A,B,C 三个城市上半年每个月的国民生产总值如见表p5.1。
试画出如图p5-3所示的三城市上半年每月生产总值的累计直方图。
表p5.1 各城市生产总值数据〔单位:亿元〕城市 1月 2月 3月 4月 5月 6月 A 170 120 180 200 190 220 B 120 100 110 180 170 180 C 70508010095120〖目的〗● 借助MATLAB 的帮助系统,学习直方图指令polar 的使用。
● bar 指令常用格式之一:bar(x,Y,'style') 。
x 是自变量列向量;Y 是与x 行数相同的矩阵,Y 的每一行被作为“一组〞数据;style 取stacked 时,同一组数据中每个元素对应的直方条被相互层叠。
第5章MATLAB绘图_习题答案

第5章MATLAB绘图习题5一、选择题1.如果x、y均为4×3矩阵,则执行plot(x,y)命令后在图形窗口中绘制()条曲线。
DA.12B.7C.4D.32.下列程序的运行结果是()。
Ax=0:pi/100:2*pi;forn=1:2:10plot(n*sin(x),n*cos(x))holdonendaxissquareA.5个同心圆B.5根平行线C.一根正弦曲线和一根余弦曲线D.5根正弦曲线和5根余弦曲线3.命令text(1,1,'{\alpha}+{\beta}')执行后,得到的标注效果是()。
C A.{\alpha}+{\beta}B.αβ}C.α+βD.αβ4.subplot(2,2,3)是指()的子图。
AA.两行两列的左下图B.两行两列的右下图C.两行两列的左上图D.两行两列的右上图x的曲线绘制成直线,应采用的绘图函数是()。
C5.要使函数y=2eA.polarB.semilogxC.semilogyD.loglog6.下列程序的运行结果是()。
B[x,y]=meshgrid(1:5);surf(x,y,5*ones(size(x)));A.z=x+y平面B.与xy平面平行的平面C.与xy平面垂直的平面D.z=5x平面7.下列函数中不能用于隐函数绘图的是()。
DA.ezmeshB.ezsurfC.ezplotD.plot38.下列程序运行后,看到的图形()。
Ct=0:pi/20:2*pi;[x,y]=meshgrid(-8:0.5:8);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);surf(x,y,z)view(0,90);axisequalA.像墨西哥帽子B.是空心的圆C.边界是正方形D.是实心的圆9.下列程序运行后得到的图形是()。
A[x,y]=meshgrid(-2:2);z=x+y;i=find(abs(x)<1&abs(y)<1);z(i)=NaN;surf(x,y,z);shadinginterpA.在一个正方形的正中心挖掉了一个小的正方形B.在一个正方形的正中心挖掉了一个小的长方形C.在一个正方形的上端挖掉了一个小的正方形D.在一个正方形的下端挖掉了一个小的正方形10.在使用MATLAB“绘图”选项卡中的命令按钮绘图之前,需要()。
MATLAB基础习题第五章习题答案

case{-1,2}
disp('您赢了');
end
%%
%是否继续判别模块
n=input('是否继续玩该游戏?否(0),是(1)');
while n~=0&n~=1
disp('您输入的不是正确数字,请正确输入');
n=input('是否继续玩该游戏?否(0),是(1)');
end
(3)程序要具有友好性
答:
%%
%该模块实现的功能是:求解六元线性方程组,方程组的系数矩阵由用户通过键盘输入;得到系数矩阵后给出方程的解;
clc;
clear all;
close all;
%%
%方程输入模块
x=inputdlg({'第一个方程系数','第二个','第三个','第四个','第五个','第六个'});
end
>>
请选择,石头(1),剪刀(2),布(3):1
您的选择是:
石头
电脑的选择是:
、剪
您赢了
是否继续玩该游戏?否(0),是(1)1
请选择,石头(1),剪刀(2),布(3):2
您的选择是:
、剪
电脑的选择是:
石头
您输了
是否继续玩该游戏?否(0),是(1)0
>>
5.编写一个日程提醒程序实现如下功能:
(1)如果当前时间为7-8点则提醒用户,该吃早饭了;
end
end
%%
%%显示模块
disp('90分以上的人有: Nhomakorabea)disp(x)
第5章MATLAB绘图99930资料

程序如下:
x=(0:pi/100:2*pi)';
y1=2*exp(-0.5*x)*[1,-1];
y2=2*exp(-0.5*x).*sin(2*pi*x);
plot(x,y1,'b:'); axis([0,2*pi,-2,2]); %设置坐标 hold on; %设置图形保持状态
plot(x,y2,'k');
参数的grid命令在两种状态之间进行切换。 box on/off命令控制是加还是不加边框线,不带
参数的box命令在两种状态之间进行切换。
例5.5 绘制分段函数曲线并添加图形标注。
程序如下:
x=linspace(0,10,100);
y=[];
for x0=x
if x0>=8
y=[y,1];
elseif x0>=6
subplot(4,4,8);
%选择4×4个区中的8号区
plot(x,ct);title('cotangent(x)');axis ([0,2*pi,-40,40]);
对图形窗口灵活分割。
5.1.3 绘制二维图形的其他函数 1. 其他形式的线性直角坐标图
在线性直角坐标系中,其他形式的图形有 条形图、阶梯图、杆图和填充图等,所 采用的函数分别是:
当前图形窗口,图形窗口原有图形将不复 存在。若希望在已存在的图形上再继续添 加新的图形,可使用图形保持命令hold。 hold on/off命令控制是保持原有图形还是刷 新原有图形,不带参数的hold命令在两种状 态之间进行切换。
5章习题答案matlab
张卫华 MATLAB课堂
例13:
分析以下程序, 并运行观察。
clf;x=3*pi*(-1:0.05:1);y=x; [X,Y]=meshgrid(x,y); R=sqrt(X.^2+Y.^2)+eps; Z=sin(R)./R; h=surf(X,Y,Z);colormap(jet); axis off n=12;mmm=moviein(n); for i=1:n rotate(h,[0 0 1],25); mmm(:,i)=getframe; end movie(mmm,5,10)
张卫华 MATLAB课堂
例5程序:
fplot('cos(tan(pi*x))',[ 0,1],1e-4)
张卫华 MATLAB课堂
例6:
绘制r=sin(t)cos(t)的极坐标图,并标记 数据点。
张卫华 MATLAB课堂例6 Nhomakorabea序:t=0:pi/50:2*pi; r=sin(t).*cos(t); polar(t,r,'-*');
title('y1=0.2e^{-0.5x}cos(4\pix) 和y2=2e^{0.5x}cos(4\pix)比较') text(2,2*exp(-1),'\fontsize{20}\bf \leftarrow x_{1}=2,y_{1}=2e^{-1}')
张卫华 MATLAB课堂
例5:
绘制f(x)=cos(tan(πx))的曲线
张卫华 MATLAB课堂
例11:
绘制柱形,剪切掉x、y小于零 的部分
张卫华 MATLAB课堂
例11程序:
t=linspace(0,2*pi,100); r=1-exp(-t/2).*cos(4*t); [X,Y,Z]=cylinder(r,60); ii=find(X<0&Y<0); Z(ii)=NaN; surf(X,Y,Z);colormap(spring), shading interp light('position',[-3,-1,3],'style','local')
第五章matlab习题
已知 h(t)=e-t,t≥0,取 N=64, 对 t 从 0~5 秒采样, 用 FFT 作快速傅里 叶变换,并绘制相应的振幅-频率图。
4.当矩阵 A 为自变量时,求 P(X)的值。其中 A 的值与题 3 中的相同
将 100 个学生的 5 门功课的成绩存入矩阵 P 中,进行以下处理 1. 2. . 4. 分别求每门课的最高分,最低分及相应学生序号 分别求每门课的平均分和标准方差。 5 门课总分的最高分,最低分及相应学生序号 将 5 门课总分按从大到小的顺序存入 ZCJ 中,相应学生序号存
入 XSXH
有 3 个多项式 P1(X)=X4+2X3+4X2+5,P2(X)=X+2,P3(X)=X2+2X+3,是进 行以下操作: 1. 求 P(X)=P1(X) +P2(X)P3(X) 2. 求 P(X)的根 3. 当 X 取矩阵 A 的每一行元素时,求 P(X)的值。其中 A= −1 1.2 −1.4 0.75 2 3.5 0 5 2.5
第5章 MATLAB绘图_习题答案
第5章 MATLAB绘图习题5一、选择题1.如果x、y均为4×3矩阵,则执行plot(x,y)命令后在图形窗口中绘制()条曲线。
DA.12 B.7 C.4 D.32.下列程序的运行结果是()。
Ax=0:pi/100:2*pi;for n=1:2:10plot(n*sin(x),n*cos(x))hold onendaxis squareA.5个同心圆B.5根平行线C.一根正弦曲线和一根余弦曲线D.5根正弦曲线和5根余弦曲线3.命令text(1,1,'{\alpha}+{\beta}')执行后,得到的标注效果是()。
CA.{\alpha}+{\beta} B.{\α}+{\β} C.α+βD.\α+\β4.subplot(2,2,3)是指()的子图。
AA.两行两列的左下图B.两行两列的右下图C.两行两列的左上图D.两行两列的右上图5.要使函数y=2e x的曲线绘制成直线,应采用的绘图函数是()。
CA.polar B.semilogx C.semilogy D.loglog6.下列程序的运行结果是()。
B[x,y]=meshgrid(1:5);surf(x,y,5*ones(size(x)));A.z=x+y平面B.与xy平面平行的平面C.与xy平面垂直的平面D.z=5x平面7.下列函数中不能用于隐函数绘图的是()。
DA.ezmesh B.ezsurf C.ezplot D.plot38.下列程序运行后,看到的图形()。
Ct=0:pi/20:2*pi;[x,y]=meshgrid(-8:0.5:8);z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);surf(x,y,z)view(0,90);axis equalA.像墨西哥帽子B.是空心的圆C.边界是正方形D.是实心的圆9.下列程序运行后得到的图形是()。
A[x,y]=meshgrid(-2:2);z=x+y;i=find(abs(x)<1 & abs(y)<1);z(i)=NaN;surf(x,y,z);shading interpA.在一个正方形的正中心挖掉了一个小的正方形B.在一个正方形的正中心挖掉了一个小的长方形C.在一个正方形的上端挖掉了一个小的正方形D.在一个正方形的下端挖掉了一个小的正方形10.在使用MA TLAB“绘图”选项卡中的命令按钮绘图之前,需要()。
matlab课后习题答案
第2章 MATLAB 矩阵运算基础2.1 在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375,并将其赋予变量a ? >> a=[5 7 3;4 9 1]2.5 计算矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535与⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242之和。
>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6]; >> a+b ans =7 7 7 9 14 13 15 12 142.6 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。
>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’ ans =4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i2.7 计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8]; >> a.*b ans =12 36 3 8 42 402.9 对于B AX =,如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ,⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ,求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;-0.5118 4.0427 1.33182.10 已知:⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a ,分别计算a 的数组平方和矩阵平方,并观察其结果。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第5章 MATLAB绘图
习题5
一、选择题
1.如果x、y均为4×3矩阵,则执行plot(x,y)命令后在图形窗口中绘制()条曲线。
D
A.12 B.7 C.4 D.3
2.下列程序的运行结果是()。
A
x=0:pi/100:2*pi;
for n=1:2:10
plot(n*sin(x),n*cos(x))
hold on
end
axis square
A.5个同心圆B.5根平行线
C.一根正弦曲线和一根余弦曲线D.5根正弦曲线和5根余弦曲线
3.命令text(1,1,'{\alpha}+{\beta}')执行后,得到的标注效果是()。
C
A.{\alpha}+{\beta} B.{\α}+{\β} C.α+βD.\α+\β
4.subplot(2,2,3)是指()的子图。
A
A.两行两列的左下图B.两行两列的右下图
C.两行两列的左上图D.两行两列的右上图
5.要使函数y=2e x的曲线绘制成直线,应采用的绘图函数是()。
C
A.polar B.semilogx C.semilogy D.loglog
6.下列程序的运行结果是()。
B
[x,y]=meshgrid(1:5);
surf(x,y,5*ones(size(x)));
A.z=x+y平面B.与xy平面平行的平面
C.与xy平面垂直的平面D.z=5x平面
7.下列函数中不能用于隐函数绘图的是()。
D
A.ezmesh B.ezsurf C.ezplot D.plot3
8.下列程序运行后,看到的图形()。
C
t=0:pi/20:2*pi;
[x,y]=meshgrid(-8:0.5:8);
z=sin(sqrt(x.^2+y.^2))./sqrt(x.^2+y.^2+eps);
surf(x,y,z)
view(0,90);axis equal
A.像墨西哥帽子B.是空心的圆
C.边界是正方形D.是实心的圆
9.下列程序运行后得到的图形是()。
A
[x,y]=meshgrid(-2:2);
z=x+y;
i=find(abs(x)<1 & abs(y)<1);
z(i)=NaN;
surf(x,y,z);shading interp
A.在一个正方形的正中心挖掉了一个小的正方形
B.在一个正方形的正中心挖掉了一个小的长方形
C.在一个正方形的上端挖掉了一个小的正方形
D.在一个正方形的下端挖掉了一个小的正方形
10.在使用MA TLAB“绘图”选项卡中的命令按钮绘图之前,需要()。
B
A.在命令行窗口中输入绘图命令B.在工作区中选择绘图变量
C.打开绘图窗口D.建立M文件
二、填空题
1.执行以下命令:
x=0:pi/20:pi;
y=sin(x);
以x为横坐标、y为纵坐标的曲线图绘制命令为,给该图形加上“正弦波”标题的命令为,给该图形的横坐标标注为“时间”,纵坐标标注为“幅度”的命令分别为和。
plot(x,y),title(‘正弦波’),xlabel(‘时间’),ylabel(‘幅度’) 2.在同一图形窗口中绘制y1和y2两条曲线,并对图形坐标轴进行控制,请补充程序。
x=-3:0.1:3;
y1=2*x+5;
y2=x.^2-3*x;
plot(x,y1) %绘制曲线y1
① ;
plot(x,y2) %绘制曲线y2
m1=max([y1,y2]);
m2=min([y1,y2]);
② ; %用axis调制坐标轴,横坐标在[-3,3]之间,纵坐标在[-5,20]之间
①hold on ②axis([-3,3,-5,20])
3.下列命令执行后得到的图形是。
(单位)圆
x=@(t) sin(t);
2
3
y=@(t) cos(t);
ezplot(x,y)
4.某工厂2015年度各季度产值分别为450.6、395.9、410.2、450.9,为了表示各季度产值占全年总产值的比例,可以绘制 ,其命令是 。
饼图,pie([450.6,395.9,410.2,450.9]) 或pie3([450.6,395.9,410.2,450.9])
三、应用题
1.绘制下列曲线。
(1)2
1100x y += (2)2221x e y -=π (3)122
=+y x (4)⎩⎨⎧==32
5t y t x (1)
x=-1:1/180:1;
y=100./(1+x.^2);
plot(x,y)
(2)
x=-1:1/180:1;
y=1/(2*pi).*exp(-x.^2./2);
plot(x,y)
(3)
ezplot('x^2+y^2-1=0')
(4)
ezplot('t^2','5*t^3',[-0.5,0.5])
2.分别用plot 和fplot 函数绘制函数x y 1sin
=的曲线,分析两曲线的差别。
3.绘制下列极坐标图。
(1)4cos 5+=θρ (2)θρ12
=
(3)7cos 5-=
θρ (4)23
θπρ= (1) theta=0:0.01:2*pi;
rho=5.*cos(theta)+4;
polar(theta,rho)
4
(2)
theta=0.1:0.1:2*pi;
rho=12./(sqrt(theta));
polar(theta,rho,'r')
(3)
theta=0:pi/90:2*pi;
t=cos(theta);
a=find(t<0.1);
t(a)=0.01;
rho=5./t-7;
polar(t,rho)
(4)
theta=0:pi/90:2*pi;
rho=pi/3.*(theta.^2);
polar(theta,rho)
4.在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。
(1)y=2x-0.5
(2)π≤≤⎩⎨⎧==t 0,sin )3sin(cos )3sin(t t y t
t x
t=0:0.01:2*pi;
x=sin(3*t).*cos(t);
y=sin(3*t).*sin(t);
plot(x,y);
hold on;
x=-1:0.01:1;
y=2*x-0.5;
plot(x,y);
hold off
5.绘制下列三维图形。
(1)⎪⎩⎪⎨⎧===t z t y t x sin cos (2)⎪⎩
⎪⎨⎧=+=+=u z v u y v u x sin sin )cos 1(cos )cos 1(
(3)5=z (4)半径为10的球面
(1)
t=0:pi/180:2*pi;
x=cos(t);
y=sin(t);
z=t;
plot3(x,y,z)
(2)
u=0:pi/50:2*pi;
v=0:pi/50:2*pi;
ezsurf('(1+cos(u))*cos(v)','(1+cos(u))*sin(v)','sin(u)') (4) 为什么用不了ezplot3('x^2+y^2+z^2-100')???????????u=0:pi/50:2*pi;
v=0:pi/50:2*pi;
ezsurf('cos(u)*cos(v)','cos(u)*sin(v)','sin(u)')
5。