k;end;end;end;end;10*.是否任意3的倍数m可以写成两个素数p1、p2、p3的和(试着写出所有的m=p1+p2+p3的可能形式)?解:functiony=fg(n);n=input('请输入n的值:');ifmod(n,3);error('n不是3的倍数.请重新运行.')elseifn<6;error('n必须不小于6.')elseform=1:n;fork=m:n;forp=k:nif(isprime(m))&(isprime(k))&(isprime(p))&(m+k+p==n);disp([num2str(n),'=',num2str(m),'+',num2str(k),'+',num2str(p)]);break;end;end;end;end;end;四、数据处理与拟合初步:11.通过测量得到一组数据:t12345678910y4.8424.3623.7543.3683.1693.0383.0343.0163.0123.005分别采用y=c1+c2e^(-t)和y=d1+d2te^(-t)进行拟合,并画出拟合曲线进行对比。解:12.计算下列定积分:Z1:解:functionf=jifenl(x)f=exp(-2*x);[z1,n]=quad(@jifen1,0,2)Z2:解:x=0:0.01:2;z2=exp(2*x);trapz(x,z2)ans=26.8000Z3:解:t=-1:0.01:1;z3=x.^2-3*x+0.5;trapz(x,z3)ans=-2.333313.微分方程组当t=0时,x1(0)=1,x2(0)=-0.5,求微分方程t在[0,25]上的解,并画出相空间轨道图像。解:t=0:0.01:25;x=1/2+1/2*exp(-t);y=1/8+1/6*exp(-t)-19/24*exp(-4*t);plot(t,x,t,y)
全文预览
