matlab符号运算实验报告
MATLAB实验报告
MATLAB实验报告江苏科技⼤学《MATLAB实践》实验报告学号：1341901111姓名：符⾸夫计算机科学与⼯程学院2016.5.31实验⼀ MATLAB的基本使⽤⽅法实验⽬的和要求：通过完成实验⼀，掌握MATLAB的基本使⽤⽅法。
实验内容：(⼀）练习数据和符号输⼊法，将前⾯的命令在命令窗⼝中执⾏通过。
1）>>52）>>x=[1 2 3 4]3）>>g=[1 2 3 4];h=[4 3 2 1]>>s1=g+h, s2=g.*h, s3=g.^hS1=5 5 5 5S2=4 6 6 4S3=1 8 9 4（⼆）输⼊A=[715;256;315],B=[111;222;333],在命令窗⼝中执⾏下列表达式，掌握其含义。
源程序：>>A=[7 1 5;2 5 6;3 1 5];B=[1 1 1;2 2 2;3 3 3];>>A(2,3)>>A(:,2)>>A(:,1:2:3) % A数组中第⼀列和第三列>>A(3,:) %A数组中第三⾏>>A(:,3).*B(:,2) %A中第三列和B中第⼆列相乘>> A(:,3)*B(2,:) %A中第三列和B中第⼆⾏相乘竖着排，形成3*3矩阵>>A*B %A 与B数组相乘>>A.*B>>A^2>>A.^2 % A中所有数平⽅>>B/A %B矩阵除以A矩阵>>B./A % B中数除以A中对应位置数(三) 输⼊C=.1：2：20，则C （i ）表⽰什么？其中i=1，2，3…10；源程序：>>c=1:2:20; >>c(1) >>c(10)实验⼆ MATLAB的数值计算实验⽬的和要求：通过完成实验⼆，掌握MATLAB的数值计算。
实验内容：(⼀) ⽤⼆分法计算多项式⽅程X^3-2X-5=0在[0,3]内的⼀个根。
MATLAB实验报告(1-4)
信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。
2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。
4.学会运用MATLAB进行连续信号时移、反折和尺度变换。
5.学会运用MATLAB进行连续时间微分、积分运算。
6.学会运用MATLAB进行连续信号相加、相乘运算。
7.学会运用MATLAB进行连续信号的奇偶分解。
二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。
三、实验内容1.MATLAB软件基本运算入门。
1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。
2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn 为结束值。
矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。
2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。
3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。
2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。
实验四MATLAB数值计算与符号计算
实验四 MATLAB数值计算与符号计算一、实验目的1.掌握数据插值和曲线拟合的方法2.掌握求数值导数和数值积分的方法3.掌握代数方程数值求解的方法4.掌握常微分方程数值求解的方法5.掌握求解优化问题的方法6.掌握求符号极限、导数和积分的方法7.掌握代数方程符号求解的方法8.掌握常微分方程符号求解的方法二、实验原理1．数据插值a) 一维数据插值 Y1=interp1(X,Y,X1,’method’)b) 二维数据插值 Z1=interp2(X,Y,Z,X1,Y1,’method’)2．曲线拟合[P,S]=polyfit(X,Y,m)3．符号对象的建立（1）符号量名=sym(符号字符串)：建立单个的符号变量或常量；（2）syms arg1 arg2,…,argn：建立n个符号变量或常量。
4．基本符号运算（1）基本四则运算：+，-，*，\，^（2）分子与分母的提取：[n,d]=numden(s)（3）因式分解与展开：factor(s),expand(s)（4）化简：simplify, simple(s)5．符号函数及其应用（1）求极限：limit(f,x,a)（2）求导数：diff(f,x,a);（3）求积分：int(f,v)三、实验内容1．按下表用3次样条方法插值计算0～900范围内整数点的正弦值和0～750范围内整数点的正切值，然后用5次多项式拟合方法计算相同的函数值，并将两种计算结果进行比较。
x2=0:75;y1=sin(pi.*x1./180);y2=tan(pi.*x2./180);;a=interp1(x1,y1,45,'cublic')b=interp1(x1,y1,45,'cublic')p1=polyfit(x1,y1,5)p2=polyfit(x2,y2,5)c1=polyval(p1,x1);c2=polyval(p2,x2);subplot(2,1,1);plot(x1,c1,':o',x1,y1,'r');subplot(2,1,2);plot(x2,c2,':o',x2,y2,'r');10203040506070802．（1）求函数33()sin cos f x x x =+在点,,,6432x ππππ=的数值导数。
实验四MATLAB符号运算
实验四MATLAB符号运算实验四MATLAB符号运算⼀、实验⽬的：1、掌握定义符号对象的⽅法；2、掌握符号表达式的运算法则以及符号矩阵运算。
3、掌握求符号函数极限及导数的⽅法。
4、掌握求符号函数定积分和不定积分的⽅法。
⼆、实验原理1、符号常量、符号变量、符号表达式的创建(1) 使⽤sym( )创建输⼊以下命令，观察Workspace 中A、B、f是什么类型的数据，占⽤多少字节的内存空间。
>>A=sym('1') %符号常量>>B=sym('x') %符号变量>>f=sym('2*x^2+3y-1') %符号表达式>>clear>>f1=sym('1+2') %有单引号，表⽰字符串>>f2=sym(1+2) %⽆单引号>>f3=sym('2*x+3')>>f4=sym(2*x+3) %为什么会出错>>x=1>>f4=sym(2*x+3)通过看MATLAB 的帮助可知，sym( )的参数可以是字符串或数值类型，⽆论是哪种类型都会⽣成符号类型数据。
(2) 使⽤syms 创建>>clear>>syms x y z %注意观察x,y,z都是什么类型的，它们的内容是什么>>x,y,z>>f1=x^2+2*x+1>>f2=exp(y)+exp(z)^2>>f3=f1+f2通过以上实验，知道⽣成符号表达式的第⼆种⽅法：由符号类型的变量经过运算(加减乘除等)得到。
⼜如：>>f1=sym('x^2+y +sin(2)')>>syms x y>>f2=x^2+y+sin(2)>>x=sym('2') , y=sym('1')>>f3=x^2+y+sin(2)>>y=sym('w')>>f4=x^2+y+sin(2)（3）符号矩阵创建>>syms a1 a2 a3 a4>>A=[a1 a2;a3 a4]>>A(1),A(3)或者>>B=sym('[ b1 b2 ;b3 b4] ')>>c1=sym('sin(x) ')>>c2=sym('x^2')>>c3=sym('3*y+z')>>c4=sym('3 ')>>C=[c1 c2; c3 c4]2、符号算术运算(1) 符号量相乘、相除符号量相乘运算和数值量相乘⼀样，分成矩阵乘和数组乘。
MATLAB实验报告
中南民族大学计算机科学学院MATLAB实验报告题目MATLAB实验年级 2010专业计算机科学与技术指导教师李波小组成员(姓名学号) 实验类型综合型2014年4月22 日一、实验安排1.实验目的1.掌握字符串的生成和操作，掌握单元数组的生成和操作，掌握结构体的生成和操作。
2.掌握MATLAB 的脚本文件及其编辑和调试方法，掌握MATLAB 程序设计和开发流程，掌握MATLAB 的关系运算,逻辑运算及函数操作，掌握MATLAB 流程控制语句。
3.掌握基本符号运算，掌握符号函数的绘制，掌握符号函数微积分的运算，掌握符号方程的求解方法，掌握符号积分变换，了解MAPLE 函数的调用方法，了解符号计算器的使用。
4.了解MATLAB 的图形窗口，掌握MATLAB 基本二维图形，三维图形的绘制，以及图形的绘制，如柱形图，饼状图，掌握图形注释的添加和管理，了解三维图形的视点控制及颜色，光照控制5.了解Matlab 的图形对象及其属性，掌握MATLAB 图形对象属性的设置及其查询，掌握MATLAB 的图形对象句柄的访问及其操作。
2.实验内容（1） 编写一个脚本，查找给定字符串中指定字符出现的次数和位（2） 创建2x2单元数组，创建 2×2 单元数组，第 1、2 个元素为字符串，第三个元素为整型变量，第四个元素为双精（double ）类型，并将其用图形表示。
（3） 创建一个结构体，用于统计学生的情况，包括学生的姓名、学号、各科成绩等。
然后使用该结构体对一个班级的学生成绩进行管理，如计算总分、平均分、排列名次等。
（4） 在MATLAB 中使用一个循环确定：如果用户最初在一个银行帐户中存储$10000，并且在每年的年终再存储$10000(银行每年支付6%的利息)，那么账户上要积累$1000000要需要多长时间。
（5）设x 为符号变量，()421f x x x =++，()32458g x x x x =+++，试进行如下运算：（1）()()f x g x + （2）()()f x g x ⨯ （3）求()g x 的反函数（4）求g 以()f x 为自变量的复合函数（6）合并同类项22325325x x x x-++--（7）因式分解将 7798666 进行因数分解，分解为素数乘积的形式（8）绘制下列函数的图像()2sinf x x x=+，[]0,2π（9）编写程序，该程序在同一窗口中绘制函数在[]0,2π之间的正弦曲线和余弦曲线，步长为/10π，线宽为 4 个象素，正弦曲线设置为蓝色实线，余弦曲线颜色设置为红色虚线，两条曲线交点处，用红色星号标记。
matlab符号计算实验报告
1. 已知x=6,y=5,利用符号表达式求z =>> syms x >> z=(x+1)/(sqrt(x+3)-sqrt(y)); >> subs(z,x,5) ans =6/(8^(1/2)-y^(1/2)) >> subs(ans,6) ans = 15.83382. 分解因式。
（1）x y -44； >> syms x y >> factor(x^4-y^4) ans =(x-y)*(x+y)*(x^2+y^2)（2）x x x +++64212575151 >> syms x >> factor(125*x^6+75*x^4+15*x^2+1) ans =(5*x^2+1)^33. 化简表达式（1）sin cos cos sin ββββ-1212；>> syms x y >> f=sin(x).*cos(y)-cos(x).*sin(y)； >> sfy1=simple(f) 结果：sfy1 =sin(x-y)（2）x x x +++248321>> syms x >> f=(4*x^2+8*x+3)/(2*x+1);sfy1=simplify(f) sfy1 =2*x+34、求下列极限，将完成实验的程序写到文件sy1.m 中：（1） （2） （3） （4）（5） （1）>> syms x >> F1=atan(x)/(x); >> w=limit(F1) w =1（2）>> syms x F2=((1+x)/(1-x))^(1/x); >> w=limit(F2) w =exp(2)（3）>> syms x F3=(x.*log(1+x))/(sin(x^2)); >> w=limit(F3) w =1（4）>> syms x F4=atan(x)/(x); >> w=limit(F4,x,inf) w =0（5）>> syms x F5=(1/(1-x)-1/(1-x^3)); >> w=limit(F5,x,1) w =NaN5、求下列函数的导数，将完成实验的程序写到文件sy2.m 中：1、 >> x = sym('x'); >> y1=(cos(x))^3-cos(3*x); >> diff(y1)ans =-3*cos(x)^2*sin(x)+3*sin(3*x)2、 >> x = sym('x'); >> y2=x.*sin(x).*(log(x)); >> diff(y2)ans =sin(x)*log(x)+x*cos(x)*log(x)+sin(x)3、>> x = sym('x'); >> y3=(x.*exp(x)-1)/sin(x); >> diff(y3)ans =(exp(x)+x*exp(x))/sin(x)-(x*exp(x)-1)/sin(x)^2*cos(x)4、 x x x x F 1011lim 2⎪⎭⎫ ⎝⎛-+=→31115lim()11x F x x →=---20sin )1ln(lim 3x x x F x +=→x x F x arctan lim 10→=arctan 4lim x x F x →∞=xx y 3cos cos 13-=xx x y ln sin 2=x xe y xsin 13-=cos x y e x =>> x = sym('x');y=cos(x).*exp(x); >> diff(y) ans =-sin(x)*exp(x)+cos(x)*exp(x)5、 >> x = sym('x');y=x^2.*sin(x); >> diff(y) ans = 2*x*sin(x)+x^2*cos(x)6、求下列函数的积分1、syms x a b c;int(sin(a*x).*sin(b*x).*sin(c*x)) ans =-1/4/(c+a-b)*cos((c+a-b)*x)+1/4/(-c+a-b)*cos((-c+a-b)*x)+1/4/(c+a+b)*cos((c+a+b)*x)-1/4/(-c +a+b)*cos((-c+a+b)*x) 2、>> syms x ;int(x^5+x^3-sqrt(x)/4) ans =1/6*x^6+1/4*x^4-1/6*x^(3/2)3、>> syms x ;int(x.*exp(x)/(x+1)^2,x,0,1) ans =1/2*exp(1)-1 4、 >> syms x y;F=int(int('x/(1+x*y)',x,0,1),y,0,1) F =2*log(2)-15、 由曲面22y x z +=，1=z ，2=z 所围成 >> syms x y z;F=int(int(int('x^2+y^2',x,-inf,inf),y,-inf,inf),z,1,2) F =Inf7、求下列级数的和（1） （2）（1）>> syms n;symsum(2*n-1/2^n,1,Inf) ans =Inf（2）syms n;symsum(1/n*(2*n-1),1,Inf) ans =Inf8、将函数 展开成2-x 的幂级数>> syms x;mtaylor(1/(x^2+5*x-3),x-2) 2sin 2y x x=⎰cxdx bx ax sin sin sin dx x x x )4(35⎰-+⎰+102)1(dx x xe x ⎰⎰+D dxdy xy x 1]1,0[]1,0[⨯=D ⎰⎰⎰Vzdxdydz 11212n n n I ∞=-=∑211(21)n I n n ∞==+∑21()53f x x x =+-。
MATLAB实验报告(1-4)
信号与系统MATLAB第一次实验报告一、实验目的1.熟悉MATLAB软件并会简单的使用运算和简单二维图的绘制。
2.学会运用MATLAB表示常用连续时间信号的方法3.观察并熟悉一些信号的波形和特性。
4.学会运用MATLAB进行连续信号时移、反折和尺度变换。
5.学会运用MATLAB进行连续时间微分、积分运算。
6.学会运用MATLAB进行连续信号相加、相乘运算。
7.学会运用MATLAB进行连续信号的奇偶分解。
二、实验任务将实验书中的例题和解析看懂，并在MATLAB软件中练习例题，最终将作业完成。
三、实验内容1.MATLAB软件基本运算入门。
1). MATLAB软件的数值计算：算数运算向量运算：1.向量元素要用”[ ]”括起来，元素之间可用空格、逗号分隔生成行向量，用分号分隔生成列向量。
2.x=x0:step:xn.其中x0位初始值，step表示步长或者增量，xn为结束值。
矩阵运算：1.矩阵”[ ]”括起来；矩阵每一行的各个元素必须用”,”或者空格分开；矩阵的不同行之间必须用分号”;”或者ENTER分开。
2.矩阵的加法或者减法运算是将矩阵的对应元素分别进行加法或者减法的运算。
3.常用的点运算包括”.*”、”./”、”.\”、”.^”等等。
举例：计算一个函数并绘制出在对应区间上对应的值。
2).MATLAB软件的符号运算：定义符号变量的语句格式为”syms 变量名”2.MATLAB软件简单二维图形绘制1).函数y=f(x)关于变量x的曲线绘制用语：>>plot(x,y)2).输出多个图像表顺序：例如m和n表示在一个窗口中显示m行n列个图像，p表示第p个区域，表达为subplot(mnp)或者subplot(m,n,p)3).表示输出表格横轴纵轴表达范围：axis([xmax,xmin,ymax,ymin])4).标上横轴纵轴的字母：xlabel(‘x’),ylabel(‘y’)5).命名图像就在subplot写在同一行或者在下一个subplot前：title(‘……’)6).输出：grid on举例1：举例2：3.matlab程序流程控制1).for循环：for循环变量=初值:增量:终值循环体End2).while循环结构：while 逻辑表达式循环体End3).If分支：(单分支表达式)if 逻辑表达式程序模块End(多分支结构的语法格式)if 逻辑表达式1程序模块1Else if 逻辑表达式2程序模块2…else 程序模块nEnd4).switch分支结构Switch 表达式Case 常量1程序模块1Case 常量2程序模块2……Otherwise 程序模块nEnd4.典型信号的MATLAB表示1).实指数信号：y=k*exp(a*t)举例：2).正弦信号：y=k*sin(w*t+phi)3).复指数信号：举例：4).抽样信号5).矩形脉冲信号:y=square(t,DUTY) (width默认为1)6).三角波脉冲信号：y=tripuls(t,width,skew)(skew的取值在-1~+1之间，若skew取值为0则对称)周期三角波信号或锯齿波:Y=sawtooth(t,width)5.单位阶跃信号的MATLAB表示6.信号的时移、反折和尺度变换:Xl=fliplr(x)实现信号的反折7.连续时间信号的微分和积分运算1).连续时间信号的微分运算：语句格式：d iff(function,’variable’,n)Function:需要进行求导运算的函数，variable:求导运算的独立变量，n：求导阶数2).连续时间信号的积分运算：语句格式：int(function,’variable’,a,b)Function:被积函数variable:积分变量a:积分下限b:积分上限(a&b默认是不定积分)8.信号的相加与相乘运算9.信号的奇偶分解四、小结这一次实验让我能够教熟悉的使用这个软件，并且能够输入简单的语句并输出相应的结果和波形图，也在一定程度上巩固了c语言的一些语法。
matlab符号计算实验总结
matlab符号计算实验总结
在这个实验中，我们使用了MATLAB中的符号计算工具，通过创建符号变量和符号表达式来进行数学计算和推导。
这些符号计算工具能够处理代数运算、微积分、方程求解等复杂的数学问题。
在实验过程中，我们首先学习了如何创建符号变量。
通过使用
'sym'函数，我们可以将常规变量转换为符号变量，从而进行符号计算。
例如，我们可以定义一个符号变量x：x = sym('x')。
接下来，我们学习了如何使用符号变量进行代数运算。
通过将符号变量组合成符号表达式，我们可以进行加减乘除等代数运算。
例如，我们可以定义一个符号表达式y = x^2 + 2*x + 1，并对其进行简化或展开操作。
除了代数运算，我们还学习了如何进行微积分计算。
通过使用符号变量和符号表达式，我们可以对函数进行求导和积分操作。
例如，我们可以对一个符号表达式y = x^3求导，并得到其导数表达式。
在实验中，我们还学习了如何使用符号计算工具解方程。
通过使用'solve'函数，我们可以求解方程的根。
例如，我们可以解一个一元二次方程，找到其根的解析解。
通过这个实验，我深刻理解了符号计算在数学问题中的重要性。
它能够帮助我们更好地理解数学概念和定理，并能够进行复杂的数学推导和计算。
MATLAB中的符号计算工具提供了强大的功能和简便的操作，使得数学问题的解决变得更加高效和准确。
总的来说，这个实验让我对MATLAB中的符号计算有了更深入的
了解和掌握。
我相信在今后的学习和研究中，这些符号计算工具将对我有很大的帮助。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
int(f)
int(g,1,2)
diff(f,x,2)
D=f/g
limit(D,x,0)
用ezplot命令绘出 在[-4,4]之间的图像。
ezplot('x^2*sin(x)',[-4,4])
试分别生成5阶的单位阵、8阶均匀分布的随机矩阵。
eye(5)
rand(8)
生成列向量x=[1, 3, 5, 7, 9, … , 29]。
实验报告
实验名称：数学建模课程设计1
院 系 部:数理专业班级：计科1102
学生姓名：张野学 号:**********
指导教师：石玉英成 绩：
华北电力大学
华北电力大学实习报告书
学生姓名
张野
性别
男
成绩
实习名称
数学建模课程设计1
实习班级
计科1102
实习周数
1
带队教师
石玉英
实习单位
数理系
地点
数理系机房
1.指出下面的M1，M2，M3分别是什么，并上机验证。
E=D([8 4 10 13],[7 1 6 9 2])
求出矩阵E的最大元素。
max(max(E))
用ezplot3命令绘出单位球的图像
s(b),sin(a)*sin(b),cos(a),[0,2*pi,0,2*pi])
（1）建立M脚本文件求1到100的和。
>> factor(‘12345678901234567890’)错误
ans =
(2)*(3)^2*(5)*(101)*(3803)*(3607)*(27961)*(3541)
syms x;f=-x+2*x^2+3*x^3+4*x^4;g=5*x^5-6*x^2+9;
A=simple(f*g)
B=simple(f/g)
disp(s);
function f=legrende(n)
syms x;
if(n==0)
f=1;
else if(n==1)
f=x;
else if(n>1)
f=((2*n-1)*legrende(n-1)-(n-1)*legrende(n-2))/n;
end
end
end
function hanshu(x)
a=[1:2:29]
x=a’
生成以题目2中向量x的元素为对角线的矩阵A，并输出A的行数。
A=diag(x)
size(A,1)
生成一个与A同阶的正态分布的随机矩阵B。
B=randn(15)
生成由矩阵A与B向量乘法得到的矩阵C。
C=A.*B
生成由矩阵B与C矩阵乘积得到的矩阵D。
D=B*C
生成一个由矩阵D的第8、4、10、13行和第7、1、6、9、2列组成的子矩阵E。
>> A=[a+b*x,sin(t)+u; x*exp(-t),log(y)+v]
A =
[ a+b*x, sin(t)+u]
[ x*exp(-t), log(y)+v]
试指出下面两条命令的结果是否相同，如果不同，哪个是正确的？为什么？
>> factor(sym('12345678901234567890'))
5050
>> sum1(150)
11325
建立M函数文件作以下计算:
1)自然数n的阶乘;
2) n中取m的组合。
（1）function f=mymutiply(n)
if(n<=1)
f=1;
else
f=n*mymutiply(n-1);
end
（2）function assemblage1(n,m)
s=mymutiply(n)/mymutiply(n-m)/mymutiply(m);
（2）建立M函数文件求
1+2+3+..+n
并求n=50,100,150的值。
（1）clear;
s=0;
for i=1:100
s=s+i;
end
（2）
function sum1(n)
s=0;
while n>=1
s=s+n;
n=n-1;
end
disp(s);
>> sum1(50)
1275
>> sum1(100)
if(x<1)
disp(x)
else if (x>=1&x<10)
disp(2*x-1);
else
disp(3*x-11);
end
end
clear;
A=[1 4 6 8;2 3 7 0;-2 -1 5 7;1 1 -1 0];
B=[-3 -7 6 -5;-2 1 3 -1;0 0 1 3;2 6 -1 0];
>> c1=(a1+b1-a1)/b1;c1=0
>> a2=sym(a1); b2=sym(b1);
>> c2=(a2+b2-a2)/b2;c2=1
c1=0
c2=1
变量类型不一样
3.设A是一个符号矩阵(定义如下)，试指出findsym(A,1)的输出结果，并由此能得出什么结论？
>> symsa b t u v x y;
A(A==-1)=-2;
B(B<0)=1;
B(B==0)=A(B==0);
A B
End
用ezplot3命令绘出单位球的图像
syms a b;
ezsurf(sin(a)*cos(b),sin(a)*sin(b),cos(a),[0,2*pi,0,2*pi])
>> a=1; b=2; c=3;d=4;
>> M1=[a,b;c,d];
>> M2='[a,b;c,d]';
>> M3=sym('[a,b;c,d]');
M1 =
1 2
3 4
M2 =
[a,b;c,d]
M3 =
[ a, b]
[ c, d]
2.下面语句计算出来的c1，c2相等吗，为什么？上机验证。
>> a1=1e10; b1=1e-10;