1. 设 y 3 e 4tsin 43t
2
3 要求以0.01秒为间隔,求出y的151个点,并求出其导数的值和曲线。 clc clear
x=0:0.01:1.5;
y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)
y1=diff(y) subplot(2,1,1) plot(x,y)
subplot(2,1,2) plot(x(1:150),y1)
2绘制极坐标系下曲线(a,b,n自定数据)
acos b n clc clear a=10; b=pi/2; n=5;
theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)
3. 列出求下列空间曲面交线的程序
z1 x2 2y2
z2 2x 3yclc clear
x=[-5:0.5:5];
[X,Y]=meshgrid(x);
z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x')ylabel('y') zlabel('z') surf(X,Y,z1) hold on
surf(X,Y,z2)
k=find(abs(z1-z2)<0.5); x1=X(k)y1=Y(k)
z3=x1.^2-2*y1.^2 hold on
plot3(x1,y1,z3,'*')
4、设 y cosx 0.5 3sinx
(1 x2) 把x=0~2π间分为101点,画出以x为横坐标,y
为纵坐标的曲线,要求有图形标注。
x=-2*pi:0.1: 2*pi;
y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x坐标'); ylabel('y坐标'); legend('原函数')
gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')
5、求下列联立方程的解 3x 4y 7z 12w 4
5x 7y 4z 2w 3
x 8z 5w 9
6x 5y 2z 10w 8clc clear
a=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];
b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)
6. 假设一曲线数据点为
x = 0:2:4*pi;
y = sin(x).*exp(-x/5);
试将x 的间距调成 0.1,采用不同插值方法进行插值,并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。 clc clear
x=0:2:4*pi;
y=sin(x).*exp(-x./5); x1=0:0.1:4*pi;
y01=interp1(x,y,x1,'spline'); y02=interp1(x,y,x1,'cubic'); y03=interp1(x,y,x1,'nearest' ); y04=interp1(x,y,x1);subplot(3,2,1); plot(x,y,'o-'); title('Ô-º¯Êý'); subplot(3,2,2); plot(x1,y01,'gx-'); title('spline'); subplot(3,2,3);plot(x1,y02,'m+-'); title('cubic'); subplot(3,2,4); plot(x1,y03,'r*-'); title('nearest'); subplot(3,2,5);plot(x1,y04,'kd-'); title('line');