Matlab讨论区 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

cycr1234 发表于 2006-3-25 23:55

[求助]关于积分的问题

我定义一个分段信号 uA1=; uA2=-uA1; uA=; 以及时间 t=linspace(0,2*pi,360); 然后我想求这个信号在一个周期2pi内的有效值 u=1/2/pi*int(uA^2,-inf,inf) 出错的原因，应该是由于求积分的话uA必须是t的函数 但不知道应该怎么样去解决这个问题

cycr1234 发表于 2006-3-26 14:37

帮帮忙啊

yqchenlolo 发表于 2006-3-26 15:25

u=1/2/pi*int(uA^2,-inf,inf) 里uA是常数向量，不太理解你想干什么？？ 说清楚点

happy 发表于 2006-3-26 16:04

回复：（cycr1234）[求助]关于积分的问题

只是一个数值积分的问题，用quad函数来实现，你的整体思路有问题，呵呵，说得不好一点有点乱七八糟，下面是实现的程序 先写一个m文件，来实现你给的分段函数 function y = myfun(x) x1=(x/pi-floor(x/pi))*180; y(x1<=30)=0; y(x1>30&x1<150)=1; y(x1>=150)=0; 然后再命令窗口运行 Q = quad(@myfun,0,2*pi)/(2*pi);
[此贴子已经被作者于2006-3-26 16:04:24编辑过]

cycr1234 发表于 2006-3-26 19:08

谢谢大虾 知道了产生周期信号的方法了

cycr1234 发表于 2006-3-26 19:25

不过还有个问题 我的函数里要是有两个变量怎么办？ 还可以用 命令quad 还求积分吗？

cycr1234 发表于 2006-3-26 19:29

就比如： function y = myfun(x，a) x1=(x/pi-floor(x/pi))*180; y(x1<=a)=0; y(x1>a&x1<(180-2a))=1; y(x1>=(180-2a))=0; quad的解释里好象没说 int的应用里就有int(y,x,x1,x2) 来区分是对哪个变量求积分

happy 发表于 2006-3-26 19:31

回复：（cycr1234）[求助]关于积分的问题

二重积分？dblquad
[此贴子已经被作者于2006-3-26 19:37:22编辑过]

aspen 发表于 2006-3-26 19:36

回复：（cycr1234）就比如：function y = myfun(x，a)...

<DIV class=quote>以下是引用cycr1234在2006-3-26 19:29:57的发言： 就比如： function y = myfun(x，a) x1=(x/pi-floor(x/pi))*180; y(x1<=a)=0; y(x1>a&x1<(180-2a))=1; y(x1>=(180-2a))=0; quad的解释里好象没说 int的应用里就有int(y,x,x1,x2) 来区分是对哪个变量求积分 </DIV> quad是数值积分，不能存在系数a
[此贴子已经被作者于2006-3-26 19:36:53编辑过]

cycr1234 发表于 2006-3-27 21:50

happy 大哥 我发现我的程序还是得用int来求积分 因为有两个变量，t,a 就比如a=; t=linspace(0,2*pi,1000); y=myfun(t,a); 然后要求出在不同 a 的情况下，y 的有效值，即 y^2 对 t 的积分 并画出 y 的有效值和 a 的的变化曲线 不知道这样表述你清不清楚

happy 发表于 2006-3-28 08:37

回复：（cycr1234）[求助]关于积分的问题

a是可以传递给函数的 q = quad(fun,a,b,tol,trace,p1,p2,...) 多个a之可以通过循环，多次调用quad，每次调用只计算一个a

cycr1234 发表于 2006-3-28 12:27

对 q = quad(fun,a,b,tol,trace,p1,p2,...) 的用法不甚了解， 可否简单解释下参数a, b tol, trace, p1, p2 的含义？

cycr1234 发表于 2006-3-28 12:40

而且我最主要的迷惑是 quad的数值积分，到底如何区分是对t或a积分 可以举个例子吗

happy 发表于 2006-3-28 15:18

回复：（cycr1234）对q = quad(fun,a,b,tol,trace,p1,...

<DIV class=quote>以下是引用cycr1234在2006-3-28 12:27:39的发言： 对 q = quad(fun,a,b,tol,trace,p1,p2,...) 的用法不甚了解， 可否简单解释下参数a, b tol, trace, p1, p2 的含义？</DIV>
如果输入的Trace值不为零，电脑会画出积分点的分布图，如果 Trace其值为零，则不会画出积分点的分布图。 其他的参数应该很清楚

happy 发表于 2006-3-28 15:22

回复：（cycr1234）而且我最主要的迷惑是quad的数值积...

<DIV class=quote>以下是引用cycr1234在2006-3-28 12:40:14的发言： 而且我最主要的迷惑是 quad的数值积分，到底如何区分是对t或a积分 可以举个例子吗</DIV>
 function y = myfun(x,a) x1=(x/pi-floor(x/pi))*180; y(x1<=a)=0; y(x1>a&x1<180-2*a)=1; y(x1>=180-2*a)=0; 命令窗口 a=30; q = quad(@myfun,0,2*pi,1.0e-6,0,a)

页: [1] 2

声振论坛's Archiver

[求助]关于积分的问题

回复：（cycr1234）[求助]关于积分的问题

回复：（cycr1234）[求助]关于积分的问题

回复：（cycr1234）就比如：function y = myfun(x，a)...

回复：（cycr1234）[求助]关于积分的问题

回复：（cycr1234）对q = quad(fun,a,b,tol,trace,p1,...

回复：（cycr1234）而且我最主要的迷惑是quad的数值积...