Matlab编程与应用习题和一些参考答案.docx

Matlab上机实验一、二 3.求下列联立方程的解 >> a=[3 4 -7 -12 ;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4．设 ， ，求 C1=A*B’;C2=A’*B;C3=A.*B, 并求 上述所有 方阵的逆阵。 >> A= [1 4 8 13;-3 6 -5 -9 ;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ，把 x=0 ~ 2 π 间分为101点，画出以 x 为横坐标， y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*s in(x)); >> plot(x,y,'r') 6．产生8 × 6阶的正态分布随机数矩阵 R 1 , 求其各列的平均值和均方差。并求 该矩阵 全体 数 的平均值和均方差。 （mean var） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g) 或者 u=reshape(a,1,48); p1= mean(u) p2= var(u) 7. 设 x=rc ost+3t,y=rsint+3, 分别令 r=2,3,4 ，画出参数 t=0~10 区间生成的 x~y 曲线。 >> t=linspace(0,10); >> r1=2; >> x1=(r1*cos(t)+3*t); >> y1=r1*sin(t)+3; >> r2=3; >> x2=(r2*cos(t)+3*t); >> y2=r2*sin(t)+3; >> r3=4; >> x3=(r3*cos(t)+3*t); >> y3=r3*sin(t)+3; >> plot(x1,y1,'r',x2,y2,'b',x3,y 3,'m') 8.设 f(x)= x 5 - 4x 4 +3x 2 - 2x+ 6 (1) 在 x=[-2,8] 之间取100个点，画出曲线，看它有几个过零点。（提示：用 polyval 函数） >> x=linspace(2,8,100); >> y=polyval([1 0 -4 3 -2 6],x); >> plot(x,y,'b',x,0,'y') (2) 用 roots 函数求此多项式的根。 t=[1 0 -4 3 -2 6] p=roots(t) 9.建立一个字符串向量，然后对该向量做如下处理： (1) 取第 1 ～ 5 个字符组成的子字符串。 (2) 将字符串倒过来重新排列。 ( 3 ) 统计字符串中小写字母的个数。 >> a='abcdef123456'; >> subch=a(1:5) subch = abcde >> revch=a(end:-1:1) revch = 654321fedcba>> k=find(a>='a'&a<='z'); >> a(k)=a(k)-('a'-'A'); >> char(a) >> length(a) >> length(k) 10.利用帮助查找 limit 函数的用法，并自己编写，验证几个函数极限的例子。 limit(sin(x)/x) returns 1 limit((x-2)/(x^2-4),2) returns 1/4 limit((1+2*t/x)^(3*x),x,inf) returns exp(6*t) limit(1/x,x,0,'right') returns inf limit(1/x,x,0,'left')