中南大学材料学院MATLAB考试重点题目老师在最后几节课勾的题目汇总Word下载.docx
- 文档编号:22362178
- 上传时间:2023-02-03
- 格式:DOCX
- 页数:23
- 大小:584.32KB
中南大学材料学院MATLAB考试重点题目老师在最后几节课勾的题目汇总Word下载.docx
《中南大学材料学院MATLAB考试重点题目老师在最后几节课勾的题目汇总Word下载.docx》由会员分享,可在线阅读,更多相关《中南大学材料学院MATLAB考试重点题目老师在最后几节课勾的题目汇总Word下载.docx(23页珍藏版)》请在冰豆网上搜索。
y=rand(10,1);
bar(x,y);
x=0:
0.1*pi:
y=x.*sin(x);
feather(x,y)
P105
x=[2,4,6,8];
pie(x,{'
math'
'
english'
chinese'
music'
});
P107
[x,y]=meshgrid(-2:
0.1:
plot3(x,y,z)
P109
x=-8:
0.5:
8;
y=x'
;
a=ones(size(y))*x;
b=y*ones(size(x));
c=sqrt(a.^2+b.^2)+eps;
z=sin(c)./c;
mesh(z)
P110
[X,Y]=meshgrid([-4:
4]);
%也可以是meshgrid(-4:
4,-4:
4)
Z=sqrt(X.^2+Y.^2);
meshc(X,Y,Z)
p113
x=0:
pi/20:
pi*3;
r=5+cos(x);
[a,b,c]=cylinder(r,30);
mesh(a,b,c)
P113
[a,b,c]=sphere(40);
t=abs(c);
surf(a,b,c,t);
axis('
equal'
square'
colormap('
hot'
P118
x=1:
y=sin(x);
plot(x,y);
xlabel('
x(0-2\pi)'
FontWeight'
bold'
ylabel('
y=sin(x)'
title('
正弦函数'
P123
y=cos(x);
z=sin(x);
plot(x,y,'
-*'
holdon
plot(x,z,'
-o'
plot(x,y+z,'
-h'
legend('
sin(x)'
cos(x)'
sin(x)+cos(x)'
0)
holdoff%这里的holdoff主要是等下次调用时,以上的两个图形将不再保留
P124
subplot(2,2,1);
plot(x,sin(x),'
);
subplot(2,2,2);
plot(x,cos(x),'
subplot(2,2,3);
sin(x).*cos(x)'
plot(x,sin(x).*cos(x));
%注意x是数组,故用sin(x)与其他表达式相乘的时候也要表达成sin(x).*cos(x)
subplot(2,2,4);
plot(x,sin(x)+cos(x),'
-k'
P222
xi=0:
.25:
yi=interp1(x,y,xi);
%注意是插值法函数interp1
o'
xi,yi)
P227
x=[0.5,1.0,1.5,2.0,2.5,3.0];
y=[1.75,2.45,3.81,4.80,8.00,8.60];
a=polyfit(x,y,2)
x1=[0.5:
0.05:
3.0];
y1=a(3)+a
(2)*x1+a
(1)*x1.^2;
*'
plot(x1,y1,'
-r'
P228
x=[1925313844];
y=[19.032.349.073.398.8];
x1=x.^2
x1=[ones(5,1),x1'
]
ab=x1\y'
x0=[19:
44];
y0=ab
(1)+ab
(2)*x0.^2;
clf%这里的clf是清除以前的图片
plot(x0,y0,'
P233
functiony=fun(t)
y=exp(-0.5*t).*sin(t+pi/6);
%注意这边有点,当数组与数组相乘时,其中一个数组应表达成有.的形式
d=pi/1000;
d:
3*pi;
nt=length(t);
y=fun(t);
sc=cumsum(y)*d;
scf=sc(nt)
scf=
0.9016
z=trapz(y)*d
z=
0.9008
P237
functionf=fun(x)
f=exp(-x/2);
%编写函数时让系统自动形成m文件
quadl('
fun'
1,3,1e-10)
ans=
0.7668
P244
x=sym('
x'
diff(sin(x^2))
2*cos(x^2)*x
P274
functionf=f(x,y)
f=[-21;
998-999]*y+[2*sin(x);
999*(cos(x)-sin(x))];
ode23('
f'
[010],[2,3]);
刚性比
a=[-21;
998-999];
b1=max(abs(real(eig(a))));
b2=min(abs(real(eig(a))));
s=b1/b2
s=
1000
P275
dsolve('
Df=f+sin(t)'
f(pi/2)=0'
)
ans=
-1/2*cos(t)-1/2*sin(t)+1/2*exp(t)/(cosh(pi)+sinh(pi))^(1/2)
D2y=-a^2*y'
y(0)=1,Dy(pi/a)=0'
cos(a*t)
S=dsolve('
Dx=y'
Dy=-x'
x(0)=0'
y(0)=1'
S=
x:
[1x1sym]
y:
S.x
sin(t)
S.y
cos(t)
几点对dsolve的解释
dsolve('
Dy=x+1'
y(0)=0'
x*t+t%此时matlab软件自动将t看为默认变量,而将x看成一个常数
P246
a=[0.40960.12340.36780.2943
0.22460.38720.40150.1129
0.36450.19200.37810.0643
0.17840.40020.27860.3927];
b=[0.4043
0.1550
0.4240
-0.2557];
a\b
-0.1819
-1.6630
2.2172
-0.4467
P265
functiony=fc(x)
y
(1)=x
(1)-0.7*sin(x
(1))-0.2*cos(x
(2));
y
(2)=x
(2)-0.7*cos(x
(1))+0.2*sin(x
(2));
x0=[0.50.5];
fsolve('
fc'
x0)%这里的x0为初始值
0.52650.5079
比较fzero函数
例如求x^2+x-1=0的解
fzero('
x^2+x-1'
0.5)%这里的0.5是初始值,初始值只是影响运行时间,对算出的结果影响不大
0.6180
P307
x=[0.2360.2380.2480.2450.243
0.2570.2530.2550.2540.261
0.2580.2640.2590.2670.262];
anova1(x'
)%anova1是一个对单因素试验的方差分析,它返回的是原假设均值相等的几率值
1.3431e-005
P308
x=[58.200056.200065.3000
52.600041.200060.8000
49.100054.100051.6000
42.800050.500048.4000
60.100070.900039.2000
58.300073.200040.7000
75.800058.200048.7000
71.500051.000041.4000];
anova2(x,2)
0.00350.02600.0001
P309
x=[100:
10:
190];
y=[45515461667074788589];
[a,b]=polyfit(x,y,1)
0.4830-2.7394
b=
R:
[2x2double]
df:
8
normr:
2.6878
y=a
(1)*x+a
(2);
plot(x,y)
Fsolve的功能
y=[y
(1)y
(2)];
x0)
0.52650.5079
P310
x=normrnd(0,1,100000,1);
normplot(x)
练习题
1.求函数在指定点的数值导数
y=[xx.^2x.^3;
12*x3*x.^2;
026*x];
x=1;
eval(diff(y))
123
026
006
x=2;
1412
0212
x=3;
1627
0218
2.求下列函数导数
(1)
y=x^10+10^x+(log(10))/log(x);
diff(y)
10*x^9+10^x*log(10)-2592480341699211/1125899906842624/log(x)^2/x
(2)
y=log(1+x);
eval(diff(y,2))%在x=1的条件下对y表达式求两次导数后导函数的值
-0.2500
3.用数值方法求下列积分
首先先讲一下trapz的用法,如下题
t=0:
0.001:
1;
y=t;
trapz(t,y)
0.5000
0.01:
5;
y=(x.^2).*sqrt(2*x.^2+3);
trapz(x,y)
232.8066
x=pi/4:
pi/3;
y=x./(sin(x).^2);
0.3810
第三题拟合曲线题
x=[0:
1];
y=[-0.4471.9783.286.167.087.347.669.569.489.3011.2];
a=polyfit(x,y,2);
x=[0.05:
1.05];
y=a(3)+a
(2)*x+a
(1)*x.^2%注意x要在y前先赋值,不然y不会运行为最新的x对呀的y值y=0.95034.38757.03988.90739.989910.2876
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 中南 大学 材料 学院 MATLAB 考试 重点 题目 老师 最后 几节课勾 汇总