MATLAB实验上机练习三.docx
《MATLAB实验上机练习三.docx》由会员分享，可在线阅读，更多相关《MATLAB实验上机练习三.docx（8页珍藏版）》请在冰豆网上搜索。
MATLAB实验上机练习三
一、在matlab的命令窗口里完成如下计算，其中t的值分别取-1,0,1，表达式如下：
t=[-101]
y=sqrt（3）/2*exp（-4.*t）.*sin（4*sqrt（3.*t）+pi/3）
t=
-101
y=
1.0e+04*
2.0897+1.2065i0.00010.0000
二、自行产生一个5行5列的数组，分别得到最中间的三行三列矩阵、右下角2行2列矩阵，奇数行矩阵、奇数列矩阵、奇数行奇数列矩阵。
>>t=rand（5）%生成矩阵
A=t（2:
4,2:
4）%中间三行散列矩阵
B=t（4:
5,4:
5）%右下角两行两列矩阵
C=t（1:
2:
end,:
）%奇数行矩阵
D=t（:
1:
2:
end）%奇数列矩阵
E=t（1:
2:
end,1:
2:
end）%奇数行列矩阵
t=
0.16220.60200.45050.82580.1067
0.79430.26300.08380.53830.9619
0.31120.65410.22900.99610.0046
0.52850.68920.91330.07820.7749
0.16560.74820.15240.44270.8173
A=
0.26300.08380.5383
0.65410.22900.9961
0.68920.91330.0782
B=
0.07820.7749
0.44270.8173
C=
0.16220.60200.45050.82580.1067
0.31120.65410.22900.99610.0046
0.16560.74820.15240.44270.8173
D=
0.16220.45050.1067
0.79430.08380.9619
0.31120.22900.0046
0.52850.91330.7749
0.16560.15240.8173
E=
0.16220.45050.1067
0.31120.22900.0046
0.16560.15240.8173
三、用magic产生一个5*5的矩阵，将这个矩阵的第二行与第三行互换位置。
>>T=magic（5）
B=T（2,:
）
T（2,:
）=T（3,:
）
T（3,:
）=B
T=
17241815
23571416
46132022
101219213
11182529
B=
23571416
T=
17241815
46132022
46132022
101219213
11182529
T=
17241815
46132022
23571416
101219213
11182529
三、求方程组的根
>>A=[14-3;25-1;161]
B=[2;11;12]
X=A\B
A=
14-3
25-1
161
B=
2
11
12
X=
4
1
2
四、电路分析常常涉及到对方程组的求解，这些方程常常是利用描述电流进入和离开节点的电流方程，或者描述电路中网络回路上的电压的电压方程得出的。
下图描述了3个回路电压的方程
方程式如下：
假设5个电阻值为已知，2个电压值也为已知，求3个电流值。
>>symsI1I2I3U1U2U3R1R2R3R4R5
[I1I2I3]=solve（'-U1+R1*I1+R2*（I1-I2）=0','R2*（I2-I1）+R3*I2+R4*（I2-I3）=0','R4*（I3-I2）+R5*I3+U2=0',I1,I2,I3）
I1=
（R2*R4*U1-R2*R4*U2+R2*R5*U1+R3*R4*U1+R3*R5*U1+R4*R5*U1）/（R1*R2*R4+R1*R2*R5+R1*R3*R4+R1*R3*R5+R2*R3*R4+R1*R4*R5+R2*R3*R5+R2*R4*R5）
I2=
-（R1*R4*U2-R2*R4*U1+R2*R4*U2-R2*R5*U1）/（R1*R2*R4+R1*R2*R5+R1*R3*R4+R1*R3*R5+R2*R3*R4+R1*R4*R5+R2*R3*R5+R2*R4*R5）
I3=
-（R1*R2*U2+R1*R3*U2+R1*R4*U2+R2*R3*U2-R2*R4*U1+R2*R4*U2）/（R1*R2*R4+R1*R2*R5+R1*R3*R4+R1*R3*R5+R2*R3*R4+R1*R4*R5+R2*R3*R5+R2*R4*R5）
五、根据a=reshape（-4:
4,3,3）做一个矩阵，然后
（1）取出所有大于0的元素构成一个向量（可推广到从一个矩阵里筛选出符合条件的元素组成一个向量）
（2）将原矩阵中大于0的元素正常显示，而小于等于0的元素全部用0来表示（可推广到将原矩阵中不符合条件的全用0来表示，符合条件的值不变）。
（思考：
大于0的正常显示，小于等于0的用-10来表示）
>>a=reshape（-4:
4,3,3）
b=a（find（a>0））
a（a<=0）=0
a=
-4-12
-303
-214
b=
1
2
3
4
a=
002
003
014
六、将矩阵
、
和
组合成两个新矩阵：
（1）组合成一个43的矩阵，第一列为按列顺序排列的a矩阵元素，第二列为按列顺序排列的b矩阵元素，第三列为按列顺序排列的c矩阵元素，即
（2）按照a、b、c的列顺序组合成一个行矢量，即
a=[42;57]
b=[71;83]
c=[59;62]
J=[a（:
1）a（:
2）b（:
1）;b（:
2）c（:
1）c（:
2）]%
（1）
K=J'
T=[K（1,:
）K（2,:
）K（3,:
）]%
（2）
a=
42
57
b=
71
83
c=
59
62
J=
427
578
159
362
K=
4513
2756
7892
T=
451327567892
7、设
和
求A-2B
>>A=[4-31;205;321]
B=[120;-103;231]
C=A-2*B
A=
4-31
205
321
B=
120
-103
231
C=
2-71
40-1
-1-4-1
升级会员 