选修课选了和自己专业毫不相干的数模,求同仁帮忙做一下吧,我是真的不会做啊
你这些题我也做过,其实还是很简单的,第四和第五没看明白
% 第一题
syms k;
symsum(k^2,k,1,100)
% 第二题
A = [3, 2, -3; 1, 2, -1; 3, 1, -2];
b = [14; 2; 10];
[rank(A), rank([A, b])] % 秩相等,而且为3,说明有唯一解
X = inv(A)*b % 第一种方法
rref([A, b]) % 第二种方法
% 第三题
fun = @(x,y)[y(2);-2*y(2)-2*y(1)];
[x ,y]=ode45(fun,[0,5],[1;0]); % 这里没说x的范围,就选取0-5
plot(x,y(:,1))
% 第六题
% 第一小问
x = [1:10];
y = [3.2, 5.3, 8.9, 14.7, 24.3, 40.1, 66.2, 109.1, 180.0, 296.8];
scatter(x,y)
% 第二小问
xi = 1:0.2:10;
yi = interp1(x,y,xi); % 分段线性插值
plot(xi, yi)
hold on;
scatter(x,y,"filled")
%1
k=1:1:100;
an=k.^2;
S = sum(an);
disp(S);
%2
syms x y z
[X,Y,Z]=solve(3*x+2*y-3*z==14,x+2*y-z==2,3*x+y-2*z==10);
vpa(X)
vpa(Y)
vpa(Z)
%3
[Y2]=dsolve('D2y+2*Dy+2*y==0','y(0)==1','Dy(0)==0','x');
vpa(Y2)
xs=-10:0.01:10;
ys=exp(-1.0*xs).*cos(xs) + exp(-1.0*xs).*sin(xs);
figure(1);
plot(xs,ys,'o-');
grid on;
%6
x=[1,2,3,4,5,6,7,8,9,10];
y=[3.2,5.3,8.9,14.7,24.3,40.1,66.2,109.1,180.0,296.8];
%6.(1)
figure(2);
plot(x,y,'o');
grid on;
%6.(2)
xint=1:0.2:10;
yint = interp1(x,y,xint);
figure(3);
plot(x,y,'ro');
hold on;
plot(xint,yint,'-gp');