Matlab数学实验报告.docx
- 文档编号:14890520
- 上传时间:2023-06-28
- 格式:DOCX
- 页数:18
- 大小:93.64KB
Matlab数学实验报告.docx
《Matlab数学实验报告.docx》由会员分享,可在线阅读,更多相关《Matlab数学实验报告.docx(18页珍藏版)》请在冰点文库上搜索。
Matlab数学实验报告
数学实验报告
姓名:
班级:
学号:
第一次实验任务
过程:
a=1+3i;b=2-i;
结果:
a+b=3.0000+2.0000i
a-b=-1.0000+4.0000i
a*b=5.0000+5.0000i
a/b=-0.2000+1.4000i
过程:
x=-4.5*pi/180;y=7.6*pi/180;
结果:
sin(abs(x)+y)/sqrt(cos(abs(x+y)))=0.2098
心得:
对于matlab中的角度计算应转为弧度。
(1)过程:
x=0:
0.01:
2*pi;y1=sin(x);y2=cos(x);y3=exp(x);y4=log(x);
plot(x,y1,x,y2,x,y3,x,y4)plot(x,y1,x,y2,x,y3,x,y4)
结果:
(2)过程:
>>subplot(2,2,1)
>>plot(x,y1)
>>subplot(2,2,2)
>>plot(x,y2)
>>subplot(2,2,3)
>>plot(x,y3)
>>subplot(2.2.4)
>>subplot(2,2,4)
>>plot(x,y4)
结果:
心得:
在matlab中,用subplot能够实现在同一页面输出多个坐标系的图像,应注意将它与holdon进行区别,后者为在同一坐标系中划出多条曲线。
5、随机生成一个3x3矩阵A及3x2矩阵B,计算
(1)AB,
(2)对B中每个元素平方后得到的矩阵C,(3)sinB,(4)A的行列式,(5)判断A是否可逆,若可逆,计算A的逆矩阵,(6)解矩阵方程AX=B,(7)矩阵A中第二行元素加1,其余元素不变,得到矩阵D,计算D。
过程:
A=fix(rand(3,3).*10);
B=fix(rand(3,3).*10);
结果:
(1)A*B=
2228
4964
76100
(2)C=B.^2
C=
14
916
2536
(3)sin(B)
ans=
0.84150.9093
0.1411-0.7568
-0.9589-0.2794
(4)a=det(A)
a=0
(5)由det(A)存在,且det(A)=0,故A不可逆
(6)X无解
(7)l=A(2,1),l1=l+1,m=A(2,2),m1=m+1,n=A(2,3),n1=n+1
D=A,D(2,1)=l1,D(2,2)=m1,D(2,3)=n1
D=
123
567
789
心得:
熟练掌握matlab中对数组的操作方法尤为重要,其中,除了要牢记一些十分方便的函数外,应注意对数组的每一个元素进行操作时,应在运算符前加“.”号,如.^.*等。
6、设y=(x2+excosx+[x])/x,分别计算x=1,3,5,7.4时y的值。
其中[x]表示x的取整函数。
过程:
x=1,y=(x^2+exp(x)*cos(x)+fix(x))/x
结果:
x=1
y=3.4687
x=3
y=-2.6282
x=5
y=14.4198
x=7.4
y=105.2995
7、已知某地区1-12月份的平均气温为5,8,13,20,25,28,31,33,27,24,19,10,请绘图表示。
过程:
:
T=[5,8,13,20,25,28,31,33,27,24,19,10];
bar(T)
结果:
第2次实验任务
1.某大学本科生就业情况如下:
458人考入研究生,60人签到外企或合资企业,184人签到国营大企业,87人签到私人企业,13人自主创业,画出饼图表示之;
过程:
x=[458601848713];
pie3(x)
答案:
2.求极限
过程:
symsxa;
limit('((x-a)/(x+a))^x',x,inf)
答案:
ans=exp(-2*a)
3.求极限
过程:
>>symsx;
>>limit('tan(x)^(1/log10(x))',x,0,'right')
答案:
ans=exp
(1)
4.设
过程:
symsxqt;
q=limit('(1+1/x)^(2*x)',x,inf);
df=diff('q^t',t)
答案:
df=q^t*log(q)
5.展开多项式
;
过程:
symsab;
F=(a+b)^3+(b-1)^2+a+2*b;
expand(F)
答案:
ans=a^3+3*a^2*b+3*a*b^2+b^3+b^2+1+a
6.分解因式
过程:
symsx;
F=x^5+3*x^4+2*x^2+x-1;
factor(F)
答案:
ans=x^5+3*x^4+2*x^2+x-1
7.求方程
的根
过程:
symsx;
C=solve(‘x^3-2*x+1=0’)
答案C=
1
1/2*5^(1/2)-1/2
-1/2-1/2*5^(1/2)
8.
过程:
>>symsx;
>>f=sqrt(x+sqrt(x+sqrt(x)));
>>df=diff(f)
答案:
df=1/2/(x+(x+x^(1/2))^(1/2))^(1/2)*(1+1/2/(x+x^(1/2))^(1/2)*(1+1/2/x^(1/2)))
9.求不定积分
过程:
symsx;
F=1/(sqrt(2*x+3)+sqrt(2*x-1));
Nf=int(F,x,)
答案:
Nf=-1/12*(2*x-1)^(3/2)+1/12*(2*x+3)^(3/2)
过程:
symsx;
F=1/(1+cos(x));
Nf=int(F)
答案:
Nf=tan(1/2*x)
10求定积分
过程:
symsx;
F=1/(x*sqrt(log(x)*(1-log(x)));
Nf=int(F,x,exp(1/2),exp(3/4))
答案:
ans=
asin(2*log(595882530372511)-96*log
(2)-1)-asin(2*log(1856295125090727)-100*log
(2)-1)
11.解方程组
过程:
symsxy;
eq1=sym(‘2*x+y=8’);
eq2=sym(‘x-3*y=1’);
[x,y]=solve(eq1,eq2)
答案:
x=25/7
y=6/7
12.求和
过程:
symsk;
symsum(1/k^2,k,1,20)
答案:
ans=17299975731542641/10838475198270720
13.
过程:
symsx;
taylor(cos(2*x),x,15,pi/6)
答案:
ans=1/2-3^(1/2)*(x-1/6*pi)-(x-1/6*pi)^2+2/3*3^(1/2)*(x-1/6*pi)^3+1/3*(x-1/6*pi)^4-2/15*3^(1/2)*(x-1/6*pi)^5-2/45*(x-1/6*pi)^6+4/315*3^(1/2)*(x-1/6*pi)^7+1/315*(x-1/6*pi)^8-2/2835*3^(1/2)*(x-1/6*pi)^9-2/14175*(x-1/6*pi)^10+4/155925*3^(1/2)*(x-1/6*pi)^11+2/467775*(x-1/6*pi)^12-4/6081075*3^(1/2)*(x-1/6*pi)^13-4/42567525*(x-1/6*pi)^14
14.编写猜数游戏程序:
首先由计算机随机产生一个[1,100]之间的一个整数,然后由用户猜测所产生的这个数。
根据用户猜测的情况给出不同的提示,如果猜测的数大于产生的数,则显示“High”,小于则显示“Low”,等于则显示“Youwon!
”,同时退出游戏。
用户最多有7次机会。
过程:
a=fix(rand
(1)*100);
fori=1:
1:
7
x=input('pleaseinputthenumberyouguess');
ifx>a
fprintf('high');
continue;
end
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Matlab 数学 实验 报告