Matlab(大一数学实验)Ⅲ
Matlab(大一数学实验)Ⅲ

Matlab(大一数学实验)Ⅲ

所以大家记得改,由于学号决定答案,故不再给出答案

*如果需要源文件可以到前文下载

x=-50:0.3:50;
syms m
m=522;
y=x;
[X Y]=meshgrid(x,y);
Z=m*X.^2+Y.^4;
mesh(X,Y,Z);
A=[9,5;2,6];
x=[1,2];
x=x.';
syms n
[P,D]=eig(A);
An=P*(D^n)*inv(P);
x=An*x
A=[9,5;2,6];
B=A/10;
x=[1,2];
x=x.';
syms n
[P,D]=eig(B);
An=P*(D^n)*inv(P);
x=An*x
A=[9,5;2,6];
for i=1:20
x=[floor(rand()*15);floor(rand()*15)]
syms n
[P,D]=eig(A);
An=P*(D^n)*inv(P);
x=An*x;
x(2,1)/x(1,1)
end

*然后可用limit(f(x),inf)求极限

暂无(我没写到啊TAT)
暂无(我没写到啊TAT)
p0=[1/2;1/4;1/4];
A=[3/4,1/2,1/4;1/8,1/4,1/2;1/8,1/4,1/4];
for i=1:20
p0=A*p0
end

可能是这么写的(?

和下面题目重合度太高了,导致不知道朴素算法怎么写了……

p0=[1/2;1/4;1/4];
A=[3/4,1/2,1/4;1/8,1/4,1/2;1/8,1/4,1/4];
[P,D]=eig(A);
for i=1:20
P*(D^i)*inv(P)*p0
end
l=0;
m=[]; 
for a=1:2000
for c=a+1:1000
b=sqrt(c^2-a^2);
if (b==floor(b))&(b+2==c)
l=l+1;
m(:,l)=[a,b,c];
end
end
end
m
for k=1:7
l=0;
m=[]; 
k
for a=1:2000
for c=a+1:1000
b=sqrt(c^2-a^2);
if (b==floor(b))&(b+k==c)
l=l+1;
m(:,l)=[a,b,c];
end
end
end
m
end

{a,b,c}(c-b=k)={u,(k^2-u^2)/(2k),(3k^2+u^2)/(2*k)}

for k=1:200
l=0;
m=[]; 
k
for a=1:2:2000
for c=a+1:1000
b=sqrt(c^2-a^2);
if (b==floor(b))&(gcd(gcd(a,b),c)==1)&(b+k==c)
l=l+1;
m(:,l)=[a,b,c];
end
end
end
m
end

*请注意两题求取Pi近似值的题目,求出的值不是Pi,需要进行计算(当时被坑了好久……)

发表评论

您的电子邮箱地址不会被公开。 必填项已用*标注