jy022706523
2009-07-08, 12:12
RT，新学MATLAB，副下列作业列表。。。有几个都不会做，高手们帮帮忙啊~~:D
1、矩阵的相关运算
设有两个矩阵：A=[11.59 12.81 15.66; 15.2 4.18 13.61; 10.59 7.59 9.22]
B=[16.00 4.41 –10.37 –21.61; 0.88 –20.04 12.86 8.56;
−1.43 10.71 18.81 –5.99; –12.48 24.35 –23.9 10.34]
分别求出这两个矩阵的逆矩阵、行列式的值、特征值，特征向量，LU和QR分解。
2、数值运算和符号运算
1）求解下列微分方程：
2）已知两个多项式f(x)=x3-3x2-3x+1，g(x)=x3-7x+6，试进行以下操作：
1）用数值的方法进行多项式的算术运算、求根、求值和求导运算；
2）用符号运算的方法进行多项式的算术运算；并分别进行因式分解。
3、二维、三维图形的绘制
1）绘出如下图所示的图形并写出绘制方法。
图片地址:[http://user.qzone.qq.com/312936686/photo/c0d6698a-0c3e-43aa-b3c1-82214b687c50/Mzh3LL37TUJAt1KhV2bRjozfuSwK0ZwAAA!!/]
4、编制MATLAB程序实现特定功能
编写程序计算y(x)([-3, 3]，步长0.01)，并画出其在[-3, 3]上的曲线。表达式如下：
第1题：
A=[11.59 12.81 15.66; 15.2 4.18 13.61; 10.59 7.59 9.22]；
B=[16.00 4.41 -10.37 -21.61;0.88 -20.04 12.86 8.56;-1.43 10.71 18.81 -5.99;-12.48 24.35 -23.9 10.34]；
inv(A)
ans =
-0.1564 0.0018 0.2630
0.0096 -0.1425 0.1940
0.1718 0.1152 -0.3533
inv(B)
ans =
0.2339 0.4104 0.1435 0.2323
0.0784 0.1326 0.0804 0.1006
0.0162 0.0530 0.0497 0.0188
0.1352 0.3055 0.0988 0.1835
det(A)
ans =
413.9417
det(B)
ans =
-4.4630e+004
[V,D]=eig(A)
V =
-0.6727 -0.5684 -0.5345
-0.5694 0.8226 -0.3679
-0.4726 -0.0142 0.7609
D =
33.4357 0 0
0 -6.5578 0
0 0 -1.8879
[V,D]=eig(B)
V =
0.3770 -0.7873 0.7363 -0.8112
-0.6637 -0.2878 -0.0532 0.1579
0.2194 -0.0752 0.1674 0.5614
0.6077 -0.5401 -0.6535 0.0421
D =
-32.6289 0 0 0
0 1.7953 0 0
0 0 32.5039 0
0 0 0 23.4397
[L,U]=lu(A)
L =
0.7625 1.0000 0
1.0000 0 0
0.6967 0.4861 1.0000
U =
15.2000 4.1800 13.6100
0 9.6227 5.2824
0 0 -2.8301
[L,U]=lu(B)
L =
1.0000 0 0 0
0.0550 -0.7299 -0.3234 1.0000
-0.0894 0.3996 1.0000 0
-0.7800 1.0000 0 0
U =
16.0000 4.4100 -10.3700 -21.6100
0 27.7898 -31.9886 -6.5158
0 0 30.6651 -5.3178
0 0 0 3.2732
[Q,R]=qr(A)
Q =
-0.5304 0.7367 -0.4195
-0.6956 -0.6610 -0.2814
-0.4846 0.1426 0.8630
R =
-21.8521 -13.3800 -22.2409
0 7.7558 3.8540
0 0 -2.4424
[Q,R]=qr(B)
Q =
-0.7858 0.4746 0.1986 0.3433
-0.0432 -0.6293 -0.0082 0.7759
0.0702 0.3164 -0.9121 0.2509
0.6129 0.5278 0.3585 0.4660
R =
-20.3610 13.0779 -5.7350 22.5286
0 30.9458 -19.6781 -12.0806
0 0 -27.8895 4.8097
0 0 0 2.5397
第2题：
（1）：
（2）：
>> p1=[1 -3 -3 1]
p1 =
1 -3 -3 1
>> roots(p1)
ans =
3.7321
-1.0000
0.2679
poly(p1)
ans =
1 4 -2 -12 9
polyder(p1)
ans =
3 -6 -3
p2=[3 0 -7 6]
roots(p2)
ans =
-1.8481
0.9241 + 0.4778i
0.9241 - 0.4778i
poly(p2)
ans =
1 -2 -45 126 0
polyder(p2)
ans =
9 0 -7
第3题：
（1）：
x=0:0.2:16;
y1=sin(x);
y2=cos(x);
plot(x,y1,'g--',x,y2,'b-')
(2)
第4题：
i=0;
for x=-3:0.01:3
i=i+1;
X(i)=x;
if and(x>=-3,x<-1)
Y(i)=(-x*x-4*x-3)/2;
elseif and(x>=-1,x<1)
Y(i)=-x*x+1;
else
Y(i)=(-x*x+4*x-3)/2;
end
end
plot(X,Y)
1、矩阵的相关运算
设有两个矩阵：A=[11.59 12.81 15.66; 15.2 4.18 13.61; 10.59 7.59 9.22]
B=[16.00 4.41 –10.37 –21.61; 0.88 –20.04 12.86 8.56;
−1.43 10.71 18.81 –5.99; –12.48 24.35 –23.9 10.34]
分别求出这两个矩阵的逆矩阵、行列式的值、特征值，特征向量，LU和QR分解。
2、数值运算和符号运算
1）求解下列微分方程：
2）已知两个多项式f(x)=x3-3x2-3x+1，g(x)=x3-7x+6，试进行以下操作：
1）用数值的方法进行多项式的算术运算、求根、求值和求导运算；
2）用符号运算的方法进行多项式的算术运算；并分别进行因式分解。
3、二维、三维图形的绘制
1）绘出如下图所示的图形并写出绘制方法。
图片地址:[http://user.qzone.qq.com/312936686/photo/c0d6698a-0c3e-43aa-b3c1-82214b687c50/Mzh3LL37TUJAt1KhV2bRjozfuSwK0ZwAAA!!/]
4、编制MATLAB程序实现特定功能
编写程序计算y(x)([-3, 3]，步长0.01)，并画出其在[-3, 3]上的曲线。表达式如下：
第1题：
A=[11.59 12.81 15.66; 15.2 4.18 13.61; 10.59 7.59 9.22]；
B=[16.00 4.41 -10.37 -21.61;0.88 -20.04 12.86 8.56;-1.43 10.71 18.81 -5.99;-12.48 24.35 -23.9 10.34]；
inv(A)
ans =
-0.1564 0.0018 0.2630
0.0096 -0.1425 0.1940
0.1718 0.1152 -0.3533
inv(B)
ans =
0.2339 0.4104 0.1435 0.2323
0.0784 0.1326 0.0804 0.1006
0.0162 0.0530 0.0497 0.0188
0.1352 0.3055 0.0988 0.1835
det(A)
ans =
413.9417
det(B)
ans =
-4.4630e+004
[V,D]=eig(A)
V =
-0.6727 -0.5684 -0.5345
-0.5694 0.8226 -0.3679
-0.4726 -0.0142 0.7609
D =
33.4357 0 0
0 -6.5578 0
0 0 -1.8879
[V,D]=eig(B)
V =
0.3770 -0.7873 0.7363 -0.8112
-0.6637 -0.2878 -0.0532 0.1579
0.2194 -0.0752 0.1674 0.5614
0.6077 -0.5401 -0.6535 0.0421
D =
-32.6289 0 0 0
0 1.7953 0 0
0 0 32.5039 0
0 0 0 23.4397
[L,U]=lu(A)
L =
0.7625 1.0000 0
1.0000 0 0
0.6967 0.4861 1.0000
U =
15.2000 4.1800 13.6100
0 9.6227 5.2824
0 0 -2.8301
[L,U]=lu(B)
L =
1.0000 0 0 0
0.0550 -0.7299 -0.3234 1.0000
-0.0894 0.3996 1.0000 0
-0.7800 1.0000 0 0
U =
16.0000 4.4100 -10.3700 -21.6100
0 27.7898 -31.9886 -6.5158
0 0 30.6651 -5.3178
0 0 0 3.2732
[Q,R]=qr(A)
Q =
-0.5304 0.7367 -0.4195
-0.6956 -0.6610 -0.2814
-0.4846 0.1426 0.8630
R =
-21.8521 -13.3800 -22.2409
0 7.7558 3.8540
0 0 -2.4424
[Q,R]=qr(B)
Q =
-0.7858 0.4746 0.1986 0.3433
-0.0432 -0.6293 -0.0082 0.7759
0.0702 0.3164 -0.9121 0.2509
0.6129 0.5278 0.3585 0.4660
R =
-20.3610 13.0779 -5.7350 22.5286
0 30.9458 -19.6781 -12.0806
0 0 -27.8895 4.8097
0 0 0 2.5397
第2题：
（1）：
（2）：
>> p1=[1 -3 -3 1]
p1 =
1 -3 -3 1
>> roots(p1)
ans =
3.7321
-1.0000
0.2679
poly(p1)
ans =
1 4 -2 -12 9
polyder(p1)
ans =
3 -6 -3
p2=[3 0 -7 6]
roots(p2)
ans =
-1.8481
0.9241 + 0.4778i
0.9241 - 0.4778i
poly(p2)
ans =
1 -2 -45 126 0
polyder(p2)
ans =
9 0 -7
第3题：
（1）：
x=0:0.2:16;
y1=sin(x);
y2=cos(x);
plot(x,y1,'g--',x,y2,'b-')
(2)
第4题：
i=0;
for x=-3:0.01:3
i=i+1;
X(i)=x;
if and(x>=-3,x<-1)
Y(i)=(-x*x-4*x-3)/2;
elseif and(x>=-1,x<1)
Y(i)=-x*x+1;
else
Y(i)=(-x*x+4*x-3)/2;
end
end
plot(X,Y)