实验二 Matlab程序设计(答案).pdf

《实验二 Matlab程序设计(答案).pdf》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验二Matlab程序设计1.编写程序，计算1+3+5+7+…+(2n+1)的值（用input语句输入n值）。san1.m文件：n=input('n=');sum=0;fori=1:2:2*n+1sum=sum+i;endsum2.编写分段函数f(x)的函数文件，存放于文件ff.m中，计算出f(-3),f(sqrt(2)),f(∞)的值。ff.m文件：functiony=ff(x)if(x>=0)&&(x<1)y=x;elseif(x>=1)&&(x<=2)y=2-x;elsey=0;end3.利用矩阵的初等变

2、换及函数rank、函数inv求下列矩阵的秩：>>a=[1-632;3-540;-1-1124];rank(a)ans=34.计算（1）>>symsta=[cos(t)-sin(t);sin(t)cos(t)];a=det(eval(a^5));simplify(a)ans=1（2）b=[1210;6241;0210;3141];inv(b)ans=1.00000-1.00000-3.00001.00003.0000-1.00006.0000-2.0000-5.00002.0000-24.00007.000020

3、.0000-6.00005.用循环语句编写函数M文件计算e^x的值，其中x,n为输入变量（本实验取x=1,n=100），?121?的近似值可用下式表示：e≈1+?+?+⋯?.2!?!San5.m文件：functiony=san5(x)y=1;a=1;forn=1:100a=a*n;b=x^n/a;y=y+b;end程序：san5(1)ans=2.71836．（1）a=[4-22;-305;153];>>b=[134;-20-3;2-11];>>det(a)ans=-158（2）2*a-b,a*b,a.*b,a*

4、inv(b),inv(a)*b,a^2,a'ans=7-70-40130115ans=1210247-14-7-30-8ans=4-6860-152-53ans=-0.0000-0.00002.0000-2.7143-8.0000-8.14292.42863.00002.2857ans=0.48730.41141.00000.3671-0.43040.0000-0.10760.24680.0000ans=2424-7319-81336ans=4-31-2052537．已知多项式f(x),g(x)求：（1）f(x

5、)的根；（2）g(x)在闭区间[-1,2]上的最小值（3）f(x)+g(x),f(x)*g(x),f(x)/g(x)（4）f(x)的导数。>>symsx>>f=3*x^5-x^4+2*x^3+x^2+3;>>g=1/3*x^3+x^2-3*x-1;>>digits(5)>>solve(f)ans=-0.895180.7416*i+0.73304-1.0066*i-0.118791.0066*i-0.118790.73304-0.7416*ih=diff(g)h=x^2+2*x-3>>solve(h)ans=-3

6、1故为1>>f+gans=3*x^5-x^4+(7*x^3)/3+2*x^2-3*x+2>>f*gans=-(-x^3/3-x^2+3*x+1)*(3*x^5-x^4+2*x^3+x^2+3)>>f/gans=-(3*x^5-x^4+2*x^3+x^2+3)/(-x^3/3-x^2+3*x+1)>>d=diff(f)d=15*x^4-4*x^3+6*x^2+2*x