西安交大数学实验报告.docx
- 文档编号:29076650
- 上传时间:2023-07-20
- 格式:DOCX
- 页数:23
- 大小:256.28KB
西安交大数学实验报告.docx
《西安交大数学实验报告.docx》由会员分享,可在线阅读,更多相关《西安交大数学实验报告.docx(23页珍藏版)》请在冰豆网上搜索。
西安交大数学实验报告
实验报告
班级:
姓名:
学号:
第一次实验
代码:
a=1+3i;b=2-i;
a+b
a-b
a*b
a/b
运行结果:
代码:
x=-4.5/180*pi;
y=7.6/180*pi;
sin(abs(x)+y)/sqrt(abs(x+y))
运行结果:
(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);
运行结果:
(2)代码:
x=0:
0.01:
3
y1=sin(x);
y2=cos(x);
y3=exp(x);
y4=log(x);
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)
plot(x,y4)
运行结果:
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=randn(3,3)
B=randn(3,2)
AB=A*B
C=B.^2
sinB=sin(B)
detA=det(A)
Ani=A^(-1)
X=A^(-1)*B
A(2,:
)=A(2,:
)+1
D=A
Det=det(D)
运行结果:
6、设y=(x2+excosx+[x])/x,分别计算x=1,3,5,7.4时y的值。
其中[x]表示x的取整函数。
代码:
X=[1,3,5,7.4];
Y=(X.^2+exp(X).*cos(X)+fix(X))./X;
y1=Y
(1)
y2=Y
(2)
y3=Y(3)
y4=Y(4)
运行结果:
7、已知某地区1-12月份的平均气温为5,8,13,20,25,28,31,33,27,24,19,10,请绘图表示。
代码:
y=[5,8,13,20,25,28,31,33,27,24,19,10];
x=[1,2,3,4,5,6,7,8,9,10,11,12];
plot(x,y)
xlabel('月份')
ylabel('平均气温')
gridon
运行结果:
代码:
axisequal
ezplot('x^2+y^2=4')
运行结果:
第二次实验
1.某大学本科生就业情况如下:
458人考入研究生,60人签到外企或合资企业,184人签到国营大企业,87人签到私人企业,13人自主创业,画出饼图表示之;
代码:
x=[601848713458];
pie(x)
pie(x,explode)
name={'考入研究生458人','签到外企或合资企业60人','签到国营大企业184人','签到私人企业87人','自主创业13人'};
pie(x,name)
运行结果:
2.求极限
代码:
symsax;
y=((x-a)/(x+a))^x;
limit(y,x,inf)
运行结果:
3.求极限
代码:
symsx;
y=(tan(x))^(1/log(x));
limit(y,x,0,'right')
运行结果:
4.设
代码:
symsxt;
f=(1+1/x)^(2*t*x);
f1=limit(f,x,inf);
diff(f1,t)
运行结果:
5.展开多项式
;
代码:
symsab;
fc=(a+1)^3+(b-1)^2+2*b;
fc1=expand(ps)
运行结果:
6.分解因式
代码:
symsx;
y=x^5-3*x^4+2*x^2+x-1;
factor(y)
运行结果:
7.求方程
的根
代码:
symx;
f=sym('x^3-2*x+1');
solve(f)
运行结果:
8.
代码:
symsx;
y=sqrt(x+sqrt(x+sqrt(x)));
diff(y,x)
运行结果:
9.求不定积分
代码:
symsx;
f=sym('1/(sqrt(2*x+3)+sqrt(2*x-1))');
int(f)
运行结果:
10.
代码:
symsx;
f=sym('1/(x*sqrt(log(x)*(1-log(x))))');
int(f,x,(exp)^(1/2),(exp)^(3/4))
运行结果:
10求定积分
代码:
symsx;
f=sym('1/(x*sqrt(log(x)*(1-log(x))))');
int(f,x,(exp
(1))^(1/2),(exp
(1))^(3/4))
运行结果:
11.解方程组
代码:
A=[2,1;1,-3];
B=[8;1];
X=A\B
运行结果:
12.求和
代码:
symsk;
A=symsum(1/(k^2),k,1,20)
运行结果:
13.
代码:
symsx;
y=taylor(cos(2*x),x,1/6*pi,'Order',15)
运行结果:
y=
(2*3^(1/2)*(x-pi/6)^3)/3-3^(1/2)*(x-pi/6)-(2*3^(1/2)*(x-pi/6)^5)/15+(4*3^(1/2)*(x-pi/6)^7)/315-(2*3^(1/2)*(x-pi/6)^9)/2835+(4*3^(1/2)*(x-pi/6)^11)/155925-(4*3^(1/2)*(x-pi/6)^13)/6081075-(x-pi/6)^2+(x-pi/6)^4/3-(2*(x-pi/6)^6)/45+(x-pi/6)^8/315-(2*(x-pi/6)^10)/14175+(2*(x-pi/6)^12)/467775-(4*(x-pi/6)^14)/42567525+1/2
14.编写猜数游戏程序:
首先由计算机随机产生一个[1,100]之间的一个整数,然后由用户猜测所产生的这个数。
根据用户猜测的情况给出不同的提示,如果猜测的数大于产生的数,则显示“High”,小于则显示“Low”,等于则显示“Youwon!
”,同时退出游戏。
用户最多有7次机会。
代码:
x=fix(rand
(1)*100);
fori=1:
7
y=input('ÇëÊäÈëÒ»¸öÕûÊý');
ifx==y
fprintf('youwon!
')
break
elseify>x
fprintf('High\n')
elseify fprintf('Low\n') end cleary end 运行结果: 第三次实验 1.每门课程考试阅卷完毕,任课教师都要对各班的考试成绩进行统计,统计内容包括: 全班人数,总得分,平均得分,不及格的人数及90分(包括90分)以上的人数.请编制程序解决这一问题,并自给一组数据验证程序的正确性.要求: 使用者在提示下通过键盘输入学生成绩,计算机自动处理后,显示需要的结果. 代码: x=input('请输入全班成绩'); l=length(x); s=sum(x); a=mean(x); n=0; m=0; fori=1: l ifx(i)>90 n=n+1; elseifx(i)<60 m=m+1; end end fprintf('学生总数=%i总得分=%.1f平均得分=%1.f不及格人数=%i90分以上人数=%i',l,s,a,m,n) 运行结果: 2.《数学实验》P67练习1第2题 代码: a=input('请输入一个大于6的偶数'); m=1; n=a-m; x=isprime(m); y=isprime(n); whilex+y<2 m=m+1; n=a-m; x=isprime(m); y=isprime(n); end fprintf('两个质数为a=%.ib=%.i',m,n) 运行结果: 3.编写猜数游戏程序: 首先由计算机随机产生一个[1,100]之间的一个整数,然后由用户猜测所产生的这个数。 根据用户猜测的情况给出不同的提示,如果猜测的数大于产生的数,则显示“High”,小于则显示“Low”,等于则显示“Youwon! ”,同时退出游戏。 用户最多有7次机会。 代码: a=rand; b=fix(a*100); forn=1: 1: 7 c=input('请输入一个1到100间的随机数c='); ifc>b disp('high') elseifc disp('low') elsec==b disp('youwin! ') end end 运行结果: 编写程序解决下列问题 4.Fibonacci数组元素满足Fibonacci规则: (1)用循环语句形成前20个分量的数组; (2)求该数组中第一个大于10000的数. 代码: (1) a (1)=1; a (2)=1; forn=1: 18 a(n+2)=a(n+1)+a(n); end a (2) a (1)=1; a (2)=1; forn=1: 22 a(n+2)=a(n+1)+a(n); end fori=1: 22 whilea(i)>10000 fprint(‘第一个大于10000的数为: %.4f\n’,a(i)) return end end 运行结果: (1) (2) 5. (1)编写一个function函数,调用该函数,可以求一组数的最大值及最小值. (2)通过键盘输入一组数,利用 (1)编写的function函数,得出最大最小数. 代码: (1)functionmaxandmin(a) max=a (1); n=length(a) fori=2: n
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 西安 交大 数学 实验 报告