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]。
实验报告
实验名称：数学建模课程设计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;