MATLAB符号运算实验报告.doc

实验报告

实验名称：数学建模课程设计1

院系部:数理专业班级：计科1102
学生姓名：张野学号:1111080223
指导教师：石玉英成绩：

华北电力大学

华北电力大学实习报告书
学生姓名张野性别男成绩实习名称数学建模课程设计1实习班级计科1102实习周数1带队教师石玉英实习单位数理系地点数理系机房1.指出下面的M1，M2，M3分别是什么，并上机验证。
>>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=
12
34
M2=
[a,b;c,d]
M3=
[a,b]
[c,d]
2.下面语句计算出来的c1，c2相等吗，为什么？上机验证。
>>a1=1e10;b1=1e-10;
>>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)的输出结果，并由此能得出什么结论？
>>symsabtuvxy;
>>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'))
>>factor(‘12345678901234567890’)错误
ans=
(2)*(3)^2*(5)*(101)*(3803)*(3607)*(27961)*(3541)

symsx;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)
C=simple(f\g)
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]。
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。
E=D([841013],[71692])
求出矩阵E的最大元素。
max(max(E))
用ezplot3命令绘出单位球的图像
symsab;
ezsurf(sin(a)*cos(b),sin(a)*sin(b),cos(a),[0,2*pi,0,2*pi])
（1）建立M脚本文件求1到100的和。
（2）建立M函数文件求
1+2+3+..+n
并求n=50,100,150的值。
（1）clear;
s=0;
fori=1:100
s=s+i;
end
（2）
functionsum1(n)
s=0;
whilen>=1
s=s+n;
n=n-1;
end
disp(s);
>>sum1(50)
1275
>>sum1(100)
5050
>>sum1(150)
11325
建立M函数文件作以下计算:
1)自然数n的阶乘;
2)n中取m的组合。
（1）functionf=mymutiply(n)
if(n<=1)
f=1;
else
f=n*mymutiply(n-1);
end
（2）functionassemblage1(n,m)
s=mymutiply(n)/mymutiply(n-m)/mymutiply(m);
disp(s);

functionf=legrende(n)
symsx;
if(n==0)
f=1;
elseif(n==1)
f=x;
elseif(n>1)
f=((2*n-1)*legrende(n-1)-(n-1)*legrende(n-2))/n;
end
end
end

functionhanshu(x)
if(x<1)
disp(x)
elseif(x>=1&x<10)
disp(2*x-1);
else
disp(3*x-11);
end
end

clear;
A=[1468;2370;-2-