2设计算:
>> A12 B46C80D35E40
>> Tatan(2*pi*A+E(2*pi*B*C))D
T
04112
>> xpi180*45
>> (sin(x)+sqrt(35))72^(15)
ans
28158
3设计算
>> xpi180*45
>> (sin(x)+sqrt(35))72^(15)
ans
28158
4设计算:
>> a567b7811
>> exp(a+b)log10(a+b)
ans
63351e+005
5计算时值
>> x3
>> ysqrt(x)6*(x+1x)+(x32)^2(x+77)^3
y
182679
6已知圆半径15求直径周长面积
>> r15d2*r
>> cpi*2*r
c
942478
>> spi*r*r
s
7068583
7已知某三角形三边边长85146184求该三角形面积
提示:
中:分三角形三边边长
>>a85b146c184
>> s(a+b+c)2
>> areasqrt(s*(sa)*(sb)*(sc))
area
606106
第二章
1设矩阵
求:
>> A[3 1 12 1 21 2 3]B[1 1 12 1 01 1 1]
(1)>> 2*A+B
ans
7 3 1
6 1 4
3 3 7
(2)>> 4*A*A3*B*B
ans
42 21 38
40 19 46
40 33 56
(3)>> A*B
ans
6 1 2
6 1 0
8 4 2
(4)>> B*A
ans
4 0 0
4 1 0
2 2 2
(5)>> A*BB*A
ans
2 1 2
2 2 0
6 6 0
2设三阶矩阵满足
中
求矩阵
>> A[13 0 0 0 14 00 0 17]
>> B6*inv(inv(A)eye(3))
B
3 0 0
0 2 0
0 0 1
3设中4阶单位矩阵4阶矩阵转置中
求矩阵
>> B[1 2 3 20 1 2 30 0 1 20 0 0 1]
>> C[1 2 0 10 1 2 00 0 1 20 0 0 1]
>> A(inv(2*eye(4)inv(C)*B)*inv(C))'
A
1 0 0 0
2 1 0 0
1 2 1 0
0 1 2 1
4设二阶矩阵满足
中
求矩阵
>> A[2 11 2]B[0 22 0]
>> X(2*A+B)2
X
2 0
2 2
5求解线性方程组
>> A[2 3 0 21 5 2 13 1 1 14 1 2 2]
>> B[82712]
>> xA\B
x
30000
00000
10000
10000
6求解元六次方程根
>> A[3 12 4 7 0 8 1]
>> xroots(A)
x
38230
05275 + 08497i
05275 08497i
05007 + 06749i
05007 06749i
01234
7求项式结果
>> a[3 12 4 7 0 8 1]A[3]
>> bpoly(A)
>> c[1 0 5 0]
>> dconv(bc)
>> [divrest]deconv(ad)
div
3 21 52
rest
0 0 0 103 55 788 1
第三章
1 求极限
>>xsym('x')
>>limit((cos(sqrt(x)))^(pix)x0'right')
ans
exp(12*pi)
2 求极限
>>xsym('x')
>>fsym('(3*sin(x)+x^2*cos(1x))((1+cos(x))*log(1+x))')
>>limit(fx0)
ans
32
3 求极限
>>xsym('x')
>>fsym('(sqrt(4*x^2+x1)+x+1)sqrt(x^2+sin(x))')
>>limit(fxinf)
ans
1
4 求极限
>>syms x y
>>fsym('(x^2+y^2)^(x^2*y^2)')
>>limit(limit(fx0)y0)
ans
1
5已知
求y’
xsym('x')
ysym('(tan(sqrt(x+sqrt(x+sqrt(2*x)))))^2')
diff(yx)
ans
tan((x+(x+2^(12)*x^(12))^(12))^(12))*(1+tan((x+(x+2^(12)*x^(12))^(12))^(12))^2)(x+(x+2^(12)*x^(12))^(12))^(12)*(1+12(x+2^(12)*x^(12))^(12)*(1+12*2^(12)x^(12)))
6 已知
求y’
>>xsym('x')
>>ysym('cos(x^2)*sin(1x)^2')
>>diff(yx)
ans
2*sin(x^2)*x*sin(1x)^22*cos(x^2)*sin(1x)*cos(1x)x^2
7求积分
>>xsym('x')
>>ysym('sqrt(sin(x)sin(x)^3)')
>>int(yx0pi)
ans
43
8 求积分
>>xsym('x')
>>ysym('sqrt((x+1)(x1))x')
>>int(yx)
ans
((1+x)(x1))^(12)*(x1)*(atan(1(1+x^2)^(12))log(x+(1+x^2)^(12)))((1+x)*(x1))^(12)
9求列微分方程通解
>>dsolve('D2y+4*Dy+4*yexp(2*x)''x')
ans
12*exp(2*x)*(2*C2+2*x*C1+x^2)
10 求微分方程
>>dsolve('x^2*Dy+x*yy^2''y(1)1')
ans
x*exp(1x)(exp(1x)+exp(1x*t)*xexp(1x*t))
第四章
1已知量请绘图表示
>> A[1 2 4 0 5 10 11 21 3 1]
>> plot(A)
2绘制单位圆
提示:axis(square’)命令保证图形横坐标刻度例相
>> ezplot('x^2+y^21')
>> axis square
3绘制伏安特性曲线:假设分151020
>>R[1 5 10 20]
>>I00110
>>UI'*R
>>plot(IU)
>>legend R1 R5 R10 R20
>>grid on
4某区年中月均气温均降雨量表4 5示请画出图形求标注出坐标轴数点位置数点等
表4 – 5 温度降雨量数
月份
湿度
降雨量
1
02
46
2
23
36
3
87
21
4
185
29
5
246
30
6
321
27
7
368
22
8
371
25
9
283
43
10
178
34
11
64
21
12
32
37
>>month112
>>degree[02 23 87 185 246 321 368 371 283 178 64 32]
>>rain[46 36 21 29 30 27 22 25 43 34 21 37]
>>plot(monthdegree'+'monthrain'p')
>> month112
>>degree[02 23 87 185 246 321 368 371 283 178 64 32]
>>rain[46 36 21 29 30 27 22 25 43 34 21 37]
>>plot(monthdegree'+'monthrain'p')
>>grid on
>>xlabel('month')
>>ylabel('degree&rainfall')
>>text(5246'\leftarrowdegree')
>>text(321'\leftarrowrainfall')
>>title('Temperature&rainfall')
5已知矩阵
请绘图表示
>>Aones(57)
>>A(2426)2
>>plot(A)
6绘制三维立体图
>> x5015
>> yx
>> [XY]meshgrid(xy)
>> ZX^2+Y^2
>> surf(XYZ)
7绘制函数形成立体图通改变观察点获该图形坐标面面投影
>>theta(00012)*pi
>>phitheta
>>xzeros(length(theta)length(phi))
>>yx
>>zx
>>for i1length(theta)
for j1length(phi)
x(ij)2*sin(theta(i))*cos(phi(j))
y(ij)3*sin(theta(i))*sin(phi(j))
z(ij)4*cos(theta(i))
end
end
>>axis square
>>subplot(221)
>>mesh(xyz)
>>title('三维')
>>subplot(222)
>>mesh(xyz)
>>view(900)
>>title('x轴视图')
>>subplot(223)
>>mesh(xyz)
>>view(00)
>>title('y轴视图')
>>subplot(224)
>>mesh(xyz)
>>view(090)
>>title('z轴视图')
第五章
1命令文件画出列分段函数表示曲线
>>x20012
>>for ii1401
if(x(ii)<0)
y(ii)x(ii)+1
elseif (x(ii)<1&x(ii)>0)
y(ii)1
elseif (x(ii)>1)
y(ii)x(ii)^3
end
end
>>plot(xy)
2计算述分段函数值求够根户值输入程序出相应结果
>>function yfenduan(x)
fenduan计算分段函数
x意实数
>>if(x<0)
yx+1
elseif(x<1)&(x>0)
y1
else
yx^3
end
3编程求中
提示:结束条件中通项公式满足求精度极值
>>function yarcsin(x)
>>if abs(x)>1
disp('输入参数绝值必须1')
return
end
>>n0
>>y0
>>ux
>>while u>eps
u(fac(2*n)*x^(2*n+1))(2^(2*n)*(fac(n)^2*(2*n+1)))
yy+u
nn+1
end
>>function yfac(n)
if(n<0)
disp('输入参数必须等0整数')
return
end
if n0|n1
y1
else
yn*fac(n1)
end
4求解鸡兔笼问题:鸡兔子关笼子里已知头36脚100求笼关少兔子少鸡?
n0
m36
while (2*n+4*m~100)
nn+1
mm1
end
nm
ex5_4
n
22
m
14
5求2~999间中时满足列条件然数
(1)该数位数字奇数
(2)该数素数
>> for i2999
arem(i10)
brem(fix(i10)10)
cfix(i100)
if rem((a+b+c)2)&isprime(i)
disp(i)
end
end
3 5 7 23 29 41 43 47 61 67 83 89 113 131 137 139 151 157 173 179 191 193 197 199 223 227 229 241 263 269 281 283 311 313 317 331 337 353 359 373 379 397 401 409 421 443 449 461 463 467 487 557 571 577 593 599 601 607 641 643 647 661 683 719 733 739 751 757 773 797 809 821 823 827 829 863 881 883 887 911 919 937 953 971 977 991 997
6编写判断意输入正整数否素数函数文件
>>function yissushu(n)
>>y1
>>for ii2fix(n2)
if (rem(nii)0)
y0
end
end
7编写阶函数调该函数生成阶表图示:
1
1
2
2
3
6
4
24
5
120
6
720
7
5040
8
40320
9
362880
10
3628800
>>function yfac(n)
>>if(n<0)
disp(请输入正整数’)
return
end
>>if n0|n1
y1
else
yn*fac(n)
end
>>for a110
bfac(a)
A(a)b
end
>>B110
>> aA'
>> bB'
>> c[b a]
c
1 1
2 2
3 6
4 24
5 120
6 720
7 5040
8 40320
9 362880
10 3628800
第六章
3利Simulink仿真两正弦信号相计算
4设计模拟口变化情况模型
根口学理表示某年口数目中表示年份年口口繁殖速率新增资源满足体数目间动力学方程差分方程描述:
面差分方程中出口变化系统非线性离散系统现假设口初始值10 000口繁殖率11新增资源满足体数目1 000 000请建立口动态变化系统模型0~200年间口数目变化趋势进行仿真
MATLAB实验报告
文档香网(httpswwwxiangdangnet)户传
《香当网》用户分享的内容,不代表《香当网》观点或立场,请自行判断内容的真实性和可靠性!
该内容是文档的文本内容,更好的格式请下载文档