资源描述：
《matlab上机学习指导3》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、七代数方程数值解例1:求方程组4xye10的实根,x4yx08symsxy;s=solve(4*x-y+exp(x)/10-1,-x+4*y+xA2/8,x,y);结果为复数，但可验证x=0.2326,y=0.0565,是方程组的解。求方程的数值解用命令：fzero,fsolvefzero的应用格式为x=fzero(@fun,xO),fun是一元函数，输出x是函数fun的一个零点。x0是一个数时，输出x0附近的零点；x0是向量[a,b]时，输岀函数在区间[a,b]中的零点。=__之间的零点。函数名
2、起为my得x=0;得x=-2.2447；例2：求函数先定义函数functiony=my(x)%y=x*sin(xA2-x-1);再解方程x=fzero(@my,0)再解x=fzero(@my,-2)都不满足，画图，观察，解x=fzero(@my,-0.6)得x=-0.6180;x=fzero(@my,-1.6)得x=-1.5956;如果用x=fzero(@my,[-2,-0.1]),则岀现???Errorusing==>fzeroThefunctionvaluesattheintervalendp
3、ointsmustdifferinsign.所以仍要通过画图观察求解。fsolve可以解方程组，应用格式为[x,f,h]=fsolve(@fun,x0),输岀x是函数fun在xO附近的一个零点，xO为迭代初值；输出f为fun在x处的函数值，应接近于0；h大于0说明计算结果可靠，否则不可靠，需另寻他法，函数fzero也有类似输出。例3：解例2,[x,f,h]=fsolve(@fun,-1.6),例4：解例1,先定义函数functiony=you(x)y(1)=4*x(1)-x(2)+exp(x(1)
4、)/10-1;y(2)=-x(1)+4*x(2)+x(1)A2/8;取迭代初值x0=[0,0],解方程组[x,f,h]=fsolve(@you,xO),或[x,f,h]=fsolve(@you,[0,0])++=练习1:fzero和fsolve只能求实根，试用它们解=210XX,看看发生了什么，再上网查函数roots的用法，用它来解方程210xx,再用函数solve求解，比较有什么不同。练习2：两个椭圆可能有0到4个交点，求下列两个椭圆的所有交点的坐标:2(x-2)+(y-3+2x)2=5；22(
5、x_3)^(y/3)2=4先用ezplot画图观察交点的大概位置，再用fsolve解出交点八数值积分符号积分intr-一+1112dxx1方法一:l=intC1/(1+sqrt(1-xA2)y)pretty(l)方法二f=1/(1+sqrt(1-xA2));l=int(f)数值积分梯形法：抛物线法symsx;a,bl=trapz(x,y),x是积分区间向量，y是被积函数值；(Simpson法)：I二quad(@fun,a,b),fun是被积函数，用法与fzero是积分上下限；相同,高精度Lobat
6、to法：l=quadl(@fun,a,b),用法同quad例2：正态分布1=e122,信号处理椭圆周长T以222asint022bcostdt12X2dx为例:x=-1:0.5:1;%y=exp(-x.A2/2);l=trapz(x,y)x1=-1:0.05:1;%y1=exp(-x1.A2/2);plot(x,y,x1,y1)I1=trapz(x1,y1)%I2二int('exp(・xA2/2)〉“)vpa(l2)%步长为0.5,画图比较A-o.1到0/积分0.05积分符号积分步长改为%符号积分
7、的结果显示为数值，注意结果仍为符号变量练习：我国第一颗人造地球卫星近地点距离h=439km,远地点距离H=2384km,绕地球一周需114分钟。取地球长半轴R=6400km。由于地球位于卫星椭圆轨道的一个焦点上，根据近地点距离和远地点距离可分别计算出椭圆长半轴、椭圆半焦距、椭圆短半轴a=(h+H+2R)/2,c=(H-h)/2,b=Vg2-c2利用椭圆周长公式:L=2」,计算卫星轨道周长，并计算卫星运动0的速度与第一宇宙速度(7.9km/秒)比较。如果取地球长半轴为R=6378(km),卫星速度会
8、如何变化？与第一宇宙速度相比如何？九线性规划例1：某工厂生产A,B两种产品，所用原料均为甲、乙、丙三种：生产一件产品所需原料和所获利润以及库存丿方,料t冃况女uF所/j、：只A丿耳悄屮1公丿I)Q加俾乙13丨)A丿耳悄内1公丿1丿A钊側1儿丿7nnn厂口口HAnpO伫Q耳伫tuuu
