p110
5>> subplot(2,2,1),ezplot('x.^2.*cos(1./x)')
> subplot(2,2,2),ezplot('cos(1./x)')
> x1=-0.01:.001:0.01;
> subplot(2,2,3),plot(x1,x1.^2.*cos(1./x1))
> subplot(2,2,4),plot(x1,cos(1./x1))
由图像可知,函数h(x)在x=0处极限是0,而函数k(x)在x=0处的极限不存在。
p110.14(a)
x=linspace(-10,5,100);
subplot(1,2,1),plot(x,x.^3+x.^2-x),grid on,hold on%画曲线。
在一个点上,取不同的步长求斜率,画出多根割线。
h=[1,0.1,0.01];
for k=1:3
x1=[1,1+h(k)];f=x1.^3+x1.^2-x1;
q=diff(f)/diff(x1),plot(x,f(1)+q*(x-x1(1)),end
q =9 4.4100 4.0401
在多个点上,去同样步长,求斜率,画出多根切线。
for k=0:2
x1=[k,k+0.1];f=x1.^3+x1.^2-x1;
q=diff(f)/diff(x1),subplot(1,2,2),plot(x,x.^3+x.^2-x),grid on,hold on
plot(x,f(1)+q*(x-x1(1)),r'),end
q =-0.8900 4.4100 15.7100
p110.14.(b)
x=linspace(-5,5,100);
subplot(1,2,2),plot(x,x.^(1/3)+x.^(2/3)),grid on,hold on %画曲线。
在一个点上,取不同的步长求斜率,画出多根割线。
h=[1,0.1,0.01];
for k=1:3
x1=[1,1+h(k)];f=x1.^(1/3)+x1.^(2/3);
q=diff(f)/diff(x1),plot(x,f(1)+q*(x-x1(1)),endq =
在多个点上,去同样步长,求斜率,画出多根切线。
for k=0:2
x1=[k,k+0.1];f=,x1.^(1/3)+x1.^(2/3;
q=diff(f)/diff(x1),subplot(1,2,2),plot(x,x.^(1/3)+x.^(2/3)),grid on,hold on
plot(x,f(1)+q*(x-x1(1)),r'),end
q =6.7960 0.9788 0.7314
p110.14.(c)
x=linspace(-6,6,100);
plot(x,sin(2*x)),grid on ,hold on %画曲线。
在一个点上,取不同的步长求斜率,画出多根割线。
h=[0.1,0.01];
for k=1:3
x1=[1.1+h(k)];f=sin(2*x1);
q=diff(f)/diff(x1),plot(x,f(1)+q*(x-x1(1)),end
q= -1.0080 -0.8504
在多个点上,去同样步长,求斜率,画出多根切线。
for k=0:2
x1=[k,k+0.1];f=sin(2*x1);
q=diff(f)/diff(x1),subplot(1,2,2),plot(x,sin(2*x)),grid on,hold on
plot(x,f(1)+q*(x-x1(1)),r'),endq =
q =q =
p110.14.(d)
x=linspace(-5,5,100);
plot(x,x.^2.*cos(x)),grid on ,hold on%画曲线。
在一个点上,取不同的步长求斜率,画出多根割线。
h=[0.1,0.01];
x1=[1.1+h];f=x1.^2.*cos(x1);
q=diff(f)/diff(x1),plot(x,f(1)+q*(x-x1(1)),q =
在多个点上,去同样步长,求斜率,画出多根切线。
for k=0:2
x1=[k,k+0.1];f=x1.^2.*cos(x1);
q=diff(f)/diff(x1),subplot(1,2,2),plot(x,x.^2.*cos(x)),grid on,hold on
plot(x,f(1)+q*(x-x1(1)),r'),endq =
p123 .1(a)
a=[1,2,3,4];
theta=linspace(0,2*pi,250);
for i=1: 4
rho(i,:)a(i)*theta;
subplot(2,2,i),polar(theta,rho(i,:)
endgtext('rho=theta');
gtext('rho=2*theta');
gtext('rho=3*theta');
gtext('rho=4*theta');
p123 .1(b)
a=[5,6,7,8];
theta=linspace(0,2*pi,250);
for i=1: 4
rho(i,:)exp(a(i)*theta);
subplot(2,2,i),polar(theta,rho(i,:)
endgtext('rho=exp(5*theta)')
gtext('rho=exp(6*theta)')
gtext('rho=exp(7*theta)')
gtext('rho=exp(8*theta)')
p123 .1(c)
a=[1,3,5,7];
theta=linspace(0,2*pi,250);
for i=1: 4
rho(i,:)a(i)*(1-cos(theta));
subplot(2,2,i),polar(theta,rho(i,:)
endgtext('rho=1*(1-cos(theta))'
gtext('rho=3*(1-cos(theta))'
gtext('rho=5*(1-cos(theta))'
gtext('rho=7*(1-cos(theta))'
p123 .1(d)
a=[2,4,6,8];
theta=linspace(0,2*pi,250);
for i=1: 4
rho(i,:)a(i)*cos(theta);
subplot(2,2,i),polar(theta,rho(i,:)
endgtext('rho=2*cos(theta)')
gtext('rho=4*cos(theta)')
gtext('rho=6*cos(theta)')
gtext('rho=8*cos(theta)')
p123.3.(a)
t=0:.1:5*pi;
x=t.*cos(t);
y=2*t.*sin(t);
z=1.5*t;
plot3(x,y,z,'*y')
p123.3(b)
[x,y]=meshgrid(-2:.1:2); 确定计算和绘图的定义域网络。
z1=x.^2+y.^2第一个曲线方程。
z2=sqrt(1-x.^2-y.^2); 第二个曲线方程。
mesh(x,y,z1);hold on;mesh(x,y,z2); 分别画出两个曲面。
r0=abs(z1-z2)<=1; %求两曲面z坐标之差小于0.1的网络矩阵。
zz=r0.*z1;yy=r0.*y;xx=r0.*x; %求这些网络上的坐标值。
plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0画出这些点。
colormap(gray),hold off
p124.5(a)
> [x,y]=meshgrid(-pi:.1:pi);
> z=exp(-(x.^2+y.^2)./4).*cos(x).^2+sin(y).^2);
> mesh(x,y,z)
p124.5.(b)
u,v]=meshgrid(0:.1:2*pi);
x=(3+cos(u)).cos(v);
y=(3+cos(u)).sin(v);
z=sin(v);
shading flat;
view(20,0);
subplot(2,2,1),mesh(x,y,z);
subplot(2,2,2),surfl(x,y,z);
subplot(2,2,3),clorbar
subplot(2,2,3),colorbar;
subplot(2,2,4),contour3(x,y,z)
x,y]=meshgrid(-2:.1:2确定计算和绘图的定义域网络。
z=x.^3-3*x.*y.^2+y.^2;
mesh(x,y,z画出曲面。
contour(x,y,z,20在xy平面上画出20条等高线。
hold on
px,py]=gradient(z,2,2);
quiver(x,y,px,py画出梯度向量。
x,y]=meshgrid(-2/3:.1:2/3); 确定计算和绘图的定义域网络。
z=2*x.^4+y.^4-2*x.^2-2*y.^2+3;
matlab作业
2011029170002王柳。a 一个问题的病态性如何,与求解它的算法有关系。错 b 无论问题是否病态,好的算法都会得到它好的近似解。错 c 计算中使用更高的精度,可以改善问题的病态性。错 d 用一个稳定的算法计算一个良态问题,一定会得到它好的近似解。对 e 浮点数在整个数轴上是均匀分布。错 f ...
matlab作业
matlab语言 第3次作业 字符串,单元数组和结构体 专业 海洋技术 海洋测绘方向 姓名 张体强学号 1026222 1 如何将一个char 数据类型的向量转化为相应的double 型数据类型的数据向量。从式1 到8,判断这些语句是否正确。如果它们正确,那么将产生什么结果?这题不要在电脑中做。1....
matlab作业
电子与通信工程学院。通信系统 实验报告。2013 2014 学年第1学期。调频 fm 系统调制解调 专业 通信工程。班级 通信111 班。学号 姓名 指导教师姓名陈多瑜。2013年 11 月日。1.频率调制或调频 fm 1 设调制信号为m t 调频信号的数学表达式为。例如 m t 的时域波形为。m ...