matlab第2章作业
matlab 部分答案
课
D=A.*B
后
答
MATLAB: A=[1,4,8,13; -3,6,-5,-9; 2,-7,-12,-8];
案
网
8 13 ⎤ ⎡ 5 4 3 −2 ⎤ ⎡1 4 ⎥ ⎢ ⎥ 10. 设 A = ⎢ −3 6 −5 −9 ⎥ ， B = ⎢ ⎢ 6 −2 3 − 8 ⎥ ⎢ ⎢ ⎣ −1 3 −9 7 ⎥ ⎦ ⎣ 2 −7 −12 −8⎥ ⎦
2
+ y2 )
，求定义域 x = [−2, 2] ， y = [−2, 2] 内的 z 值（网格取 0.1 见方） ，并绘制
15．设 x = cos(t ) ， y = sin( Nt + α ) ，若 N = 2 ， α = 0, π / 3, π / 2, π ，在 4 个子图中分别画出其
MATLAB: t=0:0.01:3; N=2;
4. 用三点公式计算 y=f(x)在 x=1.0,1.2 处的导数值， f(x)值由下表给出； x f(x)
5
1.0 0.25
4
1.1
ww w
1.2
2
.k
1.3 0.1890 1.4 0.1736
后
7. 求解下列线性常微分方程的解析解。
答
6．设方程的根为 x = [−3, −5, −8, −9] ，求它们对应的 x 多项式的系数。
π
2π
课
=∫
0
17．设 ( X , Y ) 的概率密度为
网
⎧ d 2 x dy + − x = et ⎪ ⎪ dt 2 dt ⎨ 2 ⎪ d y + dx + y = 0 ⎪ ⎩ dt 2 dt
dy 。 dx
2π
matlab课后习题解答第二章
第2章符号运算习题 2 及解答1说出以下四条指令产生的结果各属于哪一种数据种类，是“双精度”对象，仍是“符号”符号对象3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+)〖目的〗不可以从显示形式判断数据种类，而一定依赖class 指令。
〖解答〗c1=3/7+c2=sym(3/7+c3=sym('3/7+')c4=vpa(sym(3/7+)Cs1=class(c1)Cs2=class(c2)Cs3=class(c3)Cs4=class(c4)c1 =c2 =37/70c3 =c4 =Cs1 =doubleCs2 =symCs3 =symCs4 =sym2在不加特意指定的状况下，以下符号表达式中的哪一个变量被以为是自由符号变量 .sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')〖目的〗理解自由符号变量确实认规则。
〖解答〗symvar(sym('sin(w*t)'),1)ans =wsymvar(sym('a*exp(-X)'),1)ans = asymvar(sym('z*exp(j*th)'),1) ans = za11a12a 135 求符号矩阵 A a 21 a 22 a 23 的队列式值和逆， 所得结果应采纳 “子a31a32a33表达式置换”简短化。
〖目的〗理解 subexpr 指令。
〖解答〗A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A);[IAs,d]=subexpr(IA,d)A =[ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs =[ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)8（1）经过符号计算求 y(t)sin t 的导数dy。
matlab1-8章课后作业
MATLAB基础教程 1~8章作业Matlab第一章1.阐述Matlab的功能Matlab作为一种高级计算软件，是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境，已被广泛应用于不同领域。
Matlab的基本功能包括：数学计算功能、图形化显示功能、M语言编程功能、编译功能、图形用户界面开发功能、Simulink 建模仿真功能、自动代码生成功能。
Matlab第二章1．创建double的变量，并进行计算。
(1) a=87，b=190，计算 a+b、a-b、a*b。
(2)创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。
>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =>> c*dans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2，v=3，计算：(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u)ans =15.4189>> sqrt(u-3*v)/(u*v)ans =0 + 0.4410i4.计算如下表达式：(1)(2)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i>> sin(2-8*i)ans =1.3553e+003 +6.2026e+002i5.判断下面语句的运算结果。
matlab第二章习题答案
第一大题：（1）a = 7/3b = sym(7/3)c = sym(7/3,'d')d = sym('7/3')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =2.3333b =7/3c =2.3333333333333334813630699500209d =7/3v1 =0.0v2 =0.0v3 =0.00000000000000014802973661668756666666667788716（2）a = pi/3b = sym(pi/3)c = sym(pi/3,'d')d = sym('pi/3')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =1.0472b =pi/3c =1.047197551196597631317786181171d =pi/3v1 =0.0v2 =0.0v3 =0.00000000000000011483642827992216762806615818554(3)a = pi*3^(1/3)b = sym(pi*3^(1/3))c = sym(pi*3^(1/3),'d')d = sym('pi*3^(1/3)')v1=vpa(abs(a-d))v2=vpa(abs(b-d))v3=vpa(abs(c-d))a =4.5310b =1275352044764433/281474976710656c =4.5309606547207899041040946030989d =pi*3^(1/3)v1 =0.00000000000000026601114166290944374842393221638 v2 =0.00000000000000026601114166290944374842393221638 v3 =0.0000000000000002660111416629094726767991785515第二大题：（1）c1=3/7+0.1c1 =0.5286双精度（2）c2=sym(3/7+0.1)c2 =37/70符号（3）c3=vpa(sym(3/7+0.1))c3 =0.52857142857142857142857142857143完整显示精度第三大题：（1）findsym(sym('sin(w*t)'),1)ans =w(2)findsym(sym('a*exp(-X)' ) ,1)ans =a(3)findsym(sym('z*exp(j*theta)'),1)ans =z第四大题：A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]')A =[ a11, a12, a13][ a21, a22, a23][ a31, a32, a33]DA=det(A)DA =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31w=inv(A)w =[ (a22*a33 - a23*a32)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), -(a12*a33 -a13*a32)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), (a12*a23 - a13*a22)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] [ -(a21*a33 - a23*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), (a11*a33 -a13*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), -(a11*a23 - a13*a21)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] [ (a21*a32 - a22*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 +a12*a23*a31 + a13*a21*a32 - a13*a22*a31), -(a11*a32 -a12*a31)/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 +a13*a21*a32 - a13*a22*a31), (a11*a22 - a12*a21)/(a11*a22*a33 -a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)] IAs=subexpr(w,'d')d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)IAs =[ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 -a13*a22)][ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 -a13*a21)][ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 -a12*a21)]第六大题：syms ksyms x positives_s=2/(2*k+1)*((x-1)/(x+1))^(2*k+1)s_ss=simple(symsum(s_s,k,0,inf))s_s =(2*((x - 1)/(x + 1))^(2*k + 1))/(2*k + 1)警告: simple will be removed in a future release. Use simplify instead. [> In sym.simple at 41]s_ss =log(x)第八大题：syms x clearsyms xh=exp(-abs(x))*abs(sin(x))si=vpa(int(h,-5*pi,1.7*pi),64)h =abs(sin(x))*exp(-abs(x))si =1.087849417255503701102633764498941389696991336803454392428439159 第九大题：syms x y clearsyms x yr=int(int(x^2+y^2,y,1,x^2),x,1,2)r =1006/105第十大题：syms t x;f=sin(t)/t;y=int(f,t,0,x)y1=subs(y,x,sym('4.5'))ezplot(y,[0,2*pi])y =sinint(x)y1 =syms x clearsyms x ny=sin(x)^nyn=int(y,0,1/2*pi)y31=vpa(subs(yn,n,sym('1/3')))y32=vpa(subs(yn,n,1/3))y =sin(x)^nyn =piecewise([-1 < real(n), beta(1/2, n/2 + 1/2)/2], [real(n) <= -1, int(x^n/(1 - x^2)^(1/2), x, 0, 1)])y31 =1.2935547796148952674767575125656y32 =1.2935547796148952674767575125656第二十题:clearsyms y xy=dsolve('(Dy*y)/5+x/4=0','x')y =2^(1/2)*(C6 - (5*x^2)/8)^(1/2)-2^(1/2)*(C6 - (5*x^2)/8)^(1/2)y1=subs(y,'C6',1)y1 =2^(1/2)*(1 - (5*x^2)/8)^(1/2)-2^(1/2)*(1 - (5*x^2)/8)^(1/2)clfhy1=ezplot(y1(1),[-2,2,-2,2],1)set(hy1,'Color','r')grid onhold onhy2=ezplot(y1(2),[-2,2,-2,2],1)set(hy2,'Color','b')grid onxlabel('Y')ylabel('X')hold offbox onlegend('y(1)','y(2)','Location','Best')hy1 =174.0155hy2 =177.0145。
matelab作业2参考答案
matelab作业2参考答案Matlab作业2参考答案Matlab作业2是一项综合性的任务，要求学生运用Matlab编程语言解决一系列数学问题。
本文将为大家提供一份参考答案，帮助学生更好地理解和完成这项作业。
首先，我们将讨论作业的第一个问题，即给定一个矩阵A，求解其特征值和特征向量。
在Matlab中，可以使用eig函数来实现这一功能。
例如，假设我们有一个3×3的矩阵A，可以按照以下方式计算其特征值和特征向量：```A = [1 2 3; 4 5 6; 7 8 9];[eigenvectors, eigenvalues] = eig(A);```在上述代码中，变量eigenvectors将存储A的特征向量，而变量eigenvalues 将存储A的特征值。
通过打印这两个变量的值，我们可以得到矩阵A的特征值和特征向量。
接下来，我们将探讨作业的第二个问题，即求解线性方程组。
假设我们有一个3×3的系数矩阵A和一个3×1的常数向量b，我们需要求解方程组Ax=b。
在Matlab中，可以使用backslash运算符来求解线性方程组。
例如，假设我们有以下方程组：```A = [1 2 3; 4 5 6; 7 8 9];b = [10; 20; 30];x = A \ b;```在上述代码中，变量x将存储方程组的解。
通过打印变量x的值，我们可以得到方程组的解。
此外，作业的第三个问题要求学生使用Matlab绘制函数图像。
在Matlab中，可以使用plot函数来实现这一功能。
例如，假设我们要绘制函数y=sin(x)，其中x的取值范围为0到2π，可以按照以下方式绘制函数图像：```x = 0:0.1:2*pi;y = sin(x);plot(x, y);```在上述代码中，变量x将存储x的取值范围，变量y将存储对应的函数值。
通过调用plot函数，我们可以将函数y=sin(x)的图像绘制出来。
第2章--MATLAB数据及其运算-习题答案教学内容
第2章--M A T L A B数据及其运算-习题答案第2章 MATLAB数据及其运算习题2一、选择题1．下列可作为MATLAB合法变量名的是（）。
D A．合计 B．123 C．@h D．xyz_2a 2．下列数值数据表示中错误的是（）。
CA．+10 B．1.2e-5 C．2e D．2i3．使用语句t=0:7生成的是（）个元素的向量。
A A．8 B．7 C．6 D．54．执行语句A=[1,2,3;4,5,6]后，A(3)的值是（）。
B A．1 B．2 C．3 D．45．已知a为3×3矩阵，则a(:,end)是指（）。
D A．所有元素 B．第一行元素C．第三行元素 D．第三列元素6．已知a为3×3矩阵，则运行a (1)=[]后（）。
A A．a变成行向量 B．a变为2行2列C．a变为3行2列 D．a变为2行3列7．在命令行窗口输入下列命令后，x的值是（）。
B >> clear>> x=i*jA．不确定 B．-1 C．1 D．i*j 8．fix(354/100)+mod(354,10)*10的值是（）。
D A．34 B．354 C．453 D．439．下列语句中错误的是（）。
BA．x==y==3 B．x=y=3C．x=y==3 D．y=3,x=y10．find(1:2:20>15)的结果是（）。
CA．19 20 B．17 19C．9 10 D．8 911．输入字符串时，要用（）将字符括起来。
C A．[ ] B．{ } C．' ' D．" " 12．已知s='显示"hello"'，则s的元素个数是（）。
A A．9 B．11 C．7 D．1813．eval('sqrt(4)+2')的值是（）。
BA．sqrt(4)+2 B．4 C．2 D．2， 214．有3×4的结构矩阵student，每个结构有name（姓名）、scores（分数）两个成员，其中scores是以1×5矩阵表示的5门课的成绩，那么要删除第4个学生的第2门课成绩，应采用的正确命令是（）。
数字图像处理及应用（MATLAB）第2章习题答案
7.平均值说明f (x ，y )的平均值等于其傅里叶变换F (u ，v )在频率原点的值F (0,0)。
2-3证明离散傅里叶变换的频率位移和空间位移性质。
证明：)(2101),(1),(NvyM ux j M x N y e y x f MN v u F +--=-=∑∑=π),(),(1),(100)(21010)(2)(21010000v v u u F dxdy ey x f MNe ey x f MN y Nv v x M u u j M x N y N yv M x u j Nvy M ux j M x N y --==-+---=-=++--=-=∑∑∑∑πππ因为()()v u F y x f ,,⇔ 所以 ),(),(00)(200v v u u F e y x f N y v M x u j --⇔+π2-4小波变换是如何定义的？小波分析的主要优点是什么？小波之所以小，是因为它有衰减性，即是局部非零的；而称为波，则是因为它有波动性，即其取值呈正负相间的振荡形式，将)(2R L 空间的任意函数f (t )在小波基下展开，称其为函数f (t )的连续小波变换。
小波变换是时间（空间）频率的局部化分析，它通过伸缩平移运算对信号（函数）逐步进行多尺度细化，最终达到高频处时间细分，低频处频率细分，能自动适应时频信号的要求从而可聚焦到信号的任意细节，解决了Fourier 变换的困难问题。
2-5 在图像缩放中，采用最近邻域法进行放大时，如果放大倍数太大，可能会出现马赛克效应，这个问题有没有办法解决，或者有所改善。
可以利用线性插值法，当求出的分数地址与像素点不一致时，求出周围四个像素点的距离比，根据该比率， 由四个邻域的像素灰度值进行线性插值。
2-6 复合变换的矩阵等于基本变换的矩阵按顺序依次相乘得到的组合矩阵。
即，T=T N T N-1…T 1。
问矩阵顺序的改变能否影响变换的结果。
矩阵顺序的改变不会影响变换的结果。
MATLAB作业2参考答案
MATLAB作业二参考答案1、试求出如下极限。
（1）2325(2)(3)lim(5)x xxxx xx+++→∞+++，（2）23312lim()xyx y xyx y→-→++，（3）2222221cos()lim()x yxyx yx y e+→→-++【求解】极限问题可以由下面语句直接求解。
>> syms x; f=(x+2)^(x+2)*(x+3)^(x+3)/(x+5)^(2*x+5); limit(f,x,inf)ans =exp(-5)>> syms x yfa=(x^2*y+x*y^3)/(x+y)^3; limit(limit(fa,x,-1),y,2) ans =-6>> fc=(1-cos(x^2+y^2))*exp(x^2+y^2)/(x^2+y^2);limit(limit(fc,x,0),y,0)ans =2、试求出下面函数的导数。
（1）()y x=22atan ln()yx yx=+【求解】由求导函数diff() 可以直接得出如下结果，其中(2) 为隐函数，故需要用隐函数求导公式得出导数。
>> syms x;f=sqrt(x*sin(x)*sqrt(1-exp(x))); simple(diff(f))ans =1/2/(x*sin(x)*(1-exp(x))^(1/2))^(1/2)*(sin(x)*(1-exp(x))^(1/2)+x*cos(x)*(1-exp(x))^(1/2)-1/2*x*sin(x)/(1-exp(x))^(1/2)*exp(x))>> syms x,y; f=atan(y/x)-log(x^2+y^2);f1=simple(-diff(f,x)/diff(f,y))f1 =(y+2*x)/(x-2*y)3、假设1cosu-=，试验证22u ux y y x∂∂=∂∂∂∂。
【求解】证明二者相等亦可以由二者之差为零来证明，故由下面的语句直接证明。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
深圳大学实验报告课程名称：可视化计算语言
实验项目名称：Matlab矩阵及其运算
学院：医学院
专业：生物医学工程
指导教师：
报告人：学号：班级：
实验时间：2012-9-26
实验报告提交时间：
教务部制
实验目的与要求：
1．掌握MATLAB 数据对象的特点以及数据的运算规则。
2．掌握MATLAB 中建立矩阵的方法以及矩阵处理的方法。
3．掌握MATLAB 分析的方法。
实验过程及内容：
1．求下列表达式的值。
（1）w=2×(1+0.34245×10-6)
（2）a
c b abc c b a x ++-+++
=)tan(e 22ππ，其中a=3.5，b=5，c=-9.8。
（3），])4
8333.0()41[(22απβπαy ---=π其中α=3.32，β=-7.9。
（4）)1ln(2122t t e z t ++=，其中t =⎥⎦⎤⎢⎣⎡--65.05i 312。
2．已知
A = ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--7613870451，
B =⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡--023352138
求下列表达式的值：
（1）A+6B 和A 2-B+I （其中I 为单位矩阵）。
（2）A*B 、A.*B 和B*A 。
（3）A/B 及B\A 。
（4）[A,B]和[A([1,3],:);B^2]。
3．已知
⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---=14.35454
.9632053256545410778.01023Α 完成下列操作：
（1）输出A 在[10，25]范围内的全部元素。
（2）取出A 前3行构成矩阵B ，前两列构成矩阵C ，右下角3×2子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。
（3）分别求表达式E<D 、E&D 、E|D 和～E ｜～D 的值。
4．产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好，为什么？
5．已知：
⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡--=5881252018629A
求A的特征值及特征向量，并分析其数学意义。
实验内容
1.（1）w=sqrt(2)*(1+0.34245*10^(-6))
w =
1.4142
1.（2）
a=3.5
b=5;
c=-9.8;
x=(2*pi*a+((b+c)/(pi+a*b*c))-exp(2))/(tan(b+c)+a) x =
0.9829
1.（3）
>> d=3.32;
>> e=-7.9;
>> y=2*pi*d^2*((1-pi/4)*e-(0.8333-pi/4)*d)
y =
-128.4271
1.（4）
t=[2,(1-3i);5,-0.65];
z=1/2*exp(2*t)*log(t+sqrt(1+t^2))
z =
1.0e+004 *
0.0057 - 0.0007i 0.0049 - 0.0027i
1.9884 - 0.3696i 1.7706 - 1.0539i
2.（1）a=[-1,5,-4;0,7,8;3,61,7]
a =
-1 5 -4
0 7 8
3 61 7
>> b=[8,3,-1;2,5,3;-3,2,0]
b =
8 3 -1
2 5 3
-3 2 0
>> i=[1,0,0;0,1,0;0,0,1]
i =
1 0 0
0 1 0
0 0 1
>> r=a+6*b
r =
47 23 -10
12 37 26
-15 73 7
>> t=a^2-b+i
t =
-18 -217 17
22 533 109
21 867 526
2.（2）l=a*b
l =
14 14 16
-10 51 21
125 328 180
>> k=a.*b
k =
-8 15 4
0 35 24
-9 122 0
>> m=b*a
m =
-11 0 -15
7 228 53
3 -1 28
2.（3）f=a/b
f =
1.2234 -0.9255
2.9787
-0.9468 2.3511 -0.9574
4.6170 3.8723 13.8936 >> g=b\a
g =
-0.5106 -8.6170 -1.1277
0.7340 17.5745 1.8085
-0.8830 -21.2128 0.4043
2.（4）
v=[a,b]
v =
-1 5 -4 8 3 -1
0 7 8 2 5 3
3 61 7 -3 2 0
>> n=[a([1,3],:);b^2]
n =
-1 5 -4
3 61 7
73 37 1
17 37 13
-20 1 9
3（1）a=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]
a =
23.0000 10.0000 -0.7780 0
41.0000 -45.0000 65.0000 5.0000
32.0000 5.0000 0 32.0000
6.0000 -9.5400 54.0000 3.1400
>> f=find((a>=10)&(a<=25))
f =
1
5
>> a(f)
ans =
23
10
3.（2）b=a(1:3,:)
b =
23.0000 10.0000 -0.7780 0
41.0000 -45.0000 65.0000 5.0000
32.0000 5.0000 0 32.0000 c=a(:,1:2)
c =
23.0000 10.0000
41.0000 -45.0000
32.0000 5.0000
6.0000 -9.5400 d=a(2:4,3:4)
d =
65.0000 5.0000
0 32.0000
54.0000 3.1400 e=b*c
e =
1.0e+003 *
0.9141 -0.2239
1.2080
2.7123
1.1330 -0.2103
3.（3）e<d
ans =
0 1
0 0
0 1
>> e&d
ans =
1 1
0 1
1 1
e|d
ans =
1 1
1 1
1 1
>> ~e|~d
ans =
0 0
1 0
0 0
4.
format rat
>> h=hilb(5)
h =
1 1/
2 1/
3 1/
4 1/
5 1/2 1/3 1/4 1/5 1/
6 1/3 1/4 1/5 1/6 1/
7 1/4 1/5 1/6 1/7 1/
8 1/5 1/6 1/7 1/8 1/9
>> p=pascal(5)
p =
1 1 1 1 1 1
2
3
4
5 1 3
6 10 15
1 4 10 20 35
1 5 15 35 70
Hh=det(h)
Hh =
1/266716800000
>> Hp=det(p)
Hp =
1
>> th=cond(h)
th =
476607
>> tp=cond(p)
tp =
178868/21
>>
5阶帕斯卡矩阵P 比5阶希尔伯特矩阵H矩阵性能更好。
因为th>tp，所以值更小的的更稳定。
5．>> a=[-29,6,18;20,5,12;-8,8,5];
>> [v,d]=eig(a)
v =
698/979 395/1409 989/3619 -859/1412 -2430/3089 349/400 538/1543 423/769 578/1427 d =
-3595/142 0 0 0 -3755/357 0 0 0 4697/279 由上述可知特征值为
λ1= -3595/142 ； λ2 = -3755/357 ； λ3=4697/279 意义;是用来求0)(=-ζλI A 时，有非零解ζ。
即|)(I A λ-|=0.
实验结论：
指导教师批阅意见：
成绩评定：
指导教师签字：
年月日备注：
注：1、报告内的项目或内容设置，可根据实际情况加以调整和补充。
2、教师批改学生实验报告时间应在学生提交实验报告时间后10日内。