MATLAB基础操作(二)

11.求方程2x^5-3x^3+71x^2-9x+13=0的全部跟
    >> p=[2,0,-3,71,-9,13];
    >> x=roots(p);

12.求解线性方程组2x+3y-z=2
         8x+2y+3z=4
         45x+3y+9z=23
    >> a=[2,3,-1;8,2,3;45,3,9];%建立系数矩阵a
    >> b=[2,4,23]%建立列向量b
    >> x=inv(a)*b

13.求解积分0->1xln(1+x)
    >> quad('x.*log(1+x)',0,1)

14.多项式曲线拟合
    >> x=[1,2,3,4,5,6,7,8,9,10];
    >> y=[1,2,3,4,5,4.7,5,5.2,6,7.2];
    >> p1=polyfit(x,y,1)%一次多项式拟合
    >> p3=polyfit(x,y,3)%三次多项式拟合
    >> x2=1:0.1:10;
    >> y1=polyval(p1,x2);
    >> y3=polyval(p3,x2);
    >> plot(x,y,'*',x2,y1,':',x2,y3)

15.函数linspace(x1,x2,n)
    >> x=linspace(1,2,5)

    x =1.0000    1.2500    1.5000    1.7500    2.0000

16.ones(生成全1的矩阵)、zeros（生成全0的矩阵）
    >> ones(3)
    >> ones(3,4)

17.函数reshape()
    >> a=-4:4
    >> reshape(a,3,3)

18.画出y=1/(x+1)的函数曲线，x∈[0,100]
    >> x=0:100;
    >> y=1./(x+1);
    >> plot(x,y);
    >> legend('y=1/(x+1)');
    
19.生成一个信号x=sin(2*pi*t)+sin(4*pi*t)
    >> t=[0:199]./100;%采样时间点
    >> x=sin(2*pi*t)+sin(4*pi*t);%生成信号
    >> plot(t,x);
    >> legend('sin(2*pi*t)+sin(4*pi*t)');

20.判断方阵中的元素能否被3整除
    >> A=[24,35,13;23,39,47;90,41,80];
    >> P=rem(A,3)==0%被3整除，求余