春数学实验基础 作业1 电子版Word文件下载.docx
- 文档编号:21546322
- 上传时间:2023-01-31
- 格式:DOCX
- 页数:16
- 大小:672.85KB
春数学实验基础 作业1 电子版Word文件下载.docx
《春数学实验基础 作业1 电子版Word文件下载.docx》由会员分享,可在线阅读,更多相关《春数学实验基础 作业1 电子版Word文件下载.docx(16页珍藏版)》请在冰豆网上搜索。
0
7
2
(7)C(2,:
)
取C的第二行元素列成行向量
-18-372
(8)C([22],[33])
取C的第二行第三列元素构成一个2*2的矩阵
-3-3
(9)C+10
取C的每一个元素都加10得到一个新的矩阵
1112131015
91871712
151741211
(10)2*C+1
取C的每个元素乘以二在加一构成一个新的矩阵
357111
-117-5155
1115-1153
2.写出下列语句的运行结果:
(1)a=[123;
456;
789];
a([2,1],:
)=a([12],:
)
使a的第一、二行元素对换
a=
456
123
789
(2)a=[123;
a([1,3],:
)=a([22],:
使矩阵的一三行用第二行同列元素替代
(3)a=[123;
a=a([23],:
使a变成只有原二三行的新矩阵
(4)a=eye(3);
b=[123];
a(:
4)=b'
a是3*3的单位矩阵,b是一个3阶行向量,给a加上第四列为b的转置构成一个3*4的新矩阵
1001
0102
0013
3.设
写出下列命令的运行结果,若有不合法的运算,请说明理由.
(1)result=a+b;
3-3
-17
(2)result=a*d;
2-2
-15
(3)result=a.*d;
20
05
(4)result=a*c;
6
-11
(5)result=a.*c;
不合法,因为是点乘,即相同位置上的元素相乘,两矩阵行列数应该相等
(6)result=a\b;
0.6250-0.1250
0.12500.3750
(7)result=a.\b;
0.50000.5000
00.4000
(8)result=a.^b;
2.0000-0.5000
1.000025.0000
4.设
请至少用两种方法生成矩阵A.
提示:
可用MATLAB内部函数diag,eye,zeros等.
第一种:
>
A=eye(10);
A=2012*A
第二种:
a=[1,1,1,1,1,1,1,1,1,1];
a=2012*a;
diag(a)
5.作出平面曲线:
(1)标准正态分布概率密度函数(standardnormaldistributiondensityfunction):
>
ezplot('
1/sqrt(2*pi)*exp(-x^2/2)'
[-3.5,3.5])
(2)在同一直角坐标系中作出三条曲线:
>
x*sin(1/x)'
[-pi,pi]);
holdon
x'
-x'
(3)画出参数曲线:
ezplot('
(cos(t))^3'
'
(sin(t))^3'
[0,pi]);
(4)用‘subplot’命令在6个子窗口中分别画6条曲线:
subplot(2,3,1);
cos(x)'
[-3*pi,3*pi]);
subplot(2,3,2);
cos(2*x)'
subplot(2,3,3);
cos(3*x)'
subplot(2,3,4);
cos(4*x)'
subplot(2,3,5);
cos(5*x)'
subplot(2,3,6);
cos(6*x)'
(6)用‘subplot’命令在6个子窗口中分别画6条极坐标曲线:
ezpolar('
cos(t)'
);
cos(2*t)'
cos(3*t)'
cos(4*t)'
cos(5*t)'
cos(6*t)'
6.作出3-D曲线:
(1)
(2)
(1)>
ezplot3('
sin(t)'
[0,2*pi]);
(2)>
t*sin(t)'
t*cos(t)'
t'
[-12*pi,12*pi]);
7.作出下列曲面:
(meshsurfcontour3surfc)
(1)马鞍面
subplot(2,2,1)
ezcontourf('
x^2-y^2'
[-11-11])>
ezmesh('
[-11-11])
subplot(2,2,2)
ezsurf('
subplot(2,2,3)
ezcontour('
subplot(2,2,4)
ezsurfc('
(2)抛物面
x^2/3+y^2/2'
(3)
旋转曲面
sin(sqrt(x^2+y^2))/sqrt(x^2+y^2)'
(4)平面
2'
(5)圆锥面
先化为参数方程:
x=s*sin(t),y=s*cos(t),z=s*sin(
/4)
然后画出:
s*cos(t)'
s*sin(t)'
s*sin(pi/4)'
[-1,1,0,4*pi])
(6)曲面
(3+0.3*cos(u))*cos(v)'
(3+0.3*cos(u))*sin(v)'
0.3*sin(v)'
[-pi,pi,-pi,pi])
8.编写一个MATLAB文本m文件(scriptm-file)。
要求从键盘输入一个正整数,判断它是否能被3整除,若能,则输出“该整数能被3整除”,否则输出“该整数不能被3整除,余数为…”
functionzhengchu(m)
s=0;
ifmod(m,3)==0
disp('
该整数能被3整除'
else
s=mod(m,3);
disp(['
该整数不能被3整除,且余数是:
'
num2str(s)])
end
9.求满足
的最小m值。
由
故只需找出最小整数m满足
即可(这样还可以消去逐项计算中软件的舍入误差)
functions=sss
m=0;
whilelog(factorial(m+1))<
=100
m=m+1;
s=m;
运行结果:
sss
37
10.自己编写函数计算,:
(1)n!
(2)
并写出调用程序。
(1)functiony=jiecheng(x)
y=1;
fori=1:
x
y=y*i;
end
(2)
functiony=zuhe(m,n)
m1=jiecheng(m);
n1=jiecheng(n);
s=jiecheng(n-m);
y=n1/(m1*s);
11.编写函数计算
的n阶Taylor多项式Pn(x)及误差Rn(x),并画图与
比较(n=1,2,3,4,5):
。
M文件程序:
functionp=Taylor(x,n)
p=0;
r=0;
fori=0:
n
y=x^i/factorial(i);
p=p+y;
r=exp(x)-p;
r
p
Taylor(2,5)%取x=2,n取5的多项式展开
r=
0.1224%误差R
p=
7.2667%Taylor多项式值
●比较图象:
实验小结:
通过本次作业发现很多有关matlab操作的知识已经遗忘,也得到及时的巩固,其中有很多地方不熟练,许多功能平时没有去接触就感觉很陌生,作业中许多题目回答得还是很简略,遇到想优化的地方,想美观的地方都无能为力(主要是对图像的编辑),今后还要加强学习。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 春数学实验基础 作业1 电子版 数学 实验 基础 作业