连杆机构运动分析&力分析Matlab语言m文件使用说明及算例

构件上点的运动分析

函数文件（m文件）

格式：function [ 输出参数] = 函数名（输入参数）

p_crank.m function [p_Nx,p_Ny]=p_crank(Ax,Ay,theta,phi,l1)

v_crank.m function [v_Nx,v_Ny]=v_crank(l1,v_Ax,v_Ay,omiga,theta,phi)

a_crank.m function [a_Nx,a_Ny]=a_crank(l1,a_Ax,a_Ay,alpha,omiga,theta,phi)

函数中的符号说明

函数文件（m 文件）

格式： function [ 输出参数 ] = 函数名（ 输入参数 ）

p_RRR.m function [cx,cy,theta2,theta3]=p_RRR(bx,by,dx,dy,l2,l3,m)

v_RRR.m function [vcx,vcy,omiga2,omiga3]=v_RRR(vbx,vby,vdx,vdy,cx,cy,bx,by,dx,dy)

a_RRR.m function [acx,acy,alpha2,alpha3]=a_RRR(abx,aby,adx,ady,cx,cy,bx,by,dx,dy,omiga2,omiga3)

函数中的符号说明

m =1 m = -1

RRR Ⅱ级杆组运动分析

函数文件（m 文件）

格式： function [ 输出参数 ] = 函数名（ 输入参数 ）

p_RRP.m function [cx,cy,sr,theta2]=p_RRP(bx,by,px,py,theta3,l2,m)

v_RRP.m function [vcx,vcy,vr,omiga2]=v_RRP(bx,by,cx,cy,vbx,vby,vpx,vpy,theta2,theta3,l2,sr,omiga3) a_RRP.m function [acx,acy,ar,alpha2]=a_RRP(bx,by,cx,cy,px,py,abx,aby,apx,apy,theta3,vr,omiga2,omiga3,alpha3)

函数中的符号说明

1 1

∠BCP < 90?，∠BC 'P > 90?，

m =1

RRP Ⅱ级杆组运动分析

函数文件（m 文件）

格式： function [ 输出参数 ] = 函数名（ 输入参数 ）

p_RPR.m function [dx,dy,sr,theta3]=p_RPR(bx,by,cx,cy,e,l3,m)

v_RPR.m function [vdx,vdy,omiga3,vr]=v_RPR(bx,by,cx,cy,dx,dy,vcx,vcy,vbx,vby,theta3) a_RPR.m function [adx,ady,alpha3,ar]=a_RPR(bx,by,cx,cy,dx,dy,acx,acy,abx,aby,vr,omiga3,theta3)

函数中的符号说明

RRP Ⅱ级杆组运动分析

实线位置，m =1 虚线位置，m = -1

函数文件（m 文件）

格式： function [ 输出参数 ] = 函数名（ 输入参数 ）

F_RRR.m function [R12x,R12y,R23x,R23y,R34x,R34y]=F_RRR(bxy,cxy,dxy,s2,s3,m2,m3,

Js2,Js3,M2,M3,F2,F3,as2,as3,alpha2,alpha3)

RRR Ⅱ级杆组力分析

R 23x

F 2R F 3x

R 23

函数文件（m 文件）

格式： function [ 输出参数 ] = 函数名（ 输入参数 ）

F_RRP.m function [R12x,R12y,R23x,R23y,R34x,R34y,lcn]=F_RRP(bxy,cxy,s2,s3,m2,m3,

Js2,Js3,M2,M3,F2,F3,theta3,as2,as3,alpha2,alph3)

RRP Ⅱ级杆组力分析

R 34

函数文件（m 文件）

格式： function [ 输出参数 ] = 函数名（ 输入参数 ）

F_RPR.m function [R12x,R12y,R23x,R23y,R35x,R35y,lcn]=F_RRP(bxy,cxy,dxy,s2,s3,

m2,m3,Js2,Js3,M2,M3,F2,F3,R34,theta3,as2,as3,alpha3)

RPR Ⅱ级杆组力分析

23

8. 作用有平衡力的构件力分析

作用有平衡力的构件力分析

函数文件（m文件）

格式：function [ 输出参数] = 函数名（输入参数）

F_Bar.m function [R01x,R01y,Mb]=F_Bar(axy,bxy,s1,m1,Js1,M1,F1,R12,as1,alpha1)

函数中的符号说明

9. 平面连杆机构运动分析算例

例1图示曲柄摇杆机构，已知l 1=150mm ，l 2=220mm ，l 3=250mm ，l 4=300mm ，曲柄以n 1=100r/min 逆时针匀速转动，分析该机构的运动。

主程序

% 曲柄摇杆机构运动分析 clc,clear

l1=150;% 曲柄长度 l2=220;% 连杆长度 l3=250;% 摇杆长度 l4=300;% 机架长度 n=100;% 曲柄转速

m=1;% RRR II 级杆组装配模式系数

omiga1=2*pi*n/60;alpha1=0;% 曲柄角速度、角加速度

ax=0;ay=0;vax=0;vay=0;aax=0;aay=0;% A 点位置、速度、加速度 dx=l4;dy=0;vdx=0;vdy=0;adx=0;ady=0;% D 点位置、速度、加速度 phi=0;% 曲柄的结构参数

theta1=0:10:360;% 曲柄转角(每隔10°计算一次) theta1=theta1*pi/180;

% 调用crank 函数，计算B 点运动参数 [bx,by]=p_crank(ax,ay,theta1,phi,l1);

[vbx,vby]=v_crank(l1,vax,vay,omiga1,theta1,phi); [abx,aby]=a_crank(l1,aax,aay,alpha1,omiga1,theta1,phi); % 调用RRR 函数，计算BC 杆和CD 杆以及C 点运动参数 [cx,cy,theta2,theta3]=p_RRR(bx,by,dx,dy,l2,l3,m);

[vcx,vcy,omiga2,omiga3]=v_RRR(vbx,vby,vdx,vdy,cx,cy,bx,by,dx,dy);

[acx,acy,alpha2,alpha3]=a_RRR(abx,aby,adx,ady,cx,cy,bx,by,dx,dy,omiga2,omiga3); % 绘制运动线图 theta1=theta1*180/pi; figure(1) subplot(3,1,1);

plot(theta1,cx,'-',theta1,cy,':r'),grid on xlabel('曲柄转角( ° )');ylabel('位移(mm/s)'); legend('C 点x 方向位移','C 点y 方向位移');

曲柄摇杆机构

subplot(3,1,2);

plot(theta1,vcx,'-',theta1,vcy,':r'),grid on

xlabel('曲柄转角( ° )');ylabel('速度(mm/s)');

legend('C点x方向速度','C点y方向速度');

subplot(3,1,3);

plot(theta1,acx,'-',theta1,acy,':r'),grid on

xlabel('曲柄转角( ° )');ylabel('加速度(mm/s^2)'); legend('C点x方向加速度','C点y方向加速度'); figure(2)

subplot(2,1,1);

plot(theta1,theta3*180/pi),grid on

xlabel('曲柄转角( ° )');ylabel('摇杆角位移( ° )'); subplot(2,1,2);

plot(theta1,omiga3),grid on

xlabel('曲柄转角( ° )');ylabel('摇杆角速度(rad/s)'); subplot(2,1,3);

plot(theta1,alpha3),grid on

xlabel('曲柄转角( ° )');ylabel('摇杆角加速度(/rad/s^2)');

例2 图示曲柄滑块机构，已知l 1=150mm ，l 2=150mm ，e =25mm ，曲柄以n 1=955r/min 逆时针匀速转动，分析该机构的运动。

主程序

% 曲柄滑块机构运动分析 clc,clear

l1=50;% 曲柄长度 l2=150;% 连杆长度 e=25;% 偏距 n=955;% 曲柄转速

m=1;% RRP II 级杆组装配模式系数

omiga1=2*pi*n/60;alpha1=0;% 曲柄角速度、角加速度

ax=0;ay=e;vax=0;vay=0;aax=0;aay=0;% A 点位置、速度、及速度 px=0;py=0;% 滑块导路上一定点（选为O 点）的坐标 vpx=0;vpy=0;apx=0;apy=0;% 滑块导路上一定点的速度、加速度 theta3=0;omiga3=0;alpha3=0;% 滑块导路角位置、角速度、角加速度 phi=0;% 曲柄的结构参数

theta1=0:30:360;% 曲柄转角(每隔10°计算一次) theta1=theta1*pi/180;

% 调用crank 函数，计算B 点运动参数 [bx,by]=p_crank(ax,ay,theta1,phi,l1);

[vbx,vby]=v_crank(l1,vax,vay,omiga1,theta1,phi); [abx,aby]=a_crank(l1,aax,aay,alpha1,omiga1,theta1,phi); % 调用RRP 函数，计算BC 杆和滑块的运动参数 [cx,cy,sr,theta2]=p_RRP(bx,by,px,py,theta3,l2,m);

[vcx,vcy,vr,omiga2]=v_RRP(bx,by,cx,cy,vbx,vby,vpx,vpy,theta2,theta3,l2,sr,omiga3); [acx,acy,ar,alpha2]=a_RRP(bx,by,cx,cy,px,py,abx,aby,apx,apy,theta3,vr,omiga2,omiga3,alpha3); % 绘制运动线图

crankx=l1.*cos(theta1);cranky=e+l1.*sin(theta1); theta1=theta1*180/pi; figure(1) subplot(3,1,1); plot(theta1,cx,'r'),grid on

xlabel('曲柄转角( ° )');ylabel('滑块位移(mm/s)');

曲柄滑块机构

subplot(3,1,2);

plot(theta1,vcx,'r'),grid on

xlabel('曲柄转角( °)');ylabel('滑块速度(mm/s)'); subplot(3,1,3);

plot(theta1,acx,'r'),grid on

xlabel('曲柄转角( °)');ylabel('滑块加速度(mm/s^2)');

例3 图示摆动导杆机构，已知l 1=280mm ，h =380mm ，l CD =840mm ，曲柄1以ω1=18.012rad/逆时针匀速转动，分析该机构的运动。

主程序

% RPR II 级杆组运动分析程序检验 clc,clear

l1=280;% 曲柄长度 h=380;% 机架长度 lcd=840;% 导杆长度

e=0;% 内接移动副导路至C 点的偏心距 m=1;% RPR II 级杆组装配模式系数 phi=0;% 曲柄的结构参数

omiga1=18.021;alpha1=0;% 曲柄角速度、角加速度

ax=0;ay=h;vax=0;vay=0;aax=0;aay=0;% A 点位置、速度、加速度 cx=0;cy=0;vcx=0;vcy=0;acx=0;acy=0;% C 点位置、速度、加速度 theta1=(0:5:360)*pi/180;% 曲柄转角(每隔5°计算一次) % 调用crank 函数，计算B 点运动参数 [bx,by]=p_crank(ax,ay,theta1,phi,l1);

[vbx,vby]=v_crank(l1,vax,vay,omiga1,theta1,phi); [abx,aby]=a_crank(l1,aax,aay,alpha1,omiga1,theta1,phi); % 调用RPR 函数，计算BC 杆运动参数 [dx,dy,sr3,theta3]=p_RPR(cx,cy,bx,by,e,lcd,m);

[vdx,vdy,omiga3,vr3]=v_RPR(cx,cy,bx,by,dx,dy,vbx,vby,vcx,vcy,theta3); [adx,ady,alpha3]=a_RPR(cx,cy,cx,by,dx,dy,abx,aby,acx,acy,vr3,omiga3,theta3); % 绘制运动线图 theta1=theta1*180/pi; subplot(3,1,1);

plot(theta1,theta3*180/pi,'r'),grid on

xlabel('曲柄转角( ° )');ylabel('摇杆角位移( ° )'); subplot(3,1,2);

plot(theta1,omiga3,'r'),grid on

xlabel('曲柄转角( ° )');ylabel('摇杆角速度(rad/s)'); subplot(3,1,3);

plot(theta1,alpha3,'r'),grid on

摆动导杆机构

xlabel('曲柄转角( ° )');ylabel('摇杆角加速度(rad/s^2)');

例4 图示六杆机构，已知l AB =80mm ，l BC =260mm ，l DE =400mm ，l CD =300mm ，l EF =460mm ，H 1=90mm ，H 2=170mm ，曲柄1逆时针匀速转动，ω1=40rad/s ，分析滑块5的运动。

主程序

% 平面六杆机构运动分析程序检验 clc,clear

l1=80;l2=260;l3=300;l31=400;l4=460;% 各构件长度 h1=90;h2=170;

e=0;% 内接移动副导路至C 点的偏心距 m1=1;m2=1;% RRR 、RRP II 级杆组装配模式系数 phi1=0;phi2=0;% 曲柄和构件DE 的结构参数 omiga1=40;alpha1=0;% 曲柄角速度、角加速度

ax=0;ay=h1;vax=0;vay=0;aax=0;aay=0;% A 点位置、速度、加速度 dx=h2;dy=0;vdx=0;vdy=0;adx=0;ady=0;% D 点位置、速度、加速度

px=0;py=h1;vpx=0,vpy=0;apx=0;apy=0;% 滑块导路上一定点（取为A 点）运动参数 theta5=0;omiga5=0;alpha5=0;% 滑块导路角运动参数 theta1=(0:10:360)*pi/180;% 曲柄转角(每隔10°计算一次) % 调用crank 函数，计算B 点运动参数 [bx,by]=p_crank(ax,ay,theta1,phi1,l1);

[vbx,vby]=v_crank(l1,vax,vay,omiga1,theta1,phi1); [abx,aby]=a_crank(l1,aax,aay,alpha1,omiga1,theta1,phi1); % 调用RRR 函数，计算C 点运动参数 [cx,cy,theta2,theta3]=p_RRR(bx,by,dx,dy,l2,l3,m1);

[vcx,vcy,omiga2,omiga3]=v_RRR(vbx,vby,vdx,vdy,cx,cy,bx,by,dx,dy);

[acx,acy,alpha2,alpha3]=a_RRR(abx,aby,adx,ady,cx,cy,bx,by,dx,dy,omiga2,omiga3); % 调用crank 函数，计算E 点运动参数 [ex,ey]=p_crank(dx,dy,theta3,phi2,l31);

[vex,vey]=v_crank(l31,vdx,vdy,omiga3,theta3,phi2); [aex,aey]=a_crank(l31,adx,ady,alpha3,omiga3,theta3,phi2); % 调用RRP 函数，计算F 点运动参数 [fx,fy,sr,theta4]=p_RRP(ex,ey,px,py,theta5,l4,m2);

[vfx,vfy,vr,omiga4]=v_RRP(ex,ey,fx,fy,vex,vey,vpx,vpy,theta4,theta5,l4,sr,omiga5); [afx,afy,ar,alpha4]=a_RRP(ex,ey,fx,fy,px,py,aex,aey,apx,apy,theta4,vr,omiga4,omiga5,alpha5);

六杆机构

% 绘制运动线图

theta1=theta1*180/pi;

subplot(3,1,1);

plot(theta1,sr,'r'),grid on

xlabel('曲柄转角( ° )');ylabel('滑块5位移(mm)'); subplot(3,1,2);

plot(theta1,vr,'r'),grid on

xlabel('曲柄转角( ° )');ylabel('滑块5速度(mm/s)'); subplot(3,1,3);

plot(theta1,ar,'r'),grid on

xlabel('曲柄转角( ° )');ylabel('滑块5加速度(mm/s^2)');

10. 平面连杆机构力分析算例

图示平面六杆机构，已知l 1=0.1m ，l 3=0.7m ，l 4=0.45m ，H =0.65 m ，H 1=0.5 m ，曲柄1逆时针匀速转动，n 1=100 r/min ，构件质量m 3=20 kg ，m 4=15 kg ，m 5=62 kg ，构件3和构件4的质心分别位于各自构件长度的一半处，且构件3、4绕其质心的转动惯量分别为J s3=0.11 kg ?m 2，J s4=0.18 kg ?m 2，，构件1和构件2的质量忽略不计，该机构在工作行程时滑块受生产阻力F r =110 N 。分析滑块5的运动规律及在工作中应加在曲柄1上的平衡力矩。

主程序

% 六杆机构力分析 clc,clear

lab=0.1; % AB 杆长度 lcd=0.7; % CD 杆长度 lac=0.5; % 机架长度 lde=0.45; %DE 杆长度 phi=0;% 曲柄的结构参数

ax=0;ay=lac;vax=0;vay=0;aax=0;aay=0;% 转动副A 位置、速度、加速度 H=0.65; % 滑块5至铰链C 垂直距离

cx=0;cy=0;vcx=0;vcy=0;acx=0;acy=0;% C 点位置、速度、加速度

px=0;py=H;vpx=0;vpy=0;apx=0;apy=0;% 滑块导路上一定点P 的位置、速度、加速度 s1=[0,0];as1=[0,0];% 曲柄质心位置和加速度

theta5=0;omiga5=0;alpha5=0;% 滑块导路的角位置、角速度、角加速度 m1=1;% RPR II 级杆组装配模式系数 m2=1;% RRP II 级杆组装配模式系数 e=0;% 内接移动副导路至C 点的偏心距 n1=100; % 曲柄转速

omiga1=pi*n1/30;alpha1=0;% 曲柄角速度、角加速度 % 曲柄AB 的运动参数

theta1=0:10:360; 曲柄转角(每隔10°计算一次) theta1=theta1*pi/180;

% 调用crank 函数，计算B 点运动参数 [bx,by]=p_crank(ax,ay,theta1,phi,lab);

[vbx,vby]=v_crank(lab,vax,vay,omiga1,theta1,phi);

平面六杆机构

[abx,aby]=a_crank(lab,aax,aay,alpha1,omiga1,theta1,phi);

% 调用RPR函数，计算CD杆和D点运动参数

[dx,dy,sr3,theta3]=p_RPR(cx,cy,bx,by,e,lcd,m1);

[vdx,vdy,omiga3,vr3]=v_RPR(cx,cy,bx,by,dx,dy,vbx,vby,vcx,vcy,theta3);

[adx,ady,alpha3,ar3]=a_RPR(cx,cy,bx,by,dx,dy,abx,aby,acx,acy,vr3,omiga3,theta3);

% 调用RRP函数，计算E点运动参数

[ex,ey,sr,theta4]=p_RRP(dx,dy,px,py,theta5,lde,m2);

[vex,vey,vr,omiga4]=v_RRP(dx,dy,ex,ey,vdx,vdy,vpx,vpy,theta4,theta5,lde,sr,omiga5); [aex,aey,ar,alpha4]=a_RRP(dx,dy,ex,ey,px,py,adx,ady,apx,apy,theta5,vr,omiga4,omiga5,alpha5); % 调用crank函数，计算CD杆和DE杆质心点的运动参数

[s3x,s3y]=p_crank(dx,dy,theta3,phi,lcd/2);

[vs3x,vs3y]=v_crank(lcd/2,vcx,vcy,omiga3,theta3,phi);

[as3x,as3y]=a_crank(lcd/2,acx,acy,alpha3,omiga3,theta3,phi);

[s4x,s4y]=p_crank(dx,dy,theta4,phi,lcd/2);

[vs4x,vs4y]=v_crank(lcd/2,vdx,vdy,omiga4,theta4,phi);

[as4x,as4y]=a_crank(lcd/2,adx,ady,alpha4,omiga4,theta4,phi);

% 力参数

m1=0;m2=0;m3=20;m4=15;m5=62;% 各构件质量

Js1=0;Js2=0;Js3=0.11;Js4=0.18;Js5=0;% 各构件绕其质心的转动惯量

M1=0;M2=0;M3=0;M4=0;M5=0;% 作用于构件1、2、3、4、5的主矩

F1=[0,0];F2=[0,0];F3=[0,0];F4=[0,0];% 作用于构件1、2、3、4的主矢

fr=110;fr5=(omiga3>0)*fr;% 生产阻力

for i=1:length(theta1)

F5=[fr5(i),0];% 作用于构件5的主矢

% 力分析

% 调用RRP II级杆组力分析函数，计算运动副D、E中的约束反力

[R34x(i),R34y(i),R45x(i),R45y(i),R56x(i),R56y(i),lcn]=F_RRP([dx(i),dy(i)],[ex(i),ey(i)],...

[s4x(i),s4y(i)],[ex(i),ey(i)],m4,m5,Js4,Js5,M4,M5,F4,F5,theta5,...

[as4x(i),as4y(i)],[aex(i),aey(i)],alpha4(i),alpha5);

% 调用RPR II级杆组力分析函数，计算运动副B、C中的约束反力

[R12x(i),R12y(i),R23x(i),R23y(i),R36x(i),R36y(i),lcn(i)]=F_RPR([cx,cy],[bx(i),by(i)],...

[dx(i),dy(i)],[bx(i),by(i)],[s3x(i),s3y(i)],m2,m3,Js2,Js3,M2,M3,F2,F3,...

[R34x(i),R34y(i)],theta3(i),[abx(i),aby(i)],[as3x(i),as3y(i)],alpha3(i));

% 调用单杆力分析函数，计算运动副A中的约束反力及平衡力矩

[R16x(i),R16y(i),Mb(i)]=F_Bar([ax,ay],[bx(i),by(i)],s1,m1,Js1,M1,F1,...

[R12x(i),R12y(i)],as1,alpha1);

end

theta1=theta1*180/pi;

figure(1)

subplot(3,1,1)

plot(theta1,sr),grid on

xlabel('曲柄转角（°）');ylabel('滑块位移（m）');

subplot(3,1,2)

plot(theta1,vr),grid on

xlabel('曲柄转角（°）');ylabel('滑块速度（m/s）');

subplot(3,1,3)

plot(theta1,ar),grid on

xlabel('曲柄转角（°）');ylabel('滑块加速度（m/s^2）');

figure(2)

plot(theta1,Mb),grid on

xlabel('曲柄转角（°）');ylabel('平衡力矩（N.m）');

《MATLAB语言与应用》实验测验高成。.doc

《MATLAB语言与应用》实验测验1 专业学号姓名成绩 自动化200909150151 高成 要求：随机生成一个4×4矩阵A，生成一个4×4的魔术矩阵B，计算C=A.*B，D=A*B，E=A/B，F=A\B，G=A.^2，H=A^2。 矩阵C、D、E、F、G和H的实验结果分别为： >> A=rand(4) A = 0.2769 0.6948 0.4387 0.1869 0.0462 0.3171 0.3816 0.4898 0.0971 0.9502 0.7655 0.4456 0.8235 0.0344 0.7952 0.6463 >> B=magic(4) B = 16 2 3 13 5 11 10 8 9 7 6 12 4 14 1 5 1 >> C=A.*B C = 4.4308 1.3897 1.3162 2.4293 0.2309 3.4881 3.8156 3.9181 0.8742 6.6516 4.5931 5.3470 3.2938 0.4822 11.9280 0.6463 >> D=A*B D = 12.6011 13.8844 13.2146 14.6104 7.7173 13.1080 12.9453 8.2055 14.9772 22.2435 21.0705 18.4963 23.0896 16.6406 17.2807 21.1692

>> E=A/B Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND = 1.306145e-017. E = 1.0e+014 * 0.8053 2.4159 -2.4159 -0.8053 -0.5976 -1.7929 1.7929 0.5976 0.1930 0.5789 -0.5789 -0.1930 -1.9751 -5.9254 5.9254 1.9751 >> F=A\B F = 32.4732 0.3697 5.9411 15.7588 34.5598 -1.5535 0.4005 28.6980 -48.3016 -4.7808 -8.8210 -36.1811 22.4018 27.1553 26.4707 24.4556 >> G=A.^2 G = 0.0767 0.4828 0.1925 0.0349 0.0021 0.1006 0.1456 0.2399 0.0094 0.9029 0.5860 0.1985 0.6781 0.0012 0.6323 0.4177 >> H=A^2 H = 0.3053 0.8361 0.8711 0.7083 0.4678 0.5121 0.8228 0.6505 0.5120 1.1116 1.3455 1.1126 0.8391 1.3610 1.4971 0.9428

科学计算与MATLAB语言考试答案

1 单选(2分) 利用MATLAB求解科学计算问题的优势是（）。 得分/总分 ? A. 算法最优 ? B. 不需要编写程序 ? C. 程序执行效率高 ? D. 编程效率高 正确答案：D你没选择任何选项 2 单选(2分) 在MATLAB命令行窗口输入命令时，可使用续行符，其写法是（）。 得分/总分 ? A. 省略号（…） ? B. 分号（;） ? C. 三个小数点（…） ? D. 百分号（%） 正确答案：C你没选择任何选项 3

下列语句执行后，D的值为（）。 1.A=[1:3;4:6]; 2.D=sub2ind(size(A),[1,1],[2,3]) 得分/总分 ? A. 3 6 ? B. 2 5 ? C. 3 5 ? D. 4 5 正确答案：C你没选择任何选项 4 单选(2分) ceil(-2.1)+floor(-2.1)+fix(-2.1)的结果为（）。 得分/总分 ? A. -7 ? B. -6 ? C. -5 ? D. -9 正确答案：A你没选择任何选项 5

下列语句执行后，x的值是（）。 1.log=1:5; 2.x=log(1) 得分/总分 ? A. ? B. 1 ? C. 数学常数e ? D. 报错 正确答案：B你没选择任何选项 6 单选(2分) 下列语句执行后，c的值是（）。 1.ch=['abcdef';'123456']; 2.c=char(ch(2,4)-1) 得分/总分 ? A. '4' ? B. 4 ? C. '3' ? D. 3

7 单选(2分) 产生和A同样大小的全0矩阵的函数是（）。 得分/总分 ? A. zero(size(A)) ? B. zeros(size(A)) ? C. size(zero(A)) ? D. size(zeros(A)) 正确答案：B你没选择任何选项 8 单选(2分) 语句x=speye(5)==eye(5)执行后，则下列说法中正确的是（）。 得分/总分 ? A. x是5阶全1矩阵，且采用稀疏存储方式 ? B. x是5阶全1矩阵，且采用完全存储方式 ? C. x是5阶单位矩阵，且采用稀疏存储方式 ? D. x是5阶单位矩阵，且采用完全存储方式

MATLAB精通科学计算_偏微分方程求解

一、Maple V 系统 Maple V是由Waterloo大学开发的数学系统软件，它不但具有精确的数值处理功能，而且具有无以伦比的符号计算功能。Maple V的符号计算能力还是MathCAD和MATLAB等软件的符号处理的核心。Maple提供了2000余种数学函数，涉及范围包括：普通数学、高等数学、线性代数、数论、离散数学、图形学。它还提供了一套内置的编程语言，用户可以开发自己的应用程序，而且Maple自身的2000多种函数，基本上是用此语言开发的。 Maple采用字符行输入方式，输入时需要按照规定的格式输入，虽然与一般常见的数学格式不同，但灵活方便，也很容易理解。输出则可以选择字符方式和图形方式，产生的图形结果可以很方便地剪贴到Windows应用程序内。 二、MATLAB 系统 MATLAB原是矩阵实验室（Matrix Laboratory）在70年代用来提供Linpack和Eispac k软件包的接口程序，采用C语言编写。从80年代出现3.0的DOS版本，逐渐成为科技计算、视图交互系统和程序语言。MATLAB可以运行在十几个操作平台上，比较常见的有基于W indows 9X/NT、OS/2、Macintosh、Sun、Unix、Linux等平台的系统。 MATLAB程序主要由主程序和各种工具包组成，其中主程序包含数百个内部核心函数，工具包则包括复杂系统仿真、信号处理工具包、系统识别工具包、优化工具包、神经网络工具包、控制系统工具包、μ分析和综合工具包、样条工具包、符号数学工具包、图像处理工具包、统计工具包等。而且5.x版本还包含一套几十个的PDF文件，从MATLAB的使用入门到其他专题应用均有详细的介绍。 MATLAB是数值计算的先锋，它以矩阵作为基本数据单位，在应用线性代数、数理统计、自动控制、数字信号处理、动态系统仿真方面已经成为首选工具，同时也是科研工作人员和大学生、研究生进行科学研究的得力工具。MATLAB在输入方面也很方便，可以使用内部的E ditor或者其他任何字符处理器，同时它还可以与Word6.0/7.0结合在一起，在Word的页面里直接调用MATLAB的大部分功能，使Word具有特殊的计算能力。 三、MathCAD 系统 MathCAD是美国Mathsoft公司推出的一个交互式的数学系统软件。从早期的DOS下的1. 0和Windows下的4.0版本，到今日的8.0版本，功能也从简单的数值计算，直至引用Map le强大的符号计算能力，使得它发生了一个质的飞跃。 MathCAD是集文本编辑、数学计算、程序编辑和仿真于一体的软件。MathCAD7.0 Profe ssional（专业版）运行在Win9X/NT下，它的主要特点是输入格式与人们习惯的数学书写格式很近似，采用WYSWYG（所见所得）界面，特别适合一般无须进行复杂编程或要求比较特殊的计算。MathCAD 7.0 Professional 还带有一个程序编辑器，对于一般比较短小，或者要求计算速度比较低时，采用它也是可以的。这个程序编辑器的优点是语法特别简单。 MathCAD可以看作是一个功能强大的计算器，没有很复杂的规则；同时它也可以和Wor d、Lotus、WPS2000等字处理软件很好地配合使用，可以把它当作一个出色的全屏幕数学公式编辑器。 四、Mathematica 系统 Mathematica是由美国物理学家Stephen Wolfram领导的Wolfram Research开发的数学系统软件。它拥有强大的数值计算和符号计算能力，在这一方面与Maple类似，但它的符

科学计算与MATLAB语言(第四课)

第四讲绘图功能

作为一个功能强大的工具软件，Matlab 具有很强的图形处理功能，提供了大量的二维、三维图形函数。由于系统采用面向对象的技术和丰富的矩阵运算，所以在图形处理方面即常方便又高效。

4.1 二维图形 一、plot函数 函数格式：plot(x,y)其中x和y为坐标向量函数功能：以向量x、y为轴，绘制曲线。【例1】在区间0≤X≤2 内，绘制正弦曲线Y=SIN（X），其程序为： x=0:pi/100:2*pi; y=sin(x); plot(x,y)

一、plot函数 【例2】同时绘制正、余弦两条曲线Y1=SIN（X）和Y2=COS（X），其程序为： x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,x,y2) plot函数还可以为plot(x,y1,x,y2，x,y3，…)形式，其功能是以公共向量x为X轴，分别以y1，y2，y3，…为Y轴，在同一幅图内绘制出多条曲线。

一、plot函数 （一）线型与颜色 格式：plot(x,y1,’cs’,...) 其中c表示颜色，s表示线型。 【例3】用不同线型和颜色重新绘制例4.2图形，其程序为：x=0:pi/100:2*pi; y1=sin(x); y2=cos(x); plot(x,y1,'go',x,y2,'b-.') 其中参数'go'和'b-.'表示图形的颜色和线型。g表示绿色，o表示图形线型为圆圈；b表示蓝色，-.表示图形线型为点划线。

一、plot函数 （二）图形标记 在绘制图形的同时，可以对图形加上一些说明，如图形名称、图形某一部分的含义、坐标说明等，将这些操作称为添加图形标记。 title(‘加图形标题'); xlabel('加X轴标记'); ylabel('加Y轴标记'); text(X,Y,'添加文本');

中南大学材料学院科学计算与MATLAB考试题库

练习题 1.求函数在指定点的数值导数 x=sym('x'); >> y=[x x.^2 x.^3;1 2*x 3*x.^2;0 2 6*x]; >> x=1; >> eval(diff(y)) ans = 1 2 3 0 2 6 0 0 6 >> x=2; >> eval(diff(y)) ans = 1 4 12 0 2 12 0 0 6 >> x=3; >> eval(diff(y)) ans = 1 6 27 0 2 18 0 0 6 2.求下列函数导数 （1） x=sym('x'); >> y=x^10+10^x+(log(10))/log(x); >> diff(y) ans = 10*x^9+10^x*log(10)-2592480341699211/1125899906842624/log(x)^2/x （2） x=sym('x');

>> y=log(1+x); >> x=1; >> eval(diff(y,2)) %在x=1的条件下对y表达式求两次导数后导函数的值 ans = -0.2500 3.用数值方法求下列积分 首先先讲一下trapz的用法，如下题 t=0:0.001:1; >> y=t; >> trapz(t,y) ans = 0.5000 （1） >> x=1:0.01:5; >> y=(x.^2).*sqrt(2*x.^2+3); >> trapz(x,y) ans = 232.8066 (2) x=pi/4:0.01:pi/3; >> y=x./(sin(x).^2); >> trapz(x,y) ans = 0.3810 第三题拟合曲线题 x=[0:0.1:1]; >> y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; >> a=polyfit(x,y,2); >> x=[0.05:0.2:1.05]; >> y=a(3)+a(2)*x+a(1)*x.^2 %注意x要在y前先赋值，不然y不会运行为最新的x对呀的y值 y =

哈工大_Matlab__2013年春季学期《MATLAB语言及应用》试题答案

2013年春季学期 《MATLAB语言及应用》课程试卷 姓名： 学号： 学院： 专业： 必答题 1.常用的matlab界面由哪些窗口组成，各有什么主要作用？(4分) (1)菜单和工具栏功能：【File】菜单主要用于对文件的处理。【Edit】菜单主 要用于复制、粘贴等操作，与一般Windows程序的类似，在此不作详细介绍。【Debug】菜单用于调试程序。【Desktop】菜单用于设置主窗口中需要打开的窗口。【Window】菜单列出当前所有打开的窗口。【Help】菜单用于选择打开不同的帮助系统。 (2)命令窗口功能：用于输入命令并显示除图形以外的所有执行结果 (3)历史命令窗口功能：主要用于记录所有执行过的命令 (4)当前工作目录窗口功能：对文件和目录进行操作 (5)工作空间窗口功能：查看、载入和保存变量 2.如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么 区别？(2分) 方法一：在MATLAB命令窗口中输入editpath或pathtool命令或通过【File】/|【SetPath】菜单，进入“设置搜索路径”对话框，通过该对话框编辑搜索路径。 方法二：在命令窗口执行“path(path,…D:\Study ?)”，然后通过“设置搜索路径”对话查看“D:\Study”是否在搜索路径中。 方法三：在命令窗口执行“addpath D:\Study- end”，将新的目录加到整个搜索路径的末尾。如果将end改为begin，可以将新的目录加到整个搜索路径的开始。 区别：当前文件目录是正在运行的文件的目录，显示文件及文件夹的详细信息，且只有将文件设置为当前目录才能直接调用。搜索路径中的文件可以来自多个 不同目录，在调用时不用将其都设置为当前目录，为同时调用多个文件提供 方便。 3.有几种建立矩阵的方法？各有什么优点？(4分) 1.在命令窗口中直接输入优点：适合输入不规则和较小的矩阵 2.通过语句和函数生成矩阵优点：适合输入规则且较大的矩阵 3.通过M文件来建立矩阵，从外部数据文件中导入矩阵 优点：方便创建和导入大型矩阵 4.说明break语句、continue语句和return语句的用法。(3分)

《Matlab与科学计算》作业 2010010099

《Matlab与科学计算》作业 第一章MATLAB环境 1、MATLAB通用操作界面窗口包括哪些？命令窗口、历史命令窗口、当前目录窗口、工作空间窗口各有哪些功能？ 答：MATLAB通用操作界面窗口包括：命令窗口、历史命令窗口、当前目录浏览器窗口、工作空间窗口、变量编辑器窗口、M文件编辑/调试器窗口、程序性能剖析窗口、MATLAB帮助。 命令窗口是MATLAB命令操作的最主要窗口，可以把命令窗口当做高级的“草稿纸”。在命令窗口中可以输入各种MATLAB的命令、函数和表达式，并显示除图形外的所有运算结果。 历史命令窗口用来记录并显示已经运行过的命令、函数和表达式，并允许用户对它们进行选择、复制和重运行，用户可以方便地输入和修改命令，选择多行命令以产生M文件。 当前目录窗口用来设置当前目录，可以随时显示当前目录下的M、MKL等文件的信息，扬文件类型、文件名、最后个修改时间和文件的说明信息等，并可以复制、编辑和运行M文件及装载MAT数据文件。 工作空间窗口用来显示所有MATLAB工作空间中的变量名、数据结构、类型、大小和字节数。 2、熟悉课本中表格1.4、1.5、1.6、1.7、1.8的内容。 3、如何生成数据文件？如何把数据文件中的相关内容输入到工作空间中，用实例进行操作。 生成数据文件：

把数据文件中的相关内容输入到工作空间中： 结果： 4、在工作空间中可以通过哪些命令管理变量，写出每种语法的具体操作过程。答：（1）把工作空间中的数据存放到MAT数据文件。 语法：save filename 变量1 变量2 ……参数。 （2）从数据文件中取出变量存放到工作空间。 语法：load filename 变量1 变量2 ……。

MATLAB语言及应用课程教学大纲

《MATLAB语言及应用》课程教学大纲 课程编号：21311105 总学时数：32 总学分数：2 课程性质：专业必修课 适用专业：电气工程及其自动化 一、课程的任务和基本要求： 控制算法是集中现代控制系统CAD 技术的本质反映，对于控制算法系统深入地学习，在扩展计算机技术在控制理论中的应用和发展，同时培养学生运用计算机技术进行思维和开发的能力。控制系统理论、计算方法与计算机技术的结合是当代控制理论发展的标志，因此在以MATLAB 为代表的软件平台上，对控制系统进行分析、设计与仿真将成为控制工程领域工程师必须熟练掌握的重要知识和技能。因此，深入透彻地分析和理解控制算法的思想和构造就必须系统学习典型控制系统应用软件的基本原理和控制算法，将成为本课程的目的和任务。 通过本课程的学习，要求学生掌握对于控制系统的分析和综合设计的方法和基本技巧，而控制算法在控制系统CAD 技术中占有相当大的比重，本课程要求较熟练掌握控制算法的基本思想；MATLAB 是一种解释性编程语言，因此，要求熟练掌握MATLAB 的基本编程手段和模块化编程方法，消化和理解控制语言描述的图形界面的设计过程。 二、基本内容和要求： 1.自动控制系统与仿真基础知识 （1）自动控制系统基本概念 （2）自动控制系统分类 （3）控制系统仿真基本概念 （4）MA TLAB与控制系统仿真 （5）MA TLAB 7中控制相关的工具箱 要求：了解自动控制系统与仿真的基础知识，包括自动控制系统的基本概念、分类，以及控制系统仿真的基本概念和Matlab工具。 2.MA TLAB计算及仿真基础 （1）MA TLAB概述 （2）MA TLAB桌面操作环境 （3）MA TLAB数值计算 （4）关系运算和逻辑运算 （5）符号运算 （6）复数和复变函数运算 （7）MA TLAB常用绘图命令 （8）MA TLAB程序设计 要求：了解MA TLAB计算及仿真基础，包括MATLAB的安装、界面及其数值计算、函数运算、程序设计及其绘图命令。 3. Simulink仿真基础 （1）Simulink仿真概述

Matlab入门基础知识整理

MATLAB基础知识 MATLAB的主要功能 1．数值计算和符号计算功能 例如，求解线性方程组 在MATLAB命令窗口输入命令： a=[2,3,-1;3,-5,3;6,3,-8]; b=[7;8;9]; x=inv(a)*b 也可以通过符号计算来解此方程 syms x1 x2 x3 [x1,x2,x3]=solve(2*x1+3*x2-x3-7,3*x1-5*x2+3*x3-8,6*x1+3*x2-8*x3-9) 2．绘图功能 例如，分别绘制函数y=300sinx/x和y=x2的曲线 x=-20:0.1:20; plot(x,300*sin(x)./x,':',x,x.^2); 3．程序设计语言功能 MATLAB是解释性语言，程序执行速度较慢，而且不能脱离MATLAB环境而独立运行。MathWorks公司提供了将MATLAB源程序编译为独立于MATLAB集成环境运行的EXE文件以及将MATLAB程序转化为C语言程序的编译器。4．扩展功能 MATLAB包含两部分内容：基本部分和各种可选的工具箱。 基本部分构成了MATLAB的核心内容，也是使用和构造工具箱的基础。 工具箱扩展了MATLAB的功能。 功能性工具箱 学科性工具箱 MATLAB的集成开发环境 命令窗口（Command Window）用于输入命令并显示除图形以外的所有执行结果。一般来说，一个命令行输入一条命令，命令行以回车或分号结束 但一个命令行也可以输入若干条命令，各命令之间以逗号或分号隔开 如果一个命令行很长，一个物理行之内写不下，可以在第1个物理行之后加上3个小黑点并按下回车键，然后接着下一个物理行继续写命令的其他部分。3个小黑点称为续行符，即把下面的物理行看作该行的逻辑继续。例如： z=1+1/(1*2)+1/(1*2*3)+1/(1*2*3*4)+ … 1/(1*2*3*4*5) 工作空间（Workspace）是MATLAB用于存储各种变量的内存空间。 当前目录（Current Directory）是指MATLAB运行时的工作目录，只有在当前目录下的文件、函数才可以被运行或调用。 如果没有特殊指明，数据文件也将存放在当前目录下。 命令历史窗口 自动保留自安装起所有用过的命令的历史记录，并且还标明了使用时间，从而方便用户查询。而且，通过双击命令可进行历史命令的再运行。 MATLAB的帮助功能 进入MATLAB帮助界面可以通过以下方法。 ●单击MATLAB主窗口工具栏中的Help按钮。

《MATLAB与科学计算》期末论文

盐城师范学院《MATLAB与科学计算》期末论文 2016－2017学年度第一学期 用MATLAB解决解析几何的图形问题 学生姓名吴梦成 学院数学与统计学院 专业信息与计算科学 班级数15(5)信计 学号 15213542

用MATLAB 解决解析几何的图形问题 摘 要 将 MATLAB 的图形和动画功能都用于解析几何教学，可使教学形象生动。以图形问题为例，详细给出了实例的程序编写和动画实现过程 。在解析几何教学中有一定的应用价值。 【关键词】： MATLAB ； 解析几何 ；图形 ； 动 画；编程 1 引 言 在解析几何的教学中，使用传统的教学方法。许多曲线及曲面的形成过程与变换过程只通过传统的教师讲授静态图示就很难形象生动地表示出来 。在解析几何教学中使用MATLAB 软件辅助教学，不仅可以很容易绘制出复杂的立体图形。把曲线、曲面的形成和变化过程准确地模拟出来 ，而且还能够对它们进行翻转 、旋转 ，甚 至还能够轻而易举地实现图形的动画效果!这对提高教学效率和培养学生的空间想象能力可起到事半功倍的效果。下面结合实例从几个方面说明MATLAB 在解析几何画图方面的应用。 2 利用 MATLAB 绘制三维曲线 在空间解析几何中，各种曲线和曲面方程的建立都离不开图形 ，而空间曲线和曲面图形既难画又费时。借助MATLAB 的绘图功能 ，可以快捷 、 准确地绘出图形，使教学变得形象 、生动 。有利于学生观察三维空间图形的形状 ， 掌握图形的性质 。 一 般地 ，MATLAB 可用plot3，ezplot3，comet3等函数来各种三维曲线 。 例如画螺旋曲线的图形，其参数方程设为 ：t at cos x =，t b sin t y -=，ct =z 。使用 plot3语句画螺旋曲线图形的方法如下( 设a =2 ,b=4,c=3)： );*3),sin(*.*4),cos(*.*2(3;*10:50/:0t t t t t plot pi pi t -= MATLAB 用两条简单的语句就可以画出螺旋 曲线(图1)，但上述方法是静态的 ，为了体

科学计算与matlab1.5

单元测验已完成成绩：分 1 【单选题】 MATLAB一词来自（）的缩写。 A、 Mathematica Laboratory B、 Matrix Laboratory C、 MathWorks Lab D、 Matrices Lab 我的答案：B得分：分 2 【单选题】 下列选项中能反应MATLAB特点的是（）。 A、 算法最优 B、 不需要写程序 C、 程序执行效率高 D、 编程效率高

我的答案：D得分：分 单元测验已完成成绩：分 1 【单选题】 当在命令行窗口执行命令时，如果不想立即在命令行窗口中输出结果，可以在命令后加上（）。 A、 冒号（:） B、 逗号（,） C、 分号（;） D、 百分号（%） 我的答案：C得分：分 2 【单选题】 fix(264/100)+mod(264,10)*10的值是（）。 A、 86 B、 62 C、 423 D、

42 我的答案：D得分：分 3 【单选题】 在命令行窗口输入下列命令后，x的值是（）。 >> clear >> x=i*j A、 不确定 B、 -1 C、 1 D、 i*j 我的答案：B得分：分 4 【单选题】 使用语句x=linspace(0,pi,6)生成的是（）个元素的向量。 A、 8 B、 7 C、 6

D、 5 我的答案：C得分：分 5 【单选题】 ceil的结果为（）。 A、 -2 B、 -3 C、 1 D、 2 我的答案：A得分：分 6 【单选题】 eval('sqrt(4)+2')的值是（）。 A、 sqrt(4)+2 B、 4 C、 2 D、

2+2 我的答案：B得分：分 7 【单选题】 已知a为3×5矩阵，则执行完a(:,[2,4])=[]后（）。 A、 a变成行向量 B、 a变为3行2列 C、 a变为3行3列 D、 a变为2行3列 我的答案：C得分：分 8 【单选题】 在命令行窗口输入以下命令 >> A=[1:3;4:6]; >> D=sub2ind(size(A),[1,1],[2,3]) D的值为（）。 A、 3 6 B、 2 5 C、 4 5

matlab特点及应用领域

MATLAB就是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模与仿真等诸多强大功能集成在一个易于使用的视窗环境中,为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案,并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式,代表了当今国际科学计算软件的先进水平。 MATLAB与Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数与数据、实现算法、创建用户界面、连 matlab开发工作界面 接其她编程语言的程序等,主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的基本数据单位就是矩阵,它的指令表达式与数学、工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。

在新的版本中也加入了对C,FORTRAN,C++ ,JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用,此外许多的MATLAB爱好者都编写了一些经典的程序,用户可以直接进行下载就可以用。 一种语言之所以能如此迅速地普及,显示出如此旺盛的生命力,就是由于它有着不同于其她语言的特点。正如同FORTRAN与C等高级语言使人们摆脱了需要直接对计算机硬件资源进行操作一样,被称作为第四代计算机语言的MATLAB,利用其丰富的函数资源,使编程人员从繁琐的程序代码中解放出来。MATLAB的最突出的特点就就是简洁。MATLAB用更直观的、符合人们思维习惯的代码,代替了C与FORTRAN语言的冗长代码。MATLAB给用户带来的就是最直观、最简洁的程序开发环境。以下简单介绍一下MATLAB的主要特点。 ①语言简洁紧凑,使用方便灵活,库函数极其丰富。MATLAB程序书写形式自由,利用其丰富的库函数避开繁杂的子程序编程任务,压缩了一切不必要的编程工作。由于库函数都由本领域的专家编写,用户不必担心函数的可靠性。可以说,用MATLAB进行科技开发就是站在专家的肩膀上。 具有FORTRAN与C等高级计算机语言知识的读者可能已经注意到,如果用FORTRAN或C语言去编写程序,尤其当涉及矩阵运算与画图时,编程会很麻烦。例如,如果用户想求解一个线性代数方程,就得编写一个程序块读入数据,然后再使用一种求解线性方程的算法(例如追赶法)编写一个程序块来求解方程,最后再输出计算结果。在求解过程中,

Matlab与科学计算样题(加主观题答案)

Matlab 与科学计算考试样题（客观题） 1 下面的MATLAB 语句中正确的有： a) 2a ＝pi 。 b) record_1=3+4i c) a=2.0, d) c=1+6j 2. 已知水的黏度随温度的变化公式如下，其中a=0.03368，b=0.000221，计算温度t 为20，30，40度时的粘度分别是： 2 1at bt μμ=++0μ为0℃水的黏度，值为31.78510-?；a 、b 为常数，分别为0.03368、0.000221。 3. 请补充语句以画出如图所示的图形： [x,y]=meshgrid(-2:0.1:2, -2:0.1:2)。 Z=x.*exp(-x.^2-y.^2)。 。 a) Plot3(x,y,Z) b) plot3(x,y,Z) c) mesh(x,y,Z) d) plot3(x,y,z) 2 a) 0.4900 1.2501 0.8560 b) 0.8560 1.2501 0.4900 c) -0.6341 3.8189 -3.7749 d) 3.8189 -3.7749 2.8533 解释说明：

>> x=0.5:0.5:3.0。 >> y=[1.75,2.45,3.81,4.80,8.00,8.60]。 >> a=polyfit(x,y,2) a = 0.4900 1.2501 0.8560 >> x1=[0.5:0.25:3.0]。 >> y1=a(1)*x1.^2+a(2)*x1+a(3) >> plot(x,y,'*') >> hold on >> plot(x1,y1,'--r') 5. 求方程在 x=0.5附近的根. 21 x x += a) 0.6180 b) -1.1719e-25 c) －1 d) -1.6180 6. 用Newton-Cotes方法计算如下积分 1 5 x? （a）133.6625 (b)23.8600 (c) 87.9027 (d) -1.6180 7. y=ln(1+x)，求x=1时y" a) -0.25 b) 0.5 c) -0.6137 d) -1.6137 8.某公司用3台轧机来生产规格相同的铝合金薄板。取样测量薄板的 厚度，精确至‰厘M。得结果如下： 轧机1：0.236 0.238 0.248 0.245 0.243 轧机2：0.257 0.253 0.255 0.254 0.261 轧机3：0.258 0.264 0.259 0.267 0.262 计算方差分析结果，并判定各台轧机所生产的薄板的厚度有无显著的差异？ a) p＝1.3431e-005，没有显著差异。

Matlab语言及应用论文

一、MATLAB简介 MATLAB 是美国MathWorks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。 MATLAB是矩阵实验室（Matrix Laboratory）的简称，和Mathematica、Maple 并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。 MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且mathwork也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++ ，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。 二、应用 MATLAB 产品族可以用来进行以下各种工作： （1）数值分析； （2）数值和符号计算； （3）工程与科学绘图； （4）控制系统的设计与仿真； （5）数字图像处理技术； （6）数字信号处理技术； （7）通讯系统设计与仿真； （8）财务与金融工程。 MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用 MATLAB 函数集）扩展了 MATLAB 环境，以解决这些应用领域内特

《科学计算与MATLAB》期末大作业

杭州电子科技大学信息工程学院《科学计算与MATLAB》期末大作业

给出程序、图、作业分析，程序需加注释。 1. 试编写名为fun.m 的MATLAB 函数，用以计算下述的值： ?? ? ??-<->=t t n t t t n t f 的)4/sin()(si 对所有)4/sin(其他情况)sin(的)4/sin()(si 对所有)4/sin()(ππππ 绘制t 关于函数f(t)的图形，其中t 的取值范围为ππ66≤≤-t ，间距为10/π。 function y=fun()%定义函数 % t=-6*pi:pi/10:6*pi; %定义变量范围 y = (sin(pi/4)).*(sin(t)>sin(pi/4))+(sin(-pi/4)).*(sin(t)=sin(-pi/4)));%函数表示 plot(t,y); %画图 end

2.解以下线性方程组 ??? ??=+=++=--3 530 42231 321321x x x x x x x x A=[2 -1 -1;1 1 4;3 0 5];%输入矩阵 B=[2;0;3]; %输入矩阵 X = A\B %计算结果 3.已知矩阵? ? ??? ???? ???=44434241 3433323124232221 14131211A 求： (1)A(2:3,2:3) (2)A(:,1:2) (3)A(2:3,[1,3]) (4)[A,[ones(2,2)；eye(2)]]

A=[11 12 13 14;21 22 23 24;31 32 33 34;41 42 43 44];%输入矩阵A(2:3,2:3) %输出矩阵 A(:,1:2) %输出矩阵 A(2:3,[1,3]) %输出矩阵 [A,[ones(2,2);eye(2)]] %输出矩阵

科学计算与MATLAB 1.5

单元测验已完成成绩：100.0分 1 【单选题】 MATLAB一词来自（）的缩写。 ?A、 Mathematica Laboratory ?B、 Matrix Laboratory ?C、 MathWorks Lab ?D、 Matrices Lab 我的答案：B得分：50.0分 2 【单选题】 下列选项中能反应MATLAB特点的是（）。?A、 算法最优 ?B、 不需要写程序 ?C、 程序执行效率高 ?D、 编程效率高

我的答案：D得分：50.0分 单元测验已完成成绩：96.4分 1 【单选题】 当在命令行窗口执行命令时，如果不想立即在命令行窗口中输出结果，可以在命令后加上（）。 ?A、 冒号（:） ?B、 逗号（,） ?C、 分号（;） ?D、 百分号（%） 我的答案：C得分：7.1分 2 【单选题】 fix(264/100)+mod(264,10)*10的值是（）。 ?A、 86 ?B、 62 ?C、 423

?D、 42 我的答案：D得分：7.1分 3 【单选题】 在命令行窗口输入下列命令后，x的值是（）。 >> clear >> x=i*j ?A、 不确定 ?B、 -1 ?C、 1 ?D、 i*j 我的答案：B得分：7.1分 4 【单选题】 使用语句x=linspace(0,pi,6)生成的是（）个元素的向量。?A、 8 ?B、 7

?C、 6 ?D、 5 我的答案：C得分：7.1分 5 【单选题】 ceil(-2.1)的结果为（）。?A、 -2 ?B、 -3 ?C、 1 ?D、 2 我的答案：A得分：7.1分 6 【单选题】 eval('sqrt(4)+2')的值是（）。?A、 sqrt(4)+2 ?B、

(完整版)《MATLAB语言及其应用》教案

MATLAB语言及其应用 教 案 任课教师：罗靖宇 任课班级：09通信（1）（2） 教材：MATLAB程序设计与应用（第二版） 作者：刘卫国主编出版社：高等教育出版社

内容： 书籍简介：本书以MATLAB 7.0版为蓝本介绍MATLAB功能与应用强调理论和实践相结合贴近读者需要注重讲清有关数学方法和算法原理前提下介绍MATLAB功能；注重和有关学科领域结合，突出应用书中有许多应用实例些实例既是对MATLAB重点和难点诠释又可以更好地帮助读者应用MATLAB来解决实际问题具有很强代表性。 全书分为基础篇应用篇和实验篇基础篇包括MATLAB系统环境MATLAB数据及其运算MATLAB矩阵分析与处理 MATLAB程序设计MATLAB图MATLAB数值计算MATLAB符号计算应用篇包括MATLAB图形用户界面设计MATLAB Notebook使用MATLAB Simulink 仿真软件MATLAB外部程序接口技术MATLAB应用实验篇和教学内容相配合包括15实验以帮助读者更好地上机操作。本书可作为高等学校理工科各专业大学生或研究生学习教材也可供广大科技工作者参考。 第1章MATLAB系统环境 1.1 MA TLAB概貌 1.1.1 MA TLAB 发展 1.1.2 MA TLAB 主要功能 1.1.3 MA TLAB功能演示 1.2 MA TLAB环境准备 1.2.1 MA TLAB 安装 1.2.2 MA TLAB 启动与退出 1.3 MA TLAB操作界面 1.3.1 主窗口 1.3.2 命令窗口 1.3.3 工作空间窗口 1.3.4 当前目录窗口和搜索路径 1.3.5 命令历史记录窗口 1.3.6 Stalt菜单 1.4 MA TLAB帮助系统 1.4.1 帮助命令 1.4.2 帮助窗口 1.4.3 演示系统 第2章MATLAB数据及其运算 2.1 MA TLAB数据特点 2.2 变量及其操作 2.2.1 变量与赋值 2.2.2 变量管理 2.2.3 数据输出格式 2.3 MA TLAB矩阵表示 2.3.1 矩阵建立 2.3.2 冒号表达式 2.3.3 矩阵拆分 2.4 MA TLAB数据运算 2.4.1 算术运算 2.4.2 关系运算 2.4.3 逻辑运算 2.5 字符串 2.6 结构数据和单元数据 2.6.1 结构数据 2.6.2 单元数据 第3章MATLAB矩阵分析与处理 3.1 特殊矩阵 3.1.1 通用特殊矩阵 3.1.2 用专门学科特殊矩阵 3.2 矩阵结构变换 3.2.1 对角阵与三角阵 3.2.2 矩阵转置与旋转 3.3 矩阵求逆与线性方程组求解 3.3.1 矩阵逆与伪逆 3.3.2 用矩阵求逆方法求解线性方程组 3.4 矩阵求值 3.4.1 方阵行列式值 3.4.2 矩阵秩与迹 3.4.3 向量和矩阵范数 3.4.4 矩阵条件数 3.5 矩阵特征值与特征向量 3.6 矩阵超越函数 第4章MATLAB程序设计 4.1 M文件 4.1.1 M文件分类 4.1.2 M文件建立与打开 4.2 程序控制结构 顺序结构、选择结构、循环结构 4.3 函数文件 4.4 程序调试 第5章MATLAB绘图 5.1 二维图形

matlab程序设计与应用(第二版)第三章部分课后答案

第三章1. (1)A=eye(3) (2)A=100+100*rand(5,6) (3)A=1+sqrt(0.2)*randn(10,50) (4)B=ones(size(A)) (5)A+30*eye(size(A)) (6)B=diag(diag(A)) 2. B=rot90(A) C=rot90(A,-1) 3. B=inv(A) ;A的逆矩阵 C=det(A) ;A的行列式的值 D=A*B E=B*A D=E 因此A与A-1是互逆的。 4. A=[4 2 -1;3 -1 2;12 3 0]; b=[2;10;8]; x=inv(A)*b x = -6.0000 26.6667 27.3333 5. (1) diag(A) ；主对角线元素 ans = 1 1 5 9 triu(A) ；上三角阵

ans = 1 -1 2 3 0 1 -4 2 0 0 5 2 0 0 0 9 tril(A) ；下三角阵 ans = 1 0 0 0 5 1 0 0 3 0 5 0 11 15 0 9 rank(A) ；秩 ans = 4 norm(A) ；范数 ans = 21.3005 cond(A) ；条件数 ans = 11.1739 trace(A) ；迹 ans = 16 （2）略 6. A=[1 1 0.5;1 1 0.25;0.5 0.25 2] A = 1.0000 1.0000 0.5000 1.0000 1.0000 0.2500 0.5000 0.2500 2.0000

[V,D]=eig(A) V = 0.7212 0.4443 0.5315 -0.6863 0.5621 0.4615 -0.0937 -0.6976 0.7103 D = -0.0166 0 0 0 1.4801 0 0 0 2.5365

MATLAB-语言及其用-实验(答案)

M A T L A B-语言及其用- 实验(答案) -CAL-FENGHAI.-(YICAI)-Company One1

《MATLAB 语言及其用》 实验指导书 目录

实验一 Matlab 使用方法和程序设计........................ 实验二控制系统的模型及其转换............................. 实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink......................... 实验一Matlab使用方法和程序设计一、实验目的 1、掌握Matlab软件使用的基本方法； 2、熟悉Matlab的数据表示、基本运算和程序控制语句 3、熟悉Matlab绘图命令及基本绘图控制 4、熟悉Matlab程序设计的基本方法 二、实验内容： 1、帮助命令 使用help命令，查找 sqrt（开方）函数的使用方法； 在 CommandWindowL里输入help，接在在search里输入sqr即可。sqrt Square root Syntax B = sqrt(X) Description

B = sqrt(X) returns the square root of each element of the array X. For the elements of X that are negative or complex, sqrt(X) produces complex results. Tips See sqrtm for the matrix square root. Examples sqrt((-2:2)') ans = 0 + 1.4142i 0 + 1.0000i 1.0000 1.4142 See Also nthroot | realsqrt | sqrtm 2、矩阵运算 （1）矩阵的乘法 已知A=[1 2;3 4]; B=[5 5;7 8]; 求A^2*B A=[1 2;3 4];B=[5 5;7 8]; C=A^2*B >> format compact C = 105 115 229 251 （2）矩阵除法 已知 A=[1 2 3;4 5 6;7 8 9]; B=[1 0 0;0 2 0;0 0 3]; A\B,A/B A=[1 2 3;4 5 6;7 8 9];