MATLAB程序设计与应用第二版刘卫国实验部分.docx
- 文档编号:6870996
- 上传时间:2023-01-11
- 格式:DOCX
- 页数:15
- 大小:17.30KB
MATLAB程序设计与应用第二版刘卫国实验部分.docx
《MATLAB程序设计与应用第二版刘卫国实验部分.docx》由会员分享,可在线阅读,更多相关《MATLAB程序设计与应用第二版刘卫国实验部分.docx(15页珍藏版)》请在冰豆网上搜索。
MATLAB程序设计与应用第二版刘卫国实验部分
实验一
1.
(1)
z1=2*sin(85*pi/180)/(1+exp
(2))
(2)
>>x=[2,1+2i;-0.45,5];
>>z2=(1/2)*log(x+sqrt(1+x*x))
(3)
>>a=-3.0:
0.1:
3.0;
>>z3=((exp(0.3*a)-exp(-0.3*a))/2).*sin(a+0.3)+log((0.3+a)/2)
(4)
t=0:
0.5:
2.5;
>>z4=(t.*t).*(t>=0&t<1)+(t.*t-1).*(t>=1&t<2)+(t.*t-2.*t+1).*(t>=2&t<3)
2.
(1)
>>A=[12,34,-4;34,7,87;3,65,7];
>>B=[1,3,-1;2,0,3;3,-2,7];
>>A+6*B
(2)
>>A*B
(3)
>>A^3
(4)
>>A/B
(5)
>>[A,B]
3.
(1)
>>A=[1,2,3,4,5;6,7,8,9,10;11,12,13,14,15;16,17,18,19,20;21,22,23,24,25];
>>B=[3,0,16;17,-6,9;0,23,-4;9,7,0;4,13,11];
>>C=A*B
(2)
>>D=C(3:
end,2:
end)
(3)
Namevalueclass
A<5×5double>double
B<5×3double>double
C<5×3double>double
D[520397;705557;…]double
4.
(1)
>>X=100:
999;
>>length(find(~rem(X,21)))
(2)
>>y='ahHGJjskjJjsdsaSa';
>>M=y(find(y<='z'&y>='a'))
实验二
实验三
1x=input('请输入x的值')
ifx<0&x~=-3
y=x^2+x-6;
elseifx>=0&x<5&x~=2&x~=3
y=x^2-5*x+6;
else
y=x^2-x-1;
end
y
2
(1)
x=input('请输入成绩')
ifx>=90&x<=100
disp('等级为A')
elseifx>=80&x<=89
disp('等级为B')
elseifx>=70&x<=79
disp('等级为C')
elseifx>=60&x<=69
disp('等级为D')
elseifx>=0&x<=59
disp('等级为E')
else
disp('成绩有误')
end
(2)
x=input('请输入成绩')
switch(x)
case{90,100}
disp('等级为A')
case{80,89}
disp('等级为B')
case{70,79}
disp('等级为C')
case{60,69}
disp('等级为D')
case{0,59}
disp('等级为E')
otherwise
disp('成绩错误')
end
>>
3
Num=input('输入员工的工号:
');
Time=input('该工号员工的工作时数:
');
salary=input('该员工底薪:
');
ifTime>=120
salary=salary*(1+0.15);
elseifTime<=60
salary=salary-700;
elseifTime>60&Time<120
salary=Time*84;
end
salary=salary
实验四
1:
程序:
n=input('输入n:
');
fora=1:
n;
b(a)=1/(a^2);
end
pi=sqrt(6*sum(b))
2:
.
程序:
sum=0;
forn=1:
200;
b(n)=1/(2*n-1);
sum=sum+b(n);
if(sum>=3)
break;
end
end
n-1
4.
程序:
f
(1)=1;max=0;
f
(2)=0;min=0;
f(3)=1;c=0;
zhengshu=0;
fushu=0;
zero=0;
forn=4:
100
f(n)=f(n-1)-2*f(n-2)+f(n-3);
end
form=1:
100
if(f(m)>max)
max=f(m);
end
if(f(m) min=f(m); end c=sum(f); if(f>0) zhengshu=zhengshu+1; end if(f(m)<0) fushu=fushu+1; end if(f(m)==0) zero=zero+1; end if(f(m)>0) zhengshu=zhengshu+1; end if(f(m)<0) fushu=fushu+1; end if(f(m)==0) zero=zero+1; end disp('最大值') max disp('最小值') min disp('和') c disp('正数的个数') zhengshu disp('负数的个数') fushu disp('零的个数') zero 实验五 1,函数文件: 文件名: Yunsuan 程序: function[z,d,s,c]=yunsuan(a,b) z=exp(a+b*i); d=log(a+b*i); s=sin(a+b*i); c=cos(a+b*i); disp(['复数',num2str(a),'+',num2str(b),'i的指数为: ',num2str(z),',对数为: ',num2str(d),',正弦为: ',num2str(s),',余弦为: ',num2str(c)]) 函数调用: 程序: clear x=input('请输入复数实部的值'); y=input('请输入复数虚部的值'); [x,y]=yunsuan(x,y); 2函数文件: 文件名: zhu 程序: function[p,q,w]=zhu(p,q,w) A=[p*cos(w),-p,-sin(w),0;p*sin(w),0,cos(w),0;0,q,-sin(w),0;0,0,-cos(w),1]; B=[0;9.8*p;0;9.8*q]; X=A\B 函数调用: 程序: clear m1=input('m1='); m2=input('m2='); th=input('theta='); theta=pi*th/180; [m1,m2,theta]=zhu(m1,m2,theta); 4函数文件: 文件名: hanshu 程序: function[y]=hanshu(x) y=1/((x-2).^2-0.1)+1./((x-3).^4+0.01); 实验六 1、 程序: 、 x=0: 0.02*pi: 2*pi; y=(0.5+3*sin(x)./(1+x.*x)).*cos(x) plot(x,y) 2、 (1) 程序: x=(0: 0.02*pi: 2*pi)'; y1=x.*x; y2=cos(2*x); y3=y1.*y2; plot(x,y1,'y-'),holdon plot(x,y2,'c-'),holdon plot(x,y3,'k-'),holdoff (2) 程序: x=(0: 0.02*pi: 2*pi)'; y1=x.*x; y2=cos(2*x); y3=y1.*y2; subplot(2,2,1) stairs(x,y1); plot(x,y1,'r-'),holdon subplot(2,2,2) stairs(x,y2); plot(x,y2,'b-'),holdon subplot(2,2,3) stairs(x,y3); plot(x,y3,'k-'),holdoff (3) (1)绘制条形图 程序: x=(0: 0.2*pi: 2*pi)'; y1=x.*x; y2=cos(2*x); y3=y1.*y2; subplot(2,2,1) stairs(x,y1); bar(x,y1) subplot(2,2,2) stairs(x,y2); bar(x,y1) subplot(2,2,3) stairs(x,y3); bar(x,y3) (2)绘制梯形图 程序: x=(0: 0.05*pi: 2*pi)'; y1=x.*x; y2=cos(2*x); y3=y1.*y2; subplot(2,2,1) stairs(x,y1); subplot(2,2,2) stairs(x,y2); subplot(2,2,3) stairs(x,y3); (4)杆图: 程序: x=(0: 0.05*pi: 2*pi)'; y1=x.*x; y2=cos(2*x); y3=y1.*y2; subplot(2,2,1) stem(x,y1); subplot(2,2,2) stem(x,y2); subplot(2,2,3) stem(x,y3); (5)填充图 程序: x=(0: 0.05*pi: 2*pi)'; y1=x.*x; y2=cos(2*x); y3=y1.*y2; subplot(2,2,1) fill(x,y1,'r'); subplot(2,2,2) fill(x,y2,'b'); subplot(2,2,3) fill(x,y3,'g'); 3、: 程序: x=linspace(-5,5,1000); y=[]; forx0=x ifx0<=0 y=[y,(x0+sqrt(pi))./exp (2)]; else y=[y,0.5*log(x0+sqrt(1+x0.*x0))]; end end plot(x,y) 4、 程序: a=input('请输入a的值') b=input('请输入b的值') n=input('请输入n的值') theta=0: 0.01: 2*pi; rho=a.*sin(b+n.*theta) polar(theta,rho) 5、 程序: x=-linspace(-5,5,21); y=linspace(0,10,31); [x,y]=meshgrid(x,y) z=cos(x).*cos(y).*exp(-0.25*sqrt(x.^2+y.^2)) subplot(1,2,1); title('surf(x,y,z)'); surf(x,y,z); subplot(1,2,2); title('surfc(x,y,z)'); surfc(x,y,z); 实验八 1、 程序: x=rand(1,30000); EX=mean(x)%返回向量x的算术平均值 S=std(x)%返回向量x的标准方差 MAX=max(x)%返回向量x的最大元素 MIN=min(x)%返回向量x的最小元素 p=sum(x>0.5)/30000%满足条件的百分比 2、 程序: p=100*rand(100,5); [MAXOBJECT,U]=max(p)%求每门课的最高分,U记录最大元素的行号 [MAXOBJECT,V]=min(p)%求每门课的最低分,V记录最小元素的行号 EX=mean(p)%每门课的平均分 S=std(p)%每门课的标准方差 SUM_1=sum(p,2);%求每名学生的总成绩 [MAXSUM,W]=max(SUM)%求五门课总分的最高分,W记录最高分学生的序号 [MINSUM,Y]=min(SUM)%求五门课总分的最低分,Y记录最低分学生的序号 [zcj,xsxh]=sort(SUM_1); zcj(end: -1: 1); xsxh(end: -1: 1); 实验九 2、程序: I1=quad('sqrt(cos(t.^2)+4*sin(2*t).^2+1)',0,2*pi) I2=quad('(log(1+x))./(1+x.*x)',0,1) 3、 程序: A=[6,5,-2,5;9,-1,4,-1;3,4,2,-2;3,-9,0,2] b=[-4;13;1;11] x=inv(A)*b 实验十 1、 程序: x=sym('6'); y=sym('5'); z=(x+1)/(sqrt(3+x)-sqrt(y)) 2、分解因式 (1) 程序: symsxy; A=x^4-y^4; factor(A) (2) 程序: factor(sym('5135')) 3、化简表达式 (1) 程序: symsbeta1beta2 y=sin(beta1)*cos(beta2)-cos(beta1)*sin(beta2) simple(y) (2) 程序: symsx y=(4*x^2+8*x+3)/(2*x+1) simple(y) 5、用符号方法求下列极限或导数 (1) 程序: symsx f=(x*(exp(sin(x))+1)-2*(exp(tan(x))-1))/(sin(x)) limit(f) (2) 程序: symsx y=(sqrt(pi)-sqrt(acos(x)))/(sqrt(x+1)); limit(f,x,-1,'right') (3) 程序: symsx y=(1-cos(2*x))/x; y1=diff(y) y2=diff(y,x,2) 6、用符号方法求下列积分 (1) 程序: symsx f=1/(1+x^4+x^8) int(f) (2) 程序: symsx f=1/(((asin(x))^2)*sqrt(1-x^2)) int(f) (3) 程序: symsx f=(x^2+1)/(x^4+1) int(f,x,0,inf) (4) 程序: symsx f=exp(x)*(1+exp(x))^2 y=int(f,x,0,log (2)) double(y)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 程序设计 应用 第二 卫国 实验 部分
![提示](https://static.bdocx.com/images/bang_tan.gif)