matlab习题第一章
- 格式:doc
- 大小:67.50 KB
- 文档页数:3
1.输入20个数,求其中最大数和最小数。
要求分别用循环结构和调用MA TLAB的max函数、min函数来实现。
a=rand(1,20);max1=max(a) %用max函数求最大值min1=min(a) %用max函数求最大值%以下用循环方法求解:max2=a(1);min2=a(1);for i=2:20if a(i)>max2max2=a(i);endif a(i)<min2min2=a(i);endendmax2min2调用MATLAB的max函数、min函数来实现:>> a=rand(1,20)zuidazhi=max(a)zuixiaozhi=min(a)a =Columns 1 through 80.7513 0.2551 0.5060 0.6991 0.8909 0.9593 0.5472 0.1386 Columns 9 through 160.1493 0.2575 0.8407 0.2543 0.8143 0.2435 0.9293 0.3500 Columns 17 through 200.1966 0.2511 0.6160 0.4733zuidazhi =0.9593zuixiaozhi =0.13862. 求Fibonacci数列(1)大于4000的最小项。
(2)5000之内的项数。
(1)function t=fibonacc1f=[1,1];n=2;while 1f=[f,f(n-1)+f(n)];n=length(f);if f(n)>4000break;endendt=f(n);return(2)function n=fibonacc2f=[1,1];n=2;while 1f=[f,f(n-1)+f(n)];n=length(f);if f(n)>4000break;endendn=n-1;return3. 写出下列程序的输出结果:s=0;a=[12,13,14;15,16,17;18,19,20;21,22,23];for k=afor j=1:4if rem(k(j),2)~=0s=s+k(j);endendendSAns=1081、解方程组Ax =b ,分别用求逆解法与直接解法求其解。
习题 11. 执行下列指令,观察其运算结果, 理解其意义: (1) [1 2;3 4]+10-2i(2) [1 2; 3 4].*[0.1 0.2; 0.3 0.4] (3) [1 2; 3 4].\[20 10;9 2] (4) [1 2; 3 4].^2 (5) exp([1 2; 3 4]) (6)log([1 10 100]) (7)prod([1 2;3 4])(8)[a,b]=min([10 20;30 40]) (9)abs([1 2;3 4]-pi)(10) [1 2;3 4]>=[4,3;2 1](11)find([10 20;30 40]>=[40,30;20 10])(12) [a,b]=find([10 20;30 40]>=[40,30;20 10]) (提示:a 为行号,b 为列号) (13) all([1 2;3 4]>1) (14) any([1 2;3 4]>1) (15) linspace(3,4,5) (16) A=[1 2;3 4];A(:,2)2. 执行下列指令,观察其运算结果、变量类型和字节数,理解其意义: (1) clear; a=1,b=num2str(a),c=a>0, a= =b, a= =c, b= =c (2) clear; fun='abs(x)',x=-2,eval(fun),double(fun)3. 本金K 以每年n 次,每次p %的增值率(n 与p 的乘积为每年增值额的百分比)增加,当增加到rK 时所花费的时间为)01.01ln(ln p n rT +=(单位:年)用MA TLAB 表达式写出该公式并用下列数据计算:r =2, p =0.5, n =12.4.已知函数f (x )=x 4-2x 在(-2, 2)内有两个根。
取步长h =0.05, 通过计算函数值求得函数的最小值点和两个根的近似解。
阅读使人快乐,成长需要时间
江阴室内设计培训学校排名
一、江阴问鼎教育(澄江中路5号东都大厦8楼)
培训方式:小班制授课,10人以内,高端室内设计培训机构;
影响力:江阴高性价比的室内设计培训品牌;
特点:环境特优雅,教学质量上乘;
上课模式:小班化教学,不限学时,包学包会;
不足:广告做的不多,主要靠口碑推广
面向对象:成人、学生。
二、江阴逸仙教育
培训方式:面授为主,辅助在线预习等学习;
影响力:一家留学生开的机构,广告打的挺多的;
特点:课堂比较漂亮,环境比较好。
学习时间自由,要自制力比较好的学员;
不足:价格比较贵,要1万以上;课时数比较少;
上课模式:课堂面授一对4+大班课+在线学习;
面对对象:成人。
三、远大培训机构
影响力:成人培训的影响力不如上面几家培训做的那么闻名遐迩,不过品牌效用还是很好;特点:老师口才比较好,上课气氛比较好,价钱也还好;
不足:老生常谈,创新度没其他机构高,上课人数比较多;
师资:兼职老师;
上课模式:课堂面授、大班为主;
面对对象:大学生、中学生。
四、英华培训机构
培训方式:课堂面授;
影响力:打入江苏的培训机构,知名度还可以;
特点:全部学员使用3阶段学习法,自学可以借鉴这个方法。
学习时间自由,要人监督才行;不足:课程设置比较制式化,主题比较少,然后也比较贵,1万以上了;
大概描述:提供了三阶段的学习方法;
师资:中教、外教;
上课模式:课堂面授、一对4;
面对对象:企业和个人。
matlab章课后作业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.4410i 4.计算如下表达式:(1)(2)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i >> sin(2-8*i) 1.3553e+003 +6.2026e+002i5.判断下⾯语句的运算结果。
欢迎共阅习题二1.如何理解“矩阵是MATLAB 最基本的数据对象”? 答:因为向量可以看成是仅有一行或一列的矩阵,单个数据(标量)可以看成是仅含一个元素的矩阵,故向量和单个数据都可以作为矩阵的特例来处理。
(2) 删除矩阵A 的第7号元素。
答:A(7)=[](3) 将矩阵A 的每个元素值加30。
答:A=A+30;(4) 求矩阵A 的大小和维数。
答:size(A);ndims(A);(5) 将向量 t 的0元素用机器零来代替。
答:t(find(t==0))=eps; (6) 将含有12个元素的向量 x 转换成34⨯矩阵。
答:reshape(x,3,4); (7) 求一个字符串的ASCII码。
答:abs(‘123’); 或E 。
答:B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4);E=B*C;(2) 分别求E<D 、E&D 、E|D 、~E|~D 和find(A>=10&A<25)。
答:E<D=010001⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,E&D=110111⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦,11⎡⎤⎢⎥答:student(1).id='0001';student(1).name='Tom'; student(1).major='computer';student(1).grade=[89,78,67,90,86,85]; 8.建立单元矩阵B 并回答有关问题。
B{1,1}=1;B{1,2}='Brenden';B{2,1}=reshape(1:9,3,3); B{2,2}={12,34,2;54,21,3;4 ,23,67};(1)size(B)和ndims(B)的值分别是多少?答:size(B) 的值为2, 2。
ndims(B) 的值为2。
(2)B(2)和B(4)的值分别是(2)建立5×6随机矩阵A,其元素为[100,200]范围内的随机整数。
第1章一、思考题4. (1) B=A(2:5,1:2:5)(2)A(7)=[](3) A=A+30(4)size(A)(5) t(find(t==0))=eps(6)t=reshape(x,3,4)(7)abs('matlab')(8)char(93)5. A=[97,67,34,10;-78,75,65,5;32,5,-23,-59;0,-12,54,7](1) B=A(1:3,:)C=A(:,1:2)D=A(2:4,3:4)E=B*C(2)A(find(A>=50 & A<=100))二、实验题1.x=-74/180*pi;y=-27/180*pi;sin(x*x+y*y)/sqrt(tan(abs(x+y)))+pi2.a=-3.0:0.1:3.0;exp(-0.3*a).*sin(a+0.3)3.x=[2,4;-0.45,5];log(x+sqrt(1+x.*x))/24. A=[3,54,2;34,-45,7;87,90,15];B=[1,-2,67;2,8,74;9,3,0]; (1)A*Bans =129 432 41977 -407 -1052402 591 12489A.*Bans =3 -108 13468 -360 518783 270 0(2) A^3ans =-28917 240246 -4368137883 -259101 27669171333 252504 38673A.^3ans =27 157464 839304 -91125 343658503 729000 3375(3)A/Bans =-5.8845 5.3549 -0.20286.3554 -5.6596 4.3293-6.6325 6.2078 9.0241B/Aans =30.2855 19.2643 -8.561433.4394 21.1547 -9.3974-0.7443 -0.3938 0.2830(4)[A,B]ans =3 54 2 1 -2 6734 -45 7 2 8 7487 90 15 9 3 0[A([1,3],:);B^2]ans =3 54 287 90 15600 183 -81684 282 72615 6 8255.a=1+2i;b=3+4i;c=exp(pi*i/6)c =0.8660 + 0.5000ic+a*b/(a+b)ans =1.6353 + 1.8462i第2章一、思考题3.s=0;for n=0:63s=s+2^n;enddisp(s)n=0:63;s=sum(2.^n)二、实验题1.x=input('输入一个四位整数:');y=[fix(x/1000),mod(fix(x/100),10),mod(fix(x/10),10),mod(x,10)] z=mod((y+7),10)x=z(3)*1000+z(4)*100+z(1)*10+z(2)2.gh=input('输入工号');h=input('输入工时');dj=84;if h>120gz=dj*120+1.15*dj*(h-120);elseif h<60gz=dj*h-700;elsegz=dj*h;endformat bank;display([gh,gz])3.循环结构n=input('input n:');s=0;for k=1:ns=s+1/k^2;enddisplay(sqrt(s*6))向量运算n=input('input n:');k=1:n;display(sqrt(sum(1./k.^2)*6))4.y=0;k=0;while y<3k=k+1;y=y+1/(2*k-1);enddisplay([k-1,y-1/(2*k-1)])5.x0=0;x=1;k=0;a=input('a=');b=input('b=');while abs(x-x0)>=1e-5 && k<500x0=x;x=a/(b+x0);k=k+1;enddisplay([k,x]);display([(-b+sqrt(b^2+4*a))/2,(-b-sqrt(b^2+4*a))/2]);6.y=fun(40)/(fun(30)+fun(20))(1)函数文件fun.mfunction f=fun(n)f=n+log(n^2+5);(2)函数文件fun.mfunction f=fun(n)a=1:n;f=sum(a.*(a+1));第3章一、思考题4.t=0:0.01:1.5;y=sqrt(3)/2*exp(-4*t).*sin(4*sqrt(3)*t+pi/3);5.x=-10:0.01:10;y=linspace(-6,6,size(x,2))z=x.^3+3*x.*y.^2;plot3(x,y,z)6.x=100:100:400;y=100:100:400;z=[636,697,624,478;698,712,630,478;680,674,598,412;662,626,552,334]; [X,Y]=meshgrid(x,y);mesh(X,Y,z)二、实验题1.(1)x=-10:0.01:10;plot(x,x-x.^3/3/2)(2)plot(x,exp(-x.*x/2)/2/pi)(3)x=-8:0.01:8;plot(x,sqrt((64-x.*x)/2))(4)t=0:pi/100:6*pi;plot(t.*sin(t),t.*cos(t))2. (1)x1=linspace(0,1,100);y1=2*x1-0.5;t=linspace(0,pi,100);x=sin(3*t).*cos(t);y=sin(3*t).*sin(t);plot(x1,y1,'r-',x,y,'b:');text(0.4,1,'y=2x-0.5');text(-0.4,0.5,'x=sin(3t)cos(t)');text(-0.4,0.3,'y=sin(3t)sin(t)');(2)subplot(1,2,1);scatter(x1,y1,10)title('y=2x-0.5');subplot(1,2,2);scatter(x,y,10)3.subplot(1,21);x=1:1:100; y=sin(1./x);plot(x,y)subplot(1,2,2);fplot('sin(1/x)',[1,100])4.subplot(2,2,1);bar(t,y);subplot(2,2,2);stairs(t,y);subplot(2,2,3);stem(t,y)subplot(2,2,4);semilogy(t,y);5.theta=linspace(-pi,pi,100);ro=5.*cos(theta)+4;polar(theta,ro);fi=linspace(0,2*pi,100);a=1r=a.*(1+cos(fi));polar(fi,r);6.t=0:pi/20:2*pi;x=exp(-t./20).*cos(t);y=exp(-t./20).*sin(t);z=t; plot3(x,y,z);t=0:0.01:1;x=t;y=t.^2;z=t.^3;plot3(x,y,z);7.x=-30:0.1:0;y=0:0.1:30;[x,y]=meshgrid(x,y);z=10.*sin(sqrt(x.^2+y.^2))./sqrt(1+x.^2+y.^2); meshc(x,y,z);8.x=linspace(-3,3,100);y=linspace(-3,3,100);[x y]=meshgrid(x,y);fxy=-5./(1+x.^2+y.^2);i=find(abs(x)<=0.8 & abs(y)<=0.5);fxy(i)=NaN;surf(x,y,fxy)9.u=linspace(1,10,100);v=linspace(-pi,pi,100);[u,v]=meshgrid(u,v);x=3.*u.*sin(v);y=2.*u.*cos(v);z=4*u.^2;x=3*u.*sin(v);y=2*u.*cos(v);z=4*u.^2;surf(x,y,z);shading interp;light('position',[1,0,1]);10.t=0:pi/100:2*pi;y=sin(t);comet(t,y)第4章一、思考题5.(1)A=eye(3);(2)C=100+(200-100)*rand(5,6);(3)D=1+sqrt(0.2)*randn(1,500);(4)E=ones(size(A));(5)A=A+30*eye(size(A));(6)B=diag(diag(A))二、实验题1.P=pascal(5);H=hilb(5);Dp=det(P);Dh=det(H);Kp=cond(P);Kh=cond(H);P矩阵的性能更好,因为Kp较小2.A=[1,-1,2,3;0,9,3,3;7,-5,0,2;23,6,8,3]B=[3,pi/2,45;32,-76,sqrt(37);5,72,4.5e-4;exp(2),0,97]A1=diag(A);B1=diag(B);A2=triu(A);B2=triu(B);A3=tril(A);B3=tril(B);rA=rank(A);rB=rank(B);nA=norm(A);nb=norm(B);cA=cond(A);cB=cond(B);3.A=[31,1,0;-4,-1,0;4,-8,-2];[V,D]=eig(A);4.A=diag([-1,-1,-1,-1],-1)+diag([-1,-1,-1,-1],1)+diag([2,2,2,2,2]) b=[1,0,0,0,0]';x1=inv(A)*b;x2=A\b;[L,U]=lu(A);x3=U\(L\b);[Q,R]=qr(a);[Q,R]=qr(A);x4=R\(Q\b)R=chol(A);x5=R\(R'\b)5.B=sparse(A);x1=inv(B)*b;x2=B\b;[L,U]=lu(B);x3=U\(L\b);第5章一、思考题3.A=randn(10,5);mean(A)std(A)max(max(A))min(min(A))sum(A,2)sum(sum(A))sort(A,1)sort(A,2,'descend')二、实验题1.A=rand(1,30000);mean(A)std(A)max(A)min(A)size(find(A>0.5))/size(A)2.h=[466,715,950,1422,1635];w=[7.04,4.28,3.40,2.52,2.13];hh=[500,900,1500];ww=interp1(h,w,hh,'spline')3.x=linspace(1,10,50);y=log(x);f=polyfit(x,y,5);yy=polyval(f,x);plot(x,y,'r-',x,yy,'g.')4.N=64; % 采样点数T=5; % 采样时间终点t=linspace(0,T,N); % 给出N个采样时间ti(I=1:N) x=3*exp(-t); % 求各采样点样本值xdt=t(2)-t(1); % 采样周期f=1/dt; % 采样频率(Hz)X=fft(x); % 计算x的快速傅立叶变换XF=X(1:N/2+1); % F(k)=X(k)(k=1:N/2+1) f=f*(0:N/2)/N; % 使频率轴f从零开始plot(f,abs(F),'-*') % 绘制振幅-频率图xlabel('Frequency');ylabel('|F(k)|')5.(1)p1=[1 2 0 0 7];p2=[1 -2];p3=[1 0 0 5 1];p12=conv(p1,p2);p=p12+[zeros(1,size(p12,2)-size(p3,2)),p3];roots(p)(2)A=[-1,4,3;2,1,5;0,5,6];Px=polyval(p,A)Pxm=polyvalm(p,A)6.(1)z=fzero('3*x-sin(x)+1',0)(2)建立函数文件myfun.mfunction F=myfun(X)x=X(1);y=X(2);F(1)=x*x+y*y-9;F(2)=x+y-1;在命令窗口中输入以下命令:x=fsolve(@myfun,[3,0]',optimset('Display','off')) 第6章一、思考题2.fx=inline('1./(1+x.^2)');[I,n]=quad(fx,-100000,100000,1e-3);[I,n]=quadl(fx,-100000,100000,1e-3);x=-100000:0.01:100000;y=1./(1+x.*x);trapz(x,y);3.(1)fx=inline('-2*y+2*x*x+2*x');[t,y]=ode23(fx,[0,0.5],1)(2)fx=inline('y-exp(x)*cos(x)');[t,y]=ode23(fx,[0,3],1)二、实验题1.for x=1:3fx=[x,x^2,x^3;1,2*x,3*x;0,2,6*x];diff(fx)end2.(1)x=0:0.01:1;y=x.^10+10.^x+1./log10(x);dy=diff(y)/0.01;(2)x=0:0.01:1;y=log(1+x);dy=diff(y,2)/0.01;plot(x(1:99),dy)3.(1)fx=inline('x.^2.*sqrt(2*x.*x+3)');quad(fx,1,5)(2)fx=inline('x./sin(x).^2');quad(fx,pi/4,pi/3)(3)fx=inline('abs(cos(x+y))');dblquad(fx,0,pi,0,pi)(4)syms x y ;fx=x*y;int(int(fx,y^2,y+2),-1,2)x 的积分区间为【0,2】时fx=inline('x.*y');dblquad(fx,0,2,-1,2)4. x=0.3:0.2:1.5;y=[0.3895,0.6598,0.9147,1.1611,1.3971,1.6212,1.8325]; trapz(x,y)5.(1)yp=inline('-(1.2+sin(10*x))*y');[t,y]=ode23(yp,[0,5],1);(2)令''3,,21y x y x y x ='==,则可写出原方程的状态方程形式: 123233221sin 31)1(2cos 5cos x t x x t t t x x x x x +--++='='='⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡+⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡+-+-=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡'''t x x x t t tx x x cos 00)1(2cos 51sin 311000103212321 建立函数文件ztfun.m function y=ztfun(t,x)b=[0;0;cos(t)];y=[0,1,0;0,0,1;-1/(3+sin(t)),-1,5*cos(2*t)/(t+1)^2]*x+b; 解微分方程组[t,y]=ode23(@ztfun,[0,5],[1;0;2]);6.建立函数文件ztfun.mfunction yy=ztfun(t,y) yy=[y(2)*y(3);-y(1)*y(3);-0.51*y(1)*y(2)];解微分方程组 [t,y]=ode23(@ztfun,[0,5],[0;1;1])第7章一、思考题3.(1)数值积分fx=inline('exp(x).*(1+exp(x)).^2');quad(fx,0,log(2))符号积分f=sym('exp(x)*(1+exp(x))^2');v=int(f,0,log(2));eval(v)(2)略二、实验题1.A=sym('[1,2,3;x,y,z;3,2,1]')rank(A)inv(A)det(A)2.(1) y=sym('sqrt(x+sqrt(x+sqrt(x)))');y1=diff(y)y2=diff(y,'x',2)(2) syms x y;fxy=sin(x^2*y)*exp(-x^2-y);diff(diff(fxy,x),y)3.(1) syms xint(1/(1+x^4))(2) syms x tint((-2*x*x+1)/(2*x*x-3*x+1)^2,x,cos(t),exp(2*t)) 4.syms n xsymsum(1/(2*n+1)/(2*x+1)^(2*n+1),n,0,inf)symsum(1/(2*n+1)/(2*x+1)^(2*n+1),n,0,5)5.(1) syms xtaylor((exp(x)+exp(-x))/2,5,0)(2) syms a xtaylor(exp(-5*x)*sin(3*x+pi/3),5,a)6.(1)x=solve(sym('x^3+a*x+1=0'))(2)[x y]=solve('sqrt(x^2+y^2)-100=0,3*x+5*y-8=0')7. 方程转化为: π21'11't y y y y -=+=⎩⎨⎧符号解[y1,y11]=dsolve('Dy=y1,Dy1+y=1-t^2/pi','y(-2)=5,y1(-2)=5','t') 数值解编写函数文件ztfun.mfunction yy=ztfun(t,y)yy=[y(2);1-t^2/pi-y(1)];在命令窗口输入以下命令[t,y]=ode45(@ztfun,[-2,7],[-5;5]);t=linspace(-2,7,49)y2=y8.[x,y]=dsolve('Dx=3*x+4*y,Dy=-4*x+3*y','x(0)=0,y(0)=1')。
第一次上机作业1.设A和B是两个同维同大小的矩阵,问:1)A*B和A.*B的值是否相等?不相等2)A./B和B.\A的值是否相等?相等3)A/B和B\A的值是否相等?不相等4)A/B和B\A所代表的数学含义是什么?A\B等效于A的逆左乘B矩阵,也就是inv(A)*B,而B/A等效于A矩阵的逆右乘B矩阵,也就是B*inv(A)。
2.写出完成下列操作的命令:1)将矩阵A的第2~5行中的第1,3,5列元素赋给矩阵BB=A(:,1:3:5)2)删除矩阵A的第7号元素A(7)=[]3)将矩阵A的每个元素值加30A=A+304)求矩阵A的大小和维数whos和size(A)5)将含有12个元素的向量x转换成3×4矩阵y=reshape(x,3,4)3.下列命令执行后,L1、L2、L3、L4的值分别是多少?A=1:9; B=10-A;L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);L1 = 0 0 0 0 1 0 0 0 0 L2 = 1 1 1 1 1 0 0 0 0 L3 =0 0 0 1 1 1 0 0 0 L4 = 4 5 6 4.已知:A=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---14.35454.9632053256545410778.01023完成下列操作:(1)取出A 的前3行构成矩阵B ,前两列构成矩阵C ,右下角3×2子矩阵构成矩阵D ,B 与C 的乘积构成矩阵E 。
B=A(1:3,:) C=A(:,1:2) D=A(2:4,3:4) E=B*C(2)分别求E<D 、E&D 、E|D 、~E|~D 和find(A>=10&A<25)。
1)E<D 结果 0 1 0 0 0 1 2)E&D 结果 1 1 0 1 1 1 3)E|D 结果 1 1 1 1 1 1 4) ~E|~D 结果 0 0 1 0 0 05)find(A>=10&A<25)结果155.当A=[34, NaN, Inf,-Inf,-pi,eps,0]时,分析下列函数的执行结果:all(A), any(A), innan(A), isinf(A), isfinite(A)(1)all(A)=0,矩阵A不是所有的元素都非零。
Matlab程序设计课程实验目录Matlab操作实验(第1章)实验指导思考练习矩阵运算实验(第2章)实验指导思考练习程序设计实验(第3章)实验指导思考练习高级绘图实验(第4章)实验指导思考练习数据分析与多项式计算实验(第5章)实验指导思考练习解方程与最优化问题求解实验(第6章)实验指导思考练习数值运算实验(第7章)实验指导思考练习符号实验(第8章)实验指导思考练习图形句柄实验(第9章)实验指导思考练习图形用户界面设计实验(第10章)实验指导思考练习Simulink动态仿真实验(第11章)实验指导思考练习外部接口技术实验(第12章)实验指导思考练习1、Matlab操作实验一、实验内容1. 第一章实验内容1:path(path,’系统盘名:\matlab文件名\要添加的文件名’)如:path(path,’c:\matlab701\my work’)实验内容-2例1-1:x=[0:0.5:360]*pi/180;plot(x,sin(x),x,cos(x));例1-2p=[3,7,9,0,-23];x=roots(p)例1-3quad(‘x.*log(1+x)’,0,1)例1-4:实验内容3Help 函数名如:help inv help plot等实验内容4(1)x=0:pi/10:2*pi;y=sin(x);二、思考练习思考题1启动:有三种常见方法,(1)、在Windows桌面,单机任务栏上的‚开始‛按钮,选择‚所有程序‛菜单项,然后选择MATLAB程序组中的MATLAB (版本号)程序选项,即可启动。
(2)、在MATLAB的安装路径中找到MATLAB系统启动程序matlab.exe,然后运行它。
(3)、利用建立快捷方式功能,将MATLAB系统启动程序以快捷方式的形式放在Windows桌面上,以后只要在桌面上双击该图标即可启动MATLAB。
关闭:常见的也有3种方法,(1)、在MATLAB主窗口file菜单中选择exit MATLAB命令。
目录第一章MA TLAB入门 (1)1.1 MA TLAB的优点 (1)1.易用性 (1)2.平台独立性 (1)3.预定义函数 (1)4.机制独立的画图 (1)5.用户图形界面 (2)6.MA TLAB编译器 (2)1.2 MA TLAB的缺点 (2)1.3 MA TLAB的开发环境 (2)1.3.1MA TLAB桌面 (2)1.3.2命令窗口(TheCommandWindow) (2)1.3.3历史命令窗口(The History Command Window) (4)1.3.4启动平台(the launch pad) (4)1.3.5图像窗口(Figure Windows) (5)1.3.6MA TLAB工作区 (5)1.3.7MA TLAB帮助 (7)1.3.8一些重要的命令 (8)1.4把MA TLAB当作便笺薄来使用 (10)测试1.1 (10)1.5总结 (11)1.5.1MA TLAB总结 (11)1.6练习 (11)1.1 (11)1.2 (11)1.3 (11)1.4 (11)1.5 (12)1.6 (12)1.7 (12)1.8 (12)第一章MATLAB入门MATLAB(矩阵实验室的简称)是一种专业的计算机程序,用于工程科学的矩阵数学运算。
但在以后的几年内,它逐渐发展为一种极其灵活的计算体系,用于解决各种重要的技术问题。
Matlab程序执行MATLAB语言,并提供了一个极其广泛的预定义函数库,这样就使得技术工作变得简单高效。
本书将介绍matlab语言,并向大家展示如何运用它去解决经典的技术问题。
MA TLAB是一个庞大的程序,拥有难以置信的各种丰富的函数;即使基本版本的MA TLAB语言拥有的函数也比其他的工程编程语言要丰富的多。
基本的MA TLAB语言已经拥有了超过1000多个函数,而它的工具包带有更多的函数,由此扩展了它在许多专业领域的能力。
本书无意将MA TLAB的所有函数介绍给大家,而是让大家掌握编写调试和优化程序的基本功,还有一些重要函数的子集。
第一章 实验内容
1.
2.例1-1 x=[0:0.5:360]*pi/180;
plot(x,sin(x),x,cos(x));
例1-2 p=[3,7,9,0,-23];
x=roots(p)
x =
-1.8857
-0.7604 + 1.7916i
-0.7604 - 1.7916i
1.0732
例1-3 quad('x.*log(1+x)',0,1)
ans =
0.2500
例 1-4 a=[2,-3,1;8,3,2;45,1,-9];
b=[4;2;17];
x=inv(a)*b
x =
0.4784
-0.8793
0.4054
3.
(1)help inv
inv Matrix inverse.
inv(X) is the inverse of the square matrix X.
A warning message is printed if X is badly scaled or
nearly singular…
(2)help plot
plot Linear plot.
plot(X,Y) plots vector Y versus vector X. If X or Y is a matrix,
then the vector is plotted versus the rows or columns of the matrix,
whichever line up. If X is a scalar and Y is a vector, disconnected
line objects are created and plotted as discrete points vertically at
X…
(3)help max
max Largest component.
For vectors, max(X) is the largest element in X. For matrices,
max(X) is a row vector containing the maximum element from each
column. For N-D arrays, max(X) operates along the first
non-singleton dimension...
(4)help round
round Round towards nearest integer.
round(X) rounds the elements of X to the nearest integers.
4.(1) x=0:pi/10:2*pi;
y=sin(x)
y =
Columns 1 through 10
0 0.3090 0.5878 0.8090 0.9511 1.0000 0.9511
0.8090 0.5878 0.3090
Columns 11 through 20
0.0000 -0.3090 -0.5878 -0.8090 -0.9511 -1.0000 -0.9511
-0.8090 -0.5878 -0.3090
Column 21
-0.0000
(2)根据提示操作,选择y变量,并绘制y图形,生成了如下图形,图形是一个正弦函数
图像。
思考练习
1.与其他窗口程序的启动与退出类似,双击图标打开,单击X图标关闭
2. MATLAB是矩阵实验室(Matrix Laboratory)的简称用于算法开发、数据可视化、数据分析
以及数值计算的高级技术计算语言和交互式环境。它将数值分析、矩阵计算、科学数据可视
化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中,为
科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方
案。
3.如果一个命令行很长,需要分成多行输入,可以在第一个物理行之后加上三个小黑点…并
按下enter键,然后接着在下一个物理航继续输入命令的其他部分。三个小黑点称为续行符,
即把下面的物理行看作该行的逻辑继续。
4.help命令只搜索出那些与关键字完全匹配的结果,lookfor命令对搜索范围内的小文件进行
关键字搜索,条件比较宽松。
5.假如建立了一个变量fac,同时在当前目录下建立了一个M文件fac.m文件。如果在命令
窗口输入fac,按照MATLAB的搜索顺序,应该在屏幕上显示fac的值。如果没有建立fac变量,
则执行fac.m文件。