天津工业大学 电气工程仿真与计算实验报告.docx
《天津工业大学 电气工程仿真与计算实验报告.docx》由会员分享，可在线阅读，更多相关《天津工业大学 电气工程仿真与计算实验报告.docx（13页珍藏版）》请在冰豆网上搜索。
天津工业大学电气工程仿真与计算实验报告
《电气工程仿真与计算》
实验报告
班级：
学号:
姓名：
I.Matlab运算基础
A.先求下列表达式的值，然后显示matlab工作空间的使用情况并保存全部变量。
1.
程序：
%计算表达式的值
x=1+（sin（5））^3;
y=1+csc（3）;
z=log（x^（0.5）+y^（0.5））;
z1=z/（1+exp
（2））;
答案：
z1=0.1382
2.角度x=[304560]，求x的正弦、余弦、正切和余切。
程序：
%角度x=[304560]，求x的正弦、余弦、正切和余切。
x=[304560];%定义x的值
fori=1:
3{
sind（x（i））;%求sinx
cosd（x（i））;%求cosx
tand（x（i））;%求tanx
cotd（x（i））;%求cotx
}
end
答案：
x=30。
时：
ans=
[0.5000]
[0.8660]
[0.5774]
[1.7321]
x=45。
时：
ans=
[0.7071]
[0.7071]
[1]
[1]
x=60。
时：
ans=
[0.8660]
[0.5000]
[1.7321]
[0.5774]
B.已知
，求下列表达式的值：
1.求解A，B，C转置
2.（A+C）*B和A*（C+I）（其中I为单位矩阵）
3.A*B和A.*C
4.A^（0.5）和C.^3
5.A/C；A\B；A\C
6.[A;B’]；[A,C]；[A;C]
7.分别取矩阵A的1、3列与C的第2行，组成一个矩阵。
（要求C的第2行作为该矩阵的第3行）
程序：
%matlab矩阵练习
%初始矩阵定义
%DesignbyJIANGTAO
A=[414-3;12554;2117];
B=[3;6;-1];
C=[19-1i13+1i12i;1+2i4-1i7-8i;4-3i6+9i9i];
%第一问
Transposed_A=A';
Transposed_B=B';
Transposed_C=C';
%第二问
D1=（A+C）*B;
D2=A*（C+eye（3））;
%第三问
E1=A*B;
E2=A.*C;
%第四问
F1=A^（0.5）;
F2=C.^（3）;
%第五问
G1=A/C;
G2=A\B;
G3=A\C;
%第六问
H1=[A;B'];
H2=[A,C];
H3=[A;C];
%第七问
I=[A（:
1）';A（:
3）';C（2,:
）];
答案：
第一问
Transposed_A=
4122
14551
-3417
Transposed_B=
36-1
Transposed_C=
19+1i1-2i4+3i
13-1i4+1i6-9i
0-12i7+8i0-9i
第二问
D1=
234-9i
382+8i
43+36i
D2=
82+33i104-37i95-91i
311+86i455-7i389-260i
109-51i133+154i24+169i
第三问
E1=
99
362
-5
E2=
1.0e+02*
0.7600-0.0400i1.8200+0.1400i0.0000-0.3600i
0.1200+0.2400i2.2000-0.5500i0.2800-0.3200i
0.0800-0.0600i0.0600+0.0900i0.0000+1.5300i
第四问
F1=
1.42511.6139-0.6626
1.36367.26520.4296
0.34950.03824.1491
F2=
6802-1082i2158+506i-0-1728i
-11-2i52-47i-1001-664i
-44-117i-1242+243i-0-729i
第五问
G1=
0.1289+0.3962i-0.0135-1.0310i0.6180-1.1220i
0.4450+1.3294i0.6383-2.7996i2.1515-4.2091i
0.1981-0.3497i1.0412+1.3348i-0.3769+0.5737i
G2=
0.8838
-0.0722
-0.1586
G3=
14.3893-1.6295i9.7841+2.4812i-1.3242+11.8438i
-3.0283+0.3925i-2.0126-0.5793i0.4066-2.6782i
-1.2794-0.0078i-0.6797+0.2716i0.1319-0.7064i
第六问
H1=
414-3
12554
2117
36-1
H2=
4+0i14+0i-3+0i19-1i13+1i0+12i
12+0i55+0i4+0i1+2i4-1i7-8i
2+0i1+0i17+0i4-3i6+9i0+9i
H3=
4.0000+0.0000i14.0000+0.0000i-3.0000+0.0000i
12.0000+0.0000i55.0000+0.0000i4.0000+0.0000i
2.0000+0.0000i1.0000+0.0000i17.0000+0.0000i
19.0000-1.0000i13.0000+1.0000i0.0000+12.0000i
1.0000+2.0000i4.0000-1.0000i7.0000-8.0000i
4.0000-3.0000i6.0000+9.0000i0.0000+9.0000i
第七问
I=
4.0000+0.0000i12.0000+0.0000i2.0000+0.0000i
-3.0000+0.0000i4.0000+0.0000i17.0000+0.0000i
1.0000+2.0000i4.0000-1.0000i7.0000-8.0000i
II.M语言程序设计
分别用if和switch语句实现以下程序
某商场对顾客所购买的商品实行打折销售，标准如下：
price<200没有折扣
200<=price<5003%折扣
500<=price<10005%折扣
1000<=price<25008%折扣
2500<=price<500010%折扣
5000<=price14%折扣
输入所售商品的价格，求其实际销售价格
程序：
%分类计算
%DesignbyJIANGTAO
%输入
P=input（'Entertheprice:
'）;
%计算
if（P<200）
Price=P;%无折扣
elseif（P<500）
Price=P*（1-0.03）;%3%折扣
elseif（P<1000）
Price=P*（1-0.05）;%5%折扣
elseif（P<2500）;
Price=P*（1-0.08）;%8%折扣
elseif（P<5000）;
Price=P*（1-0.1）;%10%折扣
elsePrice=P*（1-0.14）;%14%折扣
end
end
end
end
end
%输出
disp（['Thespend=',num2str（Price）]）;
答案：
Entertheprice:
4700
Thespend=4230
III.Matlab绘图练习
A.绘制下列曲线（x，y）
1.
，a=3，b=4；要求椭圆的x、y轴均不少于1000个点
程序：
%绘制椭圆
t=0:
0.001*pi:
2*pi;
x=3*cos（t）;
y=4*sin（t）;
plot（x,y）;
title（'绘制x^2/9+y^2/16=1'）;%加标题
axis（[-44-55]）;%纵横坐标范围
gridon
2.
x=[0,5]，x的取值不少于500个点
程序：
x=0:
0.005:
5;
y=2.^（x/3）.*cos（x）+log（1+x.^2）/exp（x/2）;
plot（x,y）;%绘制图像
title（'绘制y=2^（^x^/^3^）*cos（x）+log（1+x^2）/e^（^x^/^2^）,x=[0,5]'）%标题
gridon
3.已知g=0.5，theta=0:
pi/20:
2*pi，其中Radius=2*g*（1+cos（theta）），用极坐标（polar）方式绘制该函数Radius（theta）曲线
程序：
%极坐标方式绘图
g=0.5;
theta=0:
pi/20:
2*pi;
Radius=2*g*（1+cos（theta））;
polar（theta,Radius）;%绘制图像
title（'¼«×ø±ê·½Ê½»æÖÆRadiusÇúÏß'）;%标题
B.在同一图形窗口采用子图形以及单图多波形的形式分别绘制sin（x），cos（x）,exp（x）
程序：
%在同一图形窗口采用子图形以及单图多波形的形式分别绘制sin（x）,cos（x）,exp（x）
%DesignbyJIANGTAO
x=-pi:
0.001*pi:
pi;
y1=sin（x）;
y2=cos（x）;
y3=exp（x）;
%绘制第1个窗口
figure
（1）;
subplot（3,1,1）;%共3行1列第一行
plot（x,y1）;%绘制y=sinx
title（'y=sinx'）;%标题
gridon%绘制虚线
subplot（3,1,2）;%第2行
plot（x,y2）;
gridon
title（'y=cosx'）;
subplot（3,1,3）;%第3行
plot（x,y3）;
title（'y=e^x'）;
gridon
%绘制第2个窗口
figure
（2）;
holdon
plot（x,y1,'b-'）;
plot（x,y2,'b--'）;
plot（x,y3,'k-'）;
title（'绘制y=sinxy=cosxy=e^x'）
gridon
legend（'sinx','cosx','e^x'）%注释
holdoff
升级会员 