数学与软件上机作业

发布 2022-09-05 10:29:28 阅读 4695

数学与软件matlab上机作业。

学院:数学与统计学院。

班级:071211班。

符号运算:

1. 求符号矩阵的行列式值和逆,所得结果应采用“子表达式置换”简洁化。

程序:a = sym ( a11 a12 a13 ; a21 a22 a23 ; a31 a32 a33] '

b = det(a);

c = inv(a);

bs,e] =subexpr( b,'e')

cs,d] =subexpr( c,'d' )

运行结果:bs =

a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31

e = empty sym ]

cs = d*(a22*a33 - a23*a32), d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)]

-d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), d*(a11*a23 - a13*a21)]

d*(a21*a32 - a22*a31), d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)]

d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)

2. 对函数,当为正实数,求。

程序:syms k a z

f = a^k * z^(-k);

s = symsum(f,k,0,inf)

运行结果:s =

piecewise([a = z, inf], a <>0 and z = 0, 0], a <>z and z <>0, limit((a^k*z)/(z^k*(a - z)),k = inf) -z/(a - z)])

3. 对,求。

程序:syms k x

f = 2*((x-1)/(x+1))^2*k+1)/(2*k+1)

s = symsum(f,k,0,inf)

运行结果:f =

2*((x - 1)/(x + 1))^2*k + 1))/2*k + 1)

s =piecewise([abs(x + 1)*(abs(x - 1) -abs(x + 1)) 0, 2*atanh((x - 1)/(x + 1))]

4. 求的导数,,

程序:syms x

f = abs(sin(x))

g = diff(f,x)

g1 = subs(g,x,-eps)

g2 = subs(g,x,pi/2)

运行结果:f =

abs(sin(x))

g =sign(sin(x))*cos(x)g1 =

g2 =6.1232e-017

5. 计算二重积分。

程序:syms x y

f = int(int((x^2+y^2),y,1,x^2),x,1,2)

vpa(f)

运行结果:f =

ans =

6. 求的fourier变换,并画出时的频谱。

程序;syms a t w tao;

ft = a*(1-abs(t)/tao)*(he**iside(t+tao)-he**iside(t-tao));

fw = fourier(ft,t,w)

h = subs(fw,a,2);

h = subs(h,tao,2);

x = linspace(-20,20,1000);

sn = subs(h,w,x);

plot(x,sn)

运行结果:fw =

a*((sign(tao)/w^2 + 1/w^2 + pi*dirac(-w, 1)*i - exp(tao*w*i)*sign(tao))/w^2 + tao*exp(tao*w*i)*sign(tao)*i)/w)/tao - 1/w^2 - sign(tao)/w^2 + pi*dirac(-w, 1)*i + 1/exp(tao*w*i))*sign(tao))/w^2 + tao*(1/exp(tao*w*i))*sign(tao)*i)/w)/tao + 1/exp(tao*w*i))*pi*dirac(-w) +i/w) -exp(tao*w*i)*(pi*dirac(-w) +i/w))

warning: imaginary parts of complex x and/or y arguments ignored

in workfive at 10

8.用solve求方程组的解。

程序:syms x y

s = solve(x^2+y^2-1,x*y-2,x,y)

x1 = x2 =

运行结果:s =

x: [4x1 sym]

y: [4x1 sym]

x1 =(1/2 + 15^(1/2)*i)/2)^(1/2)/2 - 1/2 + 15^(1/2)*i)/2)^(3/2)/2

- (1/2 + 15^(1/2)*i)/2)^(1/2)/2 + 1/2 + 15^(1/2)*i)/2)^(3/2)/2

(1/2 - 15^(1/2)*i)/2)^(1/2)/2 - 1/2 - 15^(1/2)*i)/2)^(3/2)/2

- (1/2 - 15^(1/2)*i)/2)^(1/2)/2 + 1/2 - 15^(1/2)*i)/2)^(3/2)/2

x2 =(1/2 + 15^(1/2)*i)/2)^(1/2)

-((15^(1/2)*i)/2 + 1/2)^(1/2)

(1/2 - 15^(1/2)*i)/2)^(1/2)

-(1/2 - 15^(1/2)*i)/2)^(1/2)

9. 求一阶微分方程的解。

程序:syms a b t x

s = dsolve('dx = a*t^2+b*t','x(0) =2','t')

运行程序:s =

t^2*(3*b + 2*a*t))/6 + 2

10. 求边值问题的解,并绘制解的图形。

程序:s = dsolve('df = 3*f + 4*g ',dg = 4*f+3*g','f(0) =0','g(0) =1' ,x')

f1 = g1 =

clf;hold on

f11 = ezplot(f1,[20,60]);

set(f11,'color','k')

g11 = ezplot(g1,[20,60]);

set(g11,'color','r')

grid on

运行结果:s =

f: [1x1 sym]

g: [1x1 sym]

f1 =sin(4*x)*exp(3*x)

g1 =cos(4*x)*exp(3*x)

“证明”“哥德**猜想”:

2023年6月7日哥德**写信给当时的大数学家欧拉,正式提出了以下的猜想:

a.任何一个大于 6的偶数都可以表示成两个素数之和。

b.任何一个大于9的奇数都可以表示成三个素数之和。

陈景润于2023年证明:任何充份大的偶数都是一个质数与一个自然数之和,而后者可表示为两个质数的乘积。

任务:用计算机证明106以内的数满足歌德**猜想。

程序:a(1) =2;

a(2) =3;

k = 3;

for i = 1:10^6%求出10^6以内的所有素数,记录在数组a中。

for j=2:floor(sqrt(i))

if mod(i,j) =0

break;

endif j ==floor(sqrt(i))

a(k) =i;

k = k + 1;

endend

endflag = 0; %标记对于整数n,是否找到要求的俩个或三个素数。如果找到,赋值为1;否则赋值为0

for n = 1:10^6

if mod(n,2)==0 &&n > 6%大于6的偶数。

for i = 1:length(a)

if flag==1

break;

endfor j = i:length(a)

if n==a(i) +a(j) %找到了这个偶数所对应的俩个素数。

flag = 1;

break;

endend

endif flag ==0

disp('error!')如果没有找到要求的三个素数,输出:error!

endend

flag = 0;

if mod(n,2)==1&& n > 9%大于9的奇数。

for i = 1:length(a)

if flag ==1

break;

endfor j = i:length(a)

if flag ==1

break;

endfor k = j:length(a)

if n==a(i) +a(j) +a(k) %找到了这个奇数所对应的三个素数。

flag = 1;

break;

endend

endend

if flag ==0

disp('error!')如果没有找到要求的三个素数,输出:error!

大学软件基础上机作业

商品信息管理系统,实现商品信息的输入和输出,添加和删除等操作 include 定义结构体 typedef struct goodstype int num char name 30 float price int stock goodstype goodstype goodlist 200 good...

物流软件上机实训作业

集装箱码头系统调度与优化。专业 姓名 学号 集装箱码头系统调度与优化。摘要。本文通过对闸口和堆场一体化的合理配置,以及场桥的调度规则,岸桥的动态调度规则进行简要分析,建立数学模型,结合相关 分析及 软件的应用,确定各子系统的资源调度配置情况,并对其结果进行优化,从而使得集装箱码头整体达到系统最优,提...

C语言作业与上机

作业与上机。1.以下程序的功能是计算并输出 计算精度要求 当时,停止累加。include main float s 1,t int n 1 t 1.0 2 while t 1e 6 printf f s 模仿上述程序,输入一个实数x,计算并输出 计算精度要求 当时,停止累加。2.以下程序的功能是运用...