首页 正弦函数和三角函数的积分及matlab编程[资料]

正弦函数和三角函数的积分及matlab编程[资料]

举报
开通vip

正弦函数和三角函数的积分及matlab编程[资料]正弦函数和三角函数的积分及matlab编程[资料] 正弦函数和三角函数的积分及Matlab编程 求正弦函数y = sinx从0到π的积分 当x = 0时,积分为0,画出积分的函数曲线。 定积分的结果为 ππ Sxxx,,,,sindcos2,00 不定积分的结果为 IxxxC,,,,sindcos, 其中C是积分常量,由初始条件决定。当x = 0时,积分为I = 0,必有C = 1。结果为 I = -cosx + 1 根据积分的基本概念,将积分区域分为多份,用矩形法求曲线下的近似面积表示积分的近似值 ...

正弦函数和三角函数的积分及matlab编程[资料]
正弦函数和三角函数的积分及matlab编程[资料] 正弦函数和三角函数的积分及Matlab编程 求正弦函数y = sinx从0到π的积分 当x = 0时,积分为0,画出积分的函数曲线。 定积分的结果为 ππ Sxxx,,,,sindcos2,00 不定积分的结果为 IxxxC,,,,sindcos, 其中C是积分常量,由初始条件决定。当x = 0时,积分为I = 0,必有C = 1。结果为 I = -cosx + 1 根据积分的基本概念,将积分区域分为多份,用矩形法求曲线下的近似面积 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 示积分的近似值 n Sfxx,,(),i,1i 矩形法的函数是sum(f)。 用梯形法求曲线下的近似面积表示积分的近似值 n,11 Sfxfxx,,,[()()],ii,12i,0 梯形法的函数是trapz(f)。 用数值积分的函数是quad和quadl,常用使用格式是 S = quad(f,a,b) 其中,f表示被积函数,a表示积分的下限,b表示积分的下限。 用符号的函数是int,常用使用格式是 S = int(f,a,b) 程序如下 %正弦函数的积分 clear %清除变量 x=linspace(0,pi); %自变量向量 dx=x(2); %间隔 y=sin(x); %被积函数 s1=sum(y)*dx %矩形法积分 s2=trapz(y)*dx %梯形法积分 f=inline('sin(x)'); %被积的内线函数 s3=quad(f,0,pi) %数值定积分 s4=int('sin(x)',0,pi) %符号积分 sc1=cumsum(y)*dx; %矩形法累积积分(精度稍差) sc2=cumtrapz(y)*dx; %梯形法累积积分 figure %创建图形窗口 plot(x,-cos(x)+1,x,sc1,'.',x,sc2,'o') %画解析式和矩阵法以及梯形法积分曲线 s=int('sin(x)') %符号积分 sc3=subs(s,'x',x); %替换数值求符号积分的值 C=-sc3(1) %求积分常数 hold on %保持图像 plot(x,sc3+C,'c*') %画符号法积分曲线 grid on %加网格 fs=16; %字体大小 xlabel('\itx','FontSize',fs) %横坐标 ylabel('\intsin\itx\rmd\itx','FontSize',fs)%纵坐标 title('正弦函数的积分','FontSize',fs) %标题 legend('解析解','矩形法','梯形法','符号法')%图例 2.三角函数和指数的积分 [问题]求如下函数的积分 axy = esinbx 其中a = 0.5,b = 2。积分下限为0。画出积分的函数曲线。 [数学模型] 设 11axaxaxaxIbxxbxbxbbxx,,,, esindsinde{esinecosd},,,aa 11bbaxaxaxaxax ,,,,,{esincosde}{esin[ecosesind]}bxbxbxbxbbxx,,aaaa 因此不定积分为 1axIabxbbxC,,,e(sincos) 22ab, 当x = 0时,I应该为零,所以 bC, 22ab, 因此,从0开始的积分为 1axIabxbbxb,,,e(sincos) 22ab, [算法]可用积分的解析式直接画图,也可利用被积函数通过梯形求和指令cumtrapz求积分,数值积分指令quad求积分,还可通过符号积分int求解。 程序如下。 %数值积分和符号积分方法 clear %清除变量 a=0.5; %指数的常数 b=2; %正弦函数的常数 dx=0.1; %间隔 xm=6; %上限 x=0:dx:xm; %自变量向量 s1=(exp(a*x).*(-b*cos(b*x)+a*sin(b*x)))/(a^2+b^2);%积分的解析解 C=-s1(1); %求积分常数 y=exp(a*x).*sin(b*x); %被积函数 s2=cumtrapz(y)*dx; %梯形法积分 %s2=cumsum(y)*dx; %矩形法求积分 figure %创建图形窗口 plot(x,s1+C,x,s2,'.') %画积分曲线 s=['exp(',num2str(a),'*x).*sin(',num2str(b),'*x)'];%被积分函数字符串 f=inline(s); %化为内线函数 s3=0; %第1个值 for i=2:length(x) %按自变量循环 s3=[s3,quad(f,0,x(i))]; %连接积分 end %结束循环 s=int('exp(a*x)*sin(b*x)') %对x进行符号积分 ss=subs(s,{'a','b'},{a,b}); %替换常数 s4=subs(ss,'x',x); %替换向量 hold on %保持图像 plot(x,s3,'or',x,s4-s4(1),'ch','MarkerSize',16)%画数值积分和符号积分曲线 title(['\ity\rm=e^{',num2str(a),'}\it^x\rmsin',num2str(b),'\itx\rm的积分'],'FontSize',16)%标题 legend(' 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 法','梯形法','数值法','符号法',2)%加图例 zqy4.1图 zqy4.2图
本文档为【正弦函数和三角函数的积分及matlab编程[资料]】,请使用软件OFFICE或WPS软件打开。作品中的文字与图均可以修改和编辑, 图片更改请在作品中右键图片并更换,文字修改请直接点击文字进行修改,也可以新增和删除文档中的内容。
该文档来自用户分享,如有侵权行为请发邮件ishare@vip.sina.com联系网站客服,我们会及时删除。
[版权声明] 本站所有资料为用户分享产生,若发现您的权利被侵害,请联系客服邮件isharekefu@iask.cn,我们尽快处理。
本作品所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用。
网站提供的党政主题相关内容(国旗、国徽、党徽..)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
下载需要: 免费 已有0 人下载
最新资料
资料动态
专题动态
is_841159
暂无简介~
格式:doc
大小:24KB
软件:Word
页数:4
分类:生活休闲
上传时间:2017-12-07
浏览量:31