《MATLAB及应用》实验报告2
MATLAB实验报告
MATLAB实验报告实验报告课程名称 MATLAB基础及应用专业班级电子xxxx姓名学号电气与信息学院实验二 MATLAB 数值计算（一）一实验目的：1．掌握数组的创建与运算方法；2. 掌握矩阵的创建与运算方法；3. 掌握数组的运算方法和矩阵运算方法的区别；4．掌握线性方程的求解方法二实验装置：计算机三实验内容：1．数组的创建和运算创建两个含5个元素的一维数组，并求这两个数组的四则运算。
2．矩阵的创建和运算（1）创建两个3×3的矩阵，并求这两个矩阵的四则运算。
（2）创建一个4×4的矩阵，并求这个矩阵的行列式值。
（3）线性方程的求解解方程组=6613753467294x 。
四实验要求:写出实验程序与仿真结果。
1. >> a=linspace(2,10,5) a =2 4 6 8 10>> b=linspace(1,9,5) b =1 3 5 7 9>> a+bans =3 7 11 15 19>> a-bans =1 1 1 1 1>> a.*bans =2 12 30 56 90>> a./bans =2.0000 1.3333 1.2000 1.1429 1.1111 >> a.\b ans =0.5000 0.7500 0.8333 0.8750 0.9000 2.(1) >> a=[1 1 1;2 2 2;3 3 3];>> b=[4 4 4;5 5 5;6 6 6];>> a+bans =5 5 57 7 79 9 9>> a-bans =-3 -3 -3-3 -3 -3-3 -3 -3>> a*bans =15 15 1530 30 3045 45 45>> a/bWarning: Matrix is singular to working precision. ans =NaN NaN NaNNaN NaN NaNNaN NaN NaN>> a\bWarning: Matrix is singular to working precision. ans =NaN NaN NaNNaN NaN NaNInf Inf Inf(2).>> a=[1 1 1 1;2 2 2 2;3 3 3 3;4 4 4 4];>> det(a)ans =(3).>> a=[4 9 2;7 6 4;3 5 7];>> b=[13;6;6];>> x=inv(a)*bx =-0.53081.7109-0.1374实验三 MATLAB数值计算（二）一实验目的：1．掌握多项式的创建与运算方法；2. 掌握基本的数据分析方法；二实验装置：计算机三实验内容：1．多项式创建输入系数矢量，创建多项式x^3-2*x^2+5*x+3。
MATLAB语言及其应用实验报告.doc
MATLAB语言及其应用实验报告.doc一、实验目的1、熟悉MATLAB语言及其基本操作；2、掌握利用MATLAB进行数据分析和可视化的方法；3、掌握MATLAB应用于科学计算的基本方法。
二、实验环境MATLAB R2018a三、实验内容及方法3.1 实验内容1、打开MATLAB环境；2、读入数据文件，并对数据进行分析和处理；3、通过MATLAB绘制图表，对数据进行可视化。
四、实验细节及流程4.1 读入数据文件本实验使用的数据文件为sales.xlsx，其中包括2018年各个季度的销售数据。
首先，通过以下命令导入数据：data=xlsread('sales.xlsx');4.2 数据分析和处理1、计算各季度销售总额通过以下代码计算每个季度的销售额并求和，得到每年的销售总额：Q1=data(:,2);TotalSales=sum([Q1 Q2 Q3 Q4],2);2、计算增长率根据每年的销售总额，计算出每年的增长率。
具体代码如下：GrowthRate(1)=0;for i=2:length(TotalSales)GrowthRate(i)=((TotalSales(i)-TotalSales(i-1))/TotalSales(i-1))*100; endSalesGrowth=[TotalSales GrowthRate];3、计算每个季度的均值和标准差meanQ1=mean(Q1);stdQ1=std(Q1);4、计算出每年第一个季度的销售额所占比例首先，我们将第一个季度的销售额单独提出来，具体代码如下：4.3 数据可视化1、柱形图对于销售总额，使用柱形图进行可视化，具体代码如下：结果如下图所示：图1 销售总额2、线性图3、箱形图boxplot([Q1 Q2 Q3 Q4],{'Q1','Q2','Q3','Q4'});图3 每个季度的销售额4、饼图pie(FirstQSalesRatio(:,1));五、结论本实验通过对销售数据的分析和可视化，得出以下结论：1、2018年销售总额呈逐年上升趋势，其中2017年到2018年的增长率最高；2、每年第一个季度的销售额所占比例在40%至45%之间，与其他季度相比，显著高于其他季度；3、2018年第二季度的销售额经过调整后，表现出了相对较高的波动。
matlab实验报告
Matlab实验报告实验二图像处理一、实验目的（1）通过应用MA TLAB语言编程实现对图像的处理，进一步熟悉MATLAB软件的编程及应用；（2）通过实验进一步掌握图像处理的基本技术和方法。
二、实验内容及代码㈠.应用MA TLAB语言编写显示一幅灰度图像、二值图像、索引图像及彩色图像的程序，并进行相互之间的转换首先，在matlab页面中的current directory下打开存放图像的文件夹。
1.显示各种图像⑴显示彩色图像：①代码：>> mousetif=imread('tif.TIF');>> image(mousetif)显示截图：②代码：>> mousetif=imread('tif.TIF');>> imshow(mousetif)显示截图：③代码：mousetif=imread('tif.TIF');subimage(mousetif)显示截图：显示截图：⑵显示二值图像①代码：>> I=imread('单色bmp.bmp');>> imagesc(I,[0 2])显示截图：②代码：>> I=imread('单色bmp.bmp');>> imshow(I,2)显示截图：③代码：>> I=imread('单色bmp.bmp');>> subimage(I)显示截图：⑶显示灰度图像①代码：>> I1=imread('256bmp.bmp');>> imagesc(I1,[0,256])显示截图：代码：>> I1=imread('256bmp.bmp');>> colormap(gray);>> subplot(1,2,1);>> imagesc(I1,[0,256]);>> title('灰度级为[0 256]的mouse.bmp图');>> subplot(1,2,2);>> imagesc(I1,[0,64]);>> colormap(gray);>> title('灰度级为[0 64]的mouse.bmp图');显示截图：⑷显示索引图像代码：>> mousetif=imread('tif.TIF');>> [x,map]=rgb2ind(mousetif,0.7);>> subplot(1,2,1);>> imshow(mousetif);>> title('彩色图像');>> subplot(1,2,2);>> imshow(x,map);>> title('索引图');截图显示：彩色图像索引图2.图像之间的转换⑴索引图像转灰度图像代码：>> II1=ind2gray(x,map);>> imshow(x,map);>> title('索引图像');>> imshow(II1);>> title('灰度图像');截图显示：索引图像灰度图像⑵ 索引图像转彩色图像⑶ 索引图像转二值图像⑷ 灰度图像转索引图像代码：>> I1=imread('256bmp.bmp');>> subplot(1,2,1);>> imagesc(I1,[0,256]); >> colormap(gray); >> title('灰度图');>> [x,map]=gray2ind(I1,256); >> subplot(1,2,2); >> image(x);>> colormap(map); >> title('索引图')；显示截图：灰度图2004006008001000100200300400500600索引图2004006008001000100200300400500600⑸灰度图像转彩色图像⑹灰度图像转二值图像⑺彩色图像转灰度图像⑻彩色图像转索引图像代码：>> mousetif=imread('tif.TIF');>> [x,map]=rgb2ind(mousetif,0.7);>> subplot(1,2,1);>> imshow(mousetif);>> title('彩色图像');>> subplot(1,2,2);>> imshow(x,map);>> title('索引图');显示截图：彩色图像索引图⑼彩色图像转二值图像⑽二值图像转索引图像⑾二值图像转彩色图像⑿二值图像转灰度图像㈡应用MA TLAB工具箱演示一幅图像的傅里叶变换、离散余弦变换，观察其频谱图。
MATLAB实验报告(打印版)
MATLAB实验报告班别：09电气3班姓名：李嘉明学号：200924122301实验一/二 MATLAB的基础操作3．自定义一个变量，然后分别用8种不同的数字显示格式显示查看。
>> a=[1.1,1.2,1.3;1.4,1.5,1.6;1.7,1.8,1.9]a =1.1000 1.2000 1.30001.4000 1.5000 1.60001.7000 1.8000 1.9000>> format long>> aa =1.100000000000000 1.200000000000000 1.3000000000000001.400000000000000 1.500000000000000 1.6000000000000001.700000000000000 1.800000000000000 1.900000000000000>> format short e>> aa =1.1000e+000 1.2000e+000 1.3000e+0001.4000e+000 1.5000e+000 1.6000e+0001.7000e+000 1.8000e+000 1.9000e+000>> format long e>> aa =1.100000000000000e+000 1.200000000000000e+000 1.300000000000000e+0001.400000000000000e+000 1.500000000000000e+000 1.600000000000000e+0001.700000000000000e+000 1.800000000000000e+000 1.900000000000000e+000 >> format hex>> aa =3ff199999999999a 3ff3333333333333 3ff4cccccccccccd3ff6666666666666 3ff8000000000000 3ff999999999999a3ffb333333333333 3ffccccccccccccd 3ffe666666666666>> format +>> aa =+++++++++>> format rat>> aa =11/10 6/5 13/107/5 3/2 8/517/10 9/5 19/10>> format short>> aa =1.1000 1.2000 1.30001.4000 1.5000 1.60001.7000 1.8000 1.90005、己举例比较给矩阵的全行赋值，取出矩阵中某一个或几个元素，然后组成一个新的矩阵。
Matlab实验报告
实验一：Matlab操作环境熟悉一、实验目的1．初步了解Matlab操作环境。
2．学习使用图形函数计算器命令funtool及其环境。
二、实验内容熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算：1．单函数运算操作。
求下列函数的符号导数(1)y=sin(x); (2) y=(1+x)^3*(2-x);求下列函数的符号积分(1)y=cos(x);(2)y=1/(1+x^2);(3)y=1/sqrt(1-x^2);(4)y=(x1)/(x+1)/(x+2)求反函数(1)y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2));代数式的化简(1)(x+1)*(x-1)*(x-2)/(x-3)/(x-4);(2)sin(x)^2+cos(x)^2;(3)x+sin(x)+2*x-3*cos(x)+4*x*sin(x);2．函数与参数的运算操作。
从y=x^2通过参数的选择去观察下列函数的图形变化(1)y1=(x+1)^2(2) y2=(x+2)^2(3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6)y6=x^2/23．两个函数之间的操作求和(1)sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5乘积(1)exp(-x)*sin(x) (2) sin(x)*x商(1)sin(x)/cos(x); (2) x/(1+x^2); (3) 1/(x-1)/(x-2);求复合函数(1)y=exp(u) u=sin(x) (2) y=sqrt(u) u=1+exp(x^2)(3) y=sin(u) u=asin(x) (4) y=sinh(u) u=-x实验二：MATLAB基本操作与用法一、实验目的1.掌握用MATLAB命令窗口进行简单数学运算。
matlab程序设计与应用(第二版) 刘卫国 课后实验报告
0 0
0 0
4.0000 0
0 9.0000
A,H,Trace,Rank,Norm 结果 H = -3.2509e+004 Trace =33 Rank =5 Norm =25.3208 4、 A=[-29,6,18;20,5,12;-8,8,5]; [V,D]=eig(A); 结果 A = -29 20 -8 5、 5 8 6 18 12 5
0 0 2、 H=hilb(5); P=pascal(5); Hh=det(H); Hp=det(P); Th=cond(H); Tp=cond(P); H,P,Hh,Hp,Th,Tp 结果 Hh = 3.7493e-012 Hp =1 Th = 4.7661e+005 Tp =8.5175e+003 3、 A=fix(10*rand(5)); H=det(A); Trace=trace(A); Rank=rank(A); Norm=norm(A);
b=?1 ans =1 x1 =0.6180 r10 =0.6180 r20 =-1.6180 4、 for i=1:100 if i==1 f(i)=1; elseif i==2 f(i)=0; elseif i==3 f(i)=1; else f(i)=f(i-1)-2*f(i-2)+f(i-3); end end a=max(f); b=min(f); c=sum(f); disp(['最大值',num2str(a)]); disp(['最小值',num2str(b)]); disp(['各数之和',num2str(c)]); zhengshu=length(find(f>0)); zeroshu=length(find(f==0)); fushu=length(find(f<0)); 结果 最大值；4.377633e+011 最小值：-8.994121e+011 各数之和：-7.427456e+011 zhengshu =49 zeroshu = 2 fushu =49 5、 s=0;
MATLAB 语言及其应用 实验报告
《MATLAB 语言及其应用》实验报告目录实验一Matlab 使用方法和程序设计........................实验二控制系统的模型及其转换.............................实验三控制系统的时域、频域和根轨迹分析...........实验四动态仿真集成环境-Simulink.........................实验一Matlab使用方法和程序设计一、实验目的1、掌握Matlab软件使用的基本方法；2、熟悉Matlab的数据表示、基本运算和程序控制语句3、熟悉Matlab绘图命令及基本绘图控制4、熟悉Matlab程序设计的基本方法二、实验内容：1、帮助命令使用help命令，查找 sqrt（开方）函数的使用方法；答：help sqrt2、矩阵运算（1）矩阵的乘法已知A=[1 2;3 4]; B=[5 5;7 8];求A^2*B解：>> A=[1 2;3 4];B=[5 5;7 8];A^2*Bans =105 115229 251（2）矩阵除法已知 A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];A\B,A/B解：>> A=[1 2 3;4 5 6;7 8 9];B=[1 0 0;0 2 0;0 0 3];>> A\Bans =1.0e+016 *-0.4504 1.8014 -1.35110.9007 -3.6029 2.7022-0.4504 1.8014 -1.3511>> A/Bans =1.0000 1.0000 1.00004.0000 2.5000 2.00007.0000 4.0000 3.0000（3）矩阵的转置及共轭转置已知A=[5+i,2-i,1;6*i,4,9-i];求A.', A'解：>> A=[5+i,2-i,1;6*i,4,9-i];>> A.'ans =5.0000 + 1.0000i 0 +6.0000i2.0000 - 1.0000i 4.00001.0000 9.0000 - 1.0000i>> A'ans =5.0000 - 1.0000i 0 -6.0000i2.0000 + 1.0000i 4.00001.0000 9.0000 + 1.0000i（4）使用冒号表达式选出指定元素已知： A=[1 2 3;4 5 6;7 8 9];求A中第3列前2个元素；A中所有列第2，3行的元素；>> A=[1 2 3;4 5 6;7 8 9];>> B=A(1:2,3)B =36>> B1=A(2:3,:)B1 =4 5 67 8 9用magic函数生成一个4阶魔术矩阵，删除该矩阵的第四列>> A=magic(4)A =16 2 3 135 11 10 89 7 6 124 14 15 1>> B=A(:,1:3)B =16 2 35 11 109 7 64 14 153、多项式（1）求多项式4xxp的根2=x)(3-->> p=[1 0 -2 -4];>> r=roots(p)r =2.0000-1.0000 + 1.0000i-1.0000 - 1.0000i（2）已知A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ，求矩阵A的特征多项式;>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ;>> p=poly(A)p =1.0000 -6.9000 -77.2600 -86.1300 604.5500把矩阵A作为未知数代入到多项式中；>> A=[1.2 3 5 0.9;5 1.7 5 6;3 9 0 1;1 2 3 4] ;>> syms x;P=x^3-2*x-4;>> P1=subs(P,x,A)P1 =-4.6720 17.0000 111.0000 -5.0710111.0000 -2.4870 111.0000 200.000017.0000 707.0000 -4.0000 -5.0000-5.0000 0 17.0000 52.00004、基本绘图命令（1）绘制余弦曲线 y=cos(t)，t∈[0，2π]（2）在同一坐标系中绘制余弦曲线y=cos(t-0.25)和正弦曲线y=sin(t-0.5)，t∈[0，2π]解：（1）程序为：t=[0:0.05:2*pi];y=cos(t);plot(t,y)（2）程序如下：x=[0:0.05:2*pi];y1=cos(x-0.25);y2=sin(x-0.5);plot(x,y1,x,y2)绘制曲线如下：5、基本绘图控制绘制[0，4π]区间上的x1=10sint曲线，并要求：（1）线形为点划线、颜色为红色、数据点标记为加号；（2）坐标轴控制：显示范围、刻度线、比例、网络线（3）标注控制：坐标轴名称、标题、相应文本；程序为：>> t=[0:0.05:4*pi];x1=10*sin(t);plot(t,x1,'r-.+')grid on;axis([0,15,-10,10]);title('曲线x1=10sint');xlabel('T轴');ylabel('X1轴');set(gca,'xminortick','on');set(gca,'yminortick','on');6、基本程序设计（1）编写命令文件：计算1+2+…+n<2000 时的最大n值；>> s=1;n=1;while(s<2000-n),n=n+1;s=s+n;end,[s,n]ans =1953 62（2）编写函数文件：分别用for和while循环结构编写程序，求2的0到n次幂的和。
matlab实验报告
实验二特殊函数与图形
一、实验目的及意义
掌握用Matlab软件绘制简单曲线、曲面图形，并通过绘制一些特殊函数的图形，更加深入地理解相关函数的性质，了解函数的性态。
二、实验内容
平面作图、空间作图，比较数值作图与符号作图的异同。
三、实验步骤
1.在D盘建立一个自己的文件夹;
2.开启软件平台——MATLAB，将你建立的文件夹加入到MA TLAB的搜索路径中。
3.利用帮助了解函数plot, surf, ezmesh，ezsurf等的功能和用法。
4.开启MATLAB编辑窗口,键入你编写的M文件（命令文件或函数文件）；
5.保存文件（注意将文件存入你自己的文件夹）并运行；
6.若出现错误，修改、运行直到输出正确结果；
7.写出实验报告，并浅谈学习心得体会。
四、实验要求与任务
根据实验内容和步骤，完成以下具体实验，按要求写出实验报告。
1．作出下图所示的三维图形：
图1
提示：图形为圆环面和球面的组合.
2．作出下图所示的墨西哥帽子及其剪裁图形：
图2
3．画出椭球面、双叶双曲面、单叶双曲面．
4．若要求田螺线的一条轴截面的曲边是一条抛物线：0=y 时25x z =．试重新设计田螺线的参数方程，并画出该田螺线．
5．作出下图所示的马鞍面（颜色为灰色，并有一个标题：“马鞍面”）：
图3
6．绘制黎曼函数的图形，要求分母的最大值n 的数值由键盘输入（提示：使用
input 语句）．。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
核科学技术学院实验报告实验项目名称MATLAB符号计算所属课程名称MATLAB及应用实验类型上机实验实验日期12月日指导教师谢芹班级学号姓名成绩一、实验名称MATLAB符号计算二、实验目的（1）掌握定义符号对象的方法（2）掌握符号表达式的运算法则以及符号矩阵运算（3）掌握求符号函数极限及导数的方法（4）掌握求符号函数定积分和不定积分的方法三、实验原理1. 函数极限及导数的方法（1）函数极限：limit(F,x,a) 求符号函数f(x)的极限值。
即计算当变量x趋近于常数a时，f(x)函数的极限值。
（2）limit(f)：求符号函数f(x)的极限值。
符号函数f(x)的变量为函数findsym(f)确定的默认变量；没有指定变量的目标值时，系统默认变量趋近于0，即a=0的情况。
（3）limit(f,x,a,'right')：求符号函数f的极限值。
'right'表示变量x从右边趋近于a。
（4）limit(f,x,a,‘left’)：求符号函数f的极限值。
‘left’表示变量x从左边趋近于a。
2. 微分：diff(s)：没有指定变量和导数阶数，则系统按findsym函数指示的默认变量对符号表达式s求一阶导数。
diff(s,'v')：以v为自变量，对符号表达式s求一阶导数。
diff(s,n)：按findsym函数指示的默认变量对符号表达式s求n阶导数，n为正整数。
diff(s,'v',n)：以v为自变量，对符号表达式s求n阶导数。
3. 函数定积分和不定积分的方法：int(s)：没有指定积分变量和积分阶数时，系统按findsym函数指示的默认变量对被积函数或符号表达式s求不定积分。
int(s,v)：以v为自变量，对被积函数或符号表达式s求不定积分。
int(s,v,a,b)：求定积分运算。
a,b分别表示定积分的下限和上限。
梯形法：trapz(x,y)：x为分割点构成的向量，y为被积函数在分割点上的函数值构成的向量；抛物线法：quad(f,a,b,tol)，f是被积函数，[a,b]是积分区间，tol是精度。
4. 求和及泰勒级数展开的方法：（1）求和symsum(s,v,n,m) 其中s 表示一个级数的通项，是一个符号表达式。
v 是求和变量，v 省略时使用系统的默认变量。
n 和m 是求和的开始项和末项。
（2）泰勒级数展开 taylor(f,v,n,a) 该函数将函数f 按变量v 展开为泰勒级数，展开到第n 项(即变量v 的n-1次幂)为止，n 的缺省值为6。
v 的缺省值与diff 函数相同。
参数a 指定将函数f 在自变量v=a 处展开，a 的缺省值是0。
四、实验内容1. 求下列极限： 求极限前先定义符号变量（1） 10arctan lim x x F x→= （2）1201lim()1x x x F x →+=- （3）320ln(1)limsin x x x F x →+= （4）43111lim()11x F x x→=--- （5）552lim(1)x x t F ax →∞=+ 2. 求下列函数的导数：（1） 3cos cos3y x x =- （2）1sin ln y x x x = （3）21sin x xe y x-= （4）3cos x y e x = （5） 2sin y x x = （6） 3cos ln x ae t x f t x x ⎡⎤=⎢⎥⎣⎦，求df dx ，22d f dt ，2d f dtdx 3. 求下列函数的积分（1） sin sin sin ax bx cxdx ⎰ （2）(53)4x x x dx +-⎰ （3）120(1)x xe dx x +⎰ （4）[0,1][0,1]1Dx dxdy D xy =⨯+⎰⎰ （5） V zdxdydz ⎰⎰⎰ 由曲面22y xz +=，1=z ，2=z 所围成（6） 21cos x bx dx x e x ⎡⎤⎢⎥⎢⎥⎣⎦⎰ 4. 解下列方程组。
（1）123123123257543232615x x x x x x x x x -+-=-⎧⎪+-=⎨⎪++=⎩（2）1122125sin 4cos 05cos 4sin 0x x x x x x --=⎧⎨-+=⎩ （3）()5()()()()5()dx t x t y t dt dy t y t x t dt⎧=-+⎪⎪⎨⎪=-+⎪⎩ 5. 求下列级数的和（1） 11212n n n I ∞=-=∑ （2）211(21)n I n n ∞==+∑ 6. 泰勒级数展开将函数21()53f x x x =+-展开成2-x 的幂级数 五、实验过程及结果（含源代码）（此处按题号顺序写上所写的程序语句内容以及matlab 中出现的程序运行结果） 1.syms x;f1='arctan(x)/x'F1=limit(f1,x,0)f2='((1+x)/(1-x))^(1/x)'F2=limit(f2,x,0)f3='(x*ln(1+x))/(sin(x^2))'F3=limit(f3,x,0)f4='(1/(1-x))-(1/(1-x^3))'F4=limit(f4,x,1)f5='(1+((2*t)/(a*x)))^(5*x)'F5=limit(f5,x,inf)f1 =arctan(x)/xF1 =1f2 =((1+x)/(1-x))^(1/x)F2 =exp(2)f3 =(x*ln(1+x))/(sin(x^2))F3 =1f4 =(1/(1-x))-(1/(1-x^3))F4 =NaNf5 =(1+((2*t)/(a*x)))^(5*x)F5 =exp((10*t)/a)2.syms x ln(x) cos(x) sin(x) t a; y1=(cos(x))^3-cos(3*x)Y1=diff(y1)y2=x*sin(x)*ln(x)Y2=diff(y2)y3=(x*exp(x)-1)/sin(x)Y3=diff(y3)y4=exp(x)*cos(x)Y4=diff(y4)y5=(x^2)*sin(x)Y5=diff(y5)y6=[a*exp(x) (t^3)*x;t*cos(x) ln(x) ]Y6=diff(x)Y7=diff(x,2)y1 =cos(x)^3 - cos(3*x)Y1 =3*cos(x)^2*diff(cos(x), x) - 3*D(cos)(3*x)y2 =x*ln(x)*sin(x)Y2 =ln(x)*sin(x) + x*ln(x)*diff(sin(x), x) + x*sin(x)*diff(ln(x), x)y3 =(x*exp(x) - 1)/sin(x)Y3 =(exp(x) + x*exp(x))/sin(x) - ((x*exp(x) - 1)*diff(sin(x), x))/sin(x)^2 y4 =exp(x)*cos(x)Y4 =exp(x)*cos(x) + exp(x)*diff(cos(x), x)y5 =x^2*sin(x)Y5 =2*x*sin(x) + x^2*diff(sin(x), x)y6 =[ a*exp(x), t^3*x][ t*cos(x), ln(x)]Y6 =1Y7 =3.syms x y z ;f1='sin(a*x)*sin(b*x)*sin(c*x)';f2='x^5+x^3-(x^(1/2))/4';f3='x*(exp(x))/(1+x^2)';f4='x/(1+x*y)';f='x^2+y^2';f6=sym('[1/x,b*x^2;exp(x),cos(x)]'); int(f1,x)int(f2,x)int(f3,x,0,1)int(int(f4,x,0,1),y,0,1)int(int(int(f,x),y),z,1,2)int(f6,x)ans =(a*b^2*cos(a*x)*sin(b*x)*sin(c*x) - b^3*cos(b*x)*sin(a*x)*sin(c*x) - c^3*cos(c*x)*sin(a*x)*sin(b*x) - a^3*cos(a*x)*sin(b*x)*sin(c*x) +a^2*b*cos(b*x)*sin(a*x)*sin(c*x) + a*c^2*cos(a*x)*sin(b*x)*sin(c*x) + a^2*c*cos(c*x)*sin(a*x)*sin(b*x) + b*c^2*cos(b*x)*sin(a*x)*sin(c*x) + b^2*c*cos(c*x)*sin(a*x)*sin(b*x) +2*a*b*c*cos(a*x)*cos(b*x)*cos(c*x))/(a^4 - 2*a^2*b^2 - 2*a^2*c^2 + b^4 - 2*b^2*c^2 + c^4)ans =x^4/4 - x^(3/2)/6 + x^6/6ans =- (ei(-i)*exp(i))/2 - (ei(i)*exp(-i))/2 + (ei(1 - i)*exp(i))/2 + (ei(1 + i)*exp(-i))/2ans =log(4) - 1ans =(x*y*(x^2 + y^2))/3ans =[ log(x), (b*x^3)/3][ exp(x), sin(x)]4.[x,y,z]=solve('-2*x+5*y-7*z=5','4*x+5*y-2*z=3','2*x+y+6*z=15') [x,y]=solve('x-5*sin(x)-4*cos(y)=0','y-5*cos(x)+4*sin(y)=0') [x,y]=dsolve('Dx=-5*x+y','Dy=-y+5*x')x =-212/81y =307/81z =74/27x =1.8207539129496562592168622750749y =-3.8552617061181725539710730412911x =C2/5 - C1*exp(-6*t)y =C2 + C1*exp(-6*t)5.syms nf1='(2*n-1)/2^n';f2='1/(n*(2*n+1))';I1=symsum(f1,n,1,inf)I2=symsum(f2,n,1,inf)I1 =3I2 =2 - 2*log(2)6.syms xtaylor('1/(x^2+5*x-3)',x,2)ans =(70*(x - 2)^2)/1331 - (9*x)/121 - (531*(x - 2)^3)/14641 + (4009*(x - 2)^4)/161051 - (30240*(x - 2)^5)/1771561 + 29/121六、实验总结通过本次的matlab上机实验，使我了解了符号对象定义的方法，可以较为熟练地运用计算法则以及对符号助阵的求解。