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

百度文库 - 让每个人平等地提升自我 Matlab 上机实验一、二 3x 4y 7 z 12w 4 5x 7 y 4z 2w 4 3. 求下列联立方程的解 x 8z 5w 9 6x 5 y 2z 10 w 4 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 1 4 8 13 5 4 3 2 4．设A3 6 5 9 ， B 6 2 3 8 ，求 C1=A*B’ ;C2=A’ *B;C3=A.*B, 并求上述所 2 7 12 8 1 3 9 7 有方阵的逆阵。 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．设 3sin x ，把 π间分为 点，画出以 为横坐标， 为纵坐标的 y cos x 0.5 (1 x 2 ) x=0~2 101 x y 曲线。 x=linspace(0,2*pi,101); y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); plot(x,y,r) 6．产生 8×6 阶的正态分布随机数矩阵 R1, 求其各列的平均值和均方差。并求该矩阵全体数的 平均值和均方差。 （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) 1 百度文库 - 让每个人平等地提升自我 或者 u=reshape(a,1,48); p1=mean(u) p2=var(u) 设 x=rcost+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,y3,m) 设 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) 用 roots 函数求此多项式的根。 t=[1 0 -4 3 -2 6] p=roots(t) 建立一个字符串向量，然后对该向量做如下处理： 取第 1～5 个字符组成的子字符串。 将字符串倒过来重新排列。 统计字符串中小写字母的个数。 a=abcdef123456; subch=a(1:5) subch = abcde revch=a(end:-1:1) revch = 654321fedcba k=find(a=aa=z); a(k)=a(k)-(a-A); char(a) length(a) length(k) 10. 利用帮助查找 limit 函数的用法，并自己编写，验证几个函数极限的例子。 2 百度文库 - 让每个人平等地提升自我 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) returns -inf limit((sin(x+h)-sin(x))/h,h,0) returns cos(x) v = [(1 + a/x)^x, exp(-x)]; Matlab 上机实验二 验证所授课程的课件。 2. 设 x=sint, y=sin(nt+a) ， 若 a=1，令 n =1,2,3,4, 在四个子图中分别画出其曲线。 若 n=2, 取 a=0，π /3 ，π /2 ，及π，在四个子图中分别画出其曲线。 a=1; x=sin(t);y1=sin(1*t+a);y2=sin(2*t+a);y3=sin(3*t+a);y4=sin(4*t+a); subplot(2,2,1);plot(x,y1);sub