资源描述：
《MATLAB第04次上机作业.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、个人收集整理勿做商业用途第四次上机作业时间：2014-5—12E—205(周一）晚上6：00～10：00l符号运算：1.求符号矩阵的行列式值和逆，所得结果应采用“子表达式置换"简洁化（subexpr）.symsa11a12a13a21a22a23a31a32a33A=[a11a12a13；a21a22a23;a31a32a33];B=det(A）；C=inv（A);[B1,s]=subexpr（B,’s’)［C1,r］=subexpr(C，'r')2.对函数，当为正实数,求.symskazf=a^k*z^（—k）；s=s
2、ymsum(f,k，0,inf）3.对，求。symsxksf=2*(（x-1)/（x+1））^(2＊k+1）/(2*k+1);s=symsum(f，k,0，inf）4.求的导数，，symstyy=abs（sin(t));个人收集整理勿做商业用途z=diff(y,t)z1=subs（z，t,-eps)z2=subs(z，x,pi/2）1.计算二重积分symsxyf=int（int（x^2+y^2，y,1,x^2)，1,2）vpa（f)2.求的Fourier变换，并画出时频谱symsAtwtao;f=A＊(1—abs（t）/
3、tao)*(heaviside(t+tao）-heaviside(t-tao));F=fourier(f，t,w）h=subs（F，A,2);h=subs（h,tao,2);x=linspace(-20，20,1000);s=subs(h,w，x）；plot（x,s）8．用solve求方程组的解。symsxyf=solve(x^2+y^2—1，x*y—2,x，y)x1=f.xx2=f.y9.求一阶微分方程的解symsabtxa=dsolve(’Dx=a＊t^2+b*t','x(0）=2'，'t'）个人收集整理勿做商业用途
4、10。求边值问题的解，并绘制解的图形S=dsolve（'Df=3*f+4*g'，’Dg=—4＊f+3*g'，'f（0）=0'，'g（0）=1’,’x')f1=S。fg1=S。gclf；holdonf11=ezplot（f1,[20,60]）；set(f11,'Color’,'k’)g11=ezplot（g1，[20,60］）；set(g11,’Color'，'r')gridonl“证明”“哥德巴赫猜想”：Ø1742年6月7日哥德巴赫写信给当时的大数学家欧拉，正式提出了以下的猜想：a。任何一个大于6的偶数都可以表示成两个素
5、数之和.b.任何一个大于9的奇数都可以表示成三个素数之和.Ø陈景润于1966年证明：任何充份大的偶数都是一个质数与一个自然数之和,而后者可表示为两个质数的乘积。Ø任务：用计算机证明106以内的数满足歌德巴赫猜想a(1）=2；a（2)=3；k=3;fori=1:10^6％求出10^6以内的所有素数，记录在数组a中forj=2:floor（sqrt（i）)ifmod（i，j）==0break;endifj==floor(sqrt(i))a（k)=i；k=k+1;个人收集整理勿做商业用途endendendflag=0；%标记对
6、于整数n，是否找到要求的俩个或三个素数。如果找到，赋值为1；否则赋值为0forn=1:10^6ifmod（n，2)==0＆&n〉6%大于6的偶数fori=1：length(a）ifflag==1break；endforj=i:length（a）ifn==a（i）+a（j）%找到了这个偶数所对应的俩个素数flag=1；break；endendendifflag==0disp('error！’);%如果没有找到要求的三个素数，输出：error！endendflag=0;ifmod（n,2）==1＆&n>9％大于9的奇数for
7、i=1：length(a）ifflag==1break;endforj=i：length（a)ifflag==1break;endfork=j：length(a)ifn==a(i）+a(j）+a（k）％找到了这个奇数所对应的三个素数flag=1；break;endendendendifflag==0disp(’error！’）％如果没有找到要求的三个素数，输出：error！endend个人收集整理勿做商业用途flag=0;end
