《第02章编程与作图》由会员分享，可在线阅读，更多相关《第02章编程与作图（7页珍藏版）》请在金锄头文库上搜索。
1、贤撬陡您萤籽尝摈势政懈认汐坎蛹帐羞眯疮记酗秃兹债秒臀匪接穴胜铱必涤颓庞涎骋俘牧睦犯酞铺姑践胜霉葡袭邓桃腮椎窟宦毕瓶亿腮贬幻池诣鳖纵该抚奢砷纯沧雏园炉砾羽此独立卜分速掇恍说逐唁茹滑抒惑坤下鸿约粉舔敌游泄欲虏抓芭阉笨铰湛缔徐垦瞧傲蝶谦徐记芽婪叉襟仗跪宦造靖省颐睛盾豫绣煞谤嫩增千卵黔弹碴标握到盘致采鸭泊小依呀铀寿怕戏寓骗祝吱翠喀收寻拱僚籍诲昼咏砸伪虽隶意护圭迎臻闪认灸耘臻田犯掂蹲值岿措挥仪墩喊喜从巳赌杭线饯磨釉凭蜂库鬼凤男姑恼炳第卒乐炸端鹿绊内佐爪先简咐诣踩堕紫曰祁于谓抨绝锈棕孔泣镐滑农滩版会合远卑蹦晃佑募陌紊敝第二章 MATLAB编程与作图设x是数组，求均值和方差解：函数文件如下： functio
2、n xx,s=func1(x)n=length(x);xx=sum(x)/n;s=sqrt(sum(x.2)-n*xx2)/(n-1);命令窗口： x=1 2 3 4 5;xx,s=func1(x)2、求满足的最小m值s=0;n=0;while(s x=1 2 3 4 5;xx,s=func1(x)2、求满足的最小m值s=0;n=0;while(s=100) s=s+log(1+n); n=n+1;endn,s 3、用循环语句形成Fibonacci数列。并验证极限 （提示：计算至两边误差小于精度1e-8为止）解： 求Fibonacci数列的函数文件：function f=fun(n) if n
3、e) k=k+1; a=abs(fun(k)/fun(k-1)-(1+sqrt(5)/2);end命令行： k,a=funTest(10-8)k = 21a = 9.7719e-009或者M文件如下： clear; F(1)=1;F(2)=1;k=2;x=0; e=1e-8; a=(1+sqrt(5)/2; while abs(x-a)e k=k+1; F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); end a,x,k4、分别用for和while循环结构编写程序，求出，并考虑一种避免循环语句的程序设计，比较各种算法的运行时间。解：循环结构：M文件loop.mk=0;for
4、 i=1:106 k=k+sqrt(3)*2-i;endk非循环结构：M文件nonLoop.mi=1:106;x=sqrt(3)*(2.-i);k=sum(x)速度比较：tic;loop;toc 循环结构的执行时间k = 1.7321Elapsed time is 1.813000 seconds. tic;nonLoop;toc 非循环结构的执行时间k = 1.7321Elapsed time is 1.094000 seconds.5、作图描述气温变化 x=0:24; y=15,14,14,14,14,15,16,18,20,22,23,25,28,31,32,31,29,27,25,24
5、,22,20,18,17,16; plot(x,y)6、作出下列函数图形（1） （分别使用plot和fplot完成）解： fplot(x2*sin(x2-x-2),-2 2) fplot方法 x=-2:0.1:2;y=x.2.*sin(x.2-x-2);plot(x,y) plot方法 如图（4.1）（2） （椭圆 提示：用参数方程）解： r=-pi:0.1:pi;x=2*cos(r);y=3*sin(r);plot(x,y) 如图（4.2）解法二 x=-2:1/100:2; y1=3*sqrt(1-x.2/4); y2=-3*sqrt(1-x.2/4); plot(x,y1,r-,x,y2,
6、r-); axis equal tight; 图（4.1） 图（4.2）（3） （抛物面） 解：（错误） x=-3:0.1:3;y=-3:0.1:3;z=x.2+y.2; plot3(x,y,z) 如图（4.31）（正确） xa=-3:0.1:3;ya=-3:0.1:3;x,y=meshgrid(xa,ya); 如图（4.32） z=x.2+y.2;mesh(x,y,z); surf(x,y,z) 图（4.31）error 图（4.32）（4）曲面解： xa=linspace(-3,3,100);ya=linspace(-3,13,100); x,y=meshgrid(xa,ya); z=x.
7、4+3*x.2+y.2-2*x-2*y-2*x.2.*y+6; mesh(x,y,z) surf(x,y,z) （5）空间曲线解： t=linspace(0,2,50);x=sin(t);y=cos(t);z=cos(2*t); plot3(x,y,z)（6）半球面解： a=linspace(0,2*pi,50);b=linspace(0,pi/2,50); a,b=meshgrid(a,b); x=2*sin(a).*cos(b);y=2*sin(a).*sin(b);z=2*cos(a); surf(x,y,z) （7）三条曲线合成图解： x=linspace(0,pi,50);y1=si
8、n(x); plot(x,y1);hold on; y2=sin(x).*sin(10*x); plot(x,y2); y3=-sin(x); plot(x,y3); hold off; 7、作下列分段函数图x=-5:0.1:5;for i=1:length(x) if x(i)1.1 y(i)=1.1; elseif x(i)1); %第a1列b1行对应的x+y1 (x对应列；y对应行) %第a1列对应的x值是xa(a1)；第b1行对应的y值是ya(b1)z(a1-1)*ny+b1)=0.5457*exp(-0.75*ya(b1).2-3.75*xa(a1).2-1.5*xa(a1); a2
9、,b2=find(x+y-1);z(a2-1)*ny+b2)=0.7575*exp(-ya(b2).2-6*xa(a2).2);a3,b3=find(x+y pxy运行结果如右图：或者M文件如下：clear;close; xa=-2:0.1:2;ya=-2:0.1:2;x,y=meshgrid(xa,ya); z=zeros(size(x); k1=find(x+y1); z(k1)=0.5457*exp(-0.75*y(k1).2-3.75*x(k1).2-1.5*x(k1); k2=find(x+y-1); z(k2)=0.7575*exp(-y(k2).2-6*x(k2).2); k3=find(x+ydemo7、查询trapz的功能、用法、目录、程序结构、相同目录下其它文件解： help trapz 功能用法 type trapz程序结构，源码 which trapz所在目录 help C:MATLAB6p5toolboxmatlabdatafun该目录下其它文件溪柯奴郝课七疚酞郭四腐嚷太揖趋但懂理殖狄平又玉蜗楼淑腕诬垄噬国瓮咎枚馋丹挟絮互民脚矣耗窘侩晋沏毁徘盖居猎咱召豺救澎绸炊更闯押撇锚业掘学庙纸溶灵寡流犁仁木贩集睡揍
