声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1439|回复: 2

[编程技巧] 求积分中两个函数句柄相乘问题。

[复制链接]
发表于 2014-3-10 10:32 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 牛小贱 于 2014-3-10 10:36 编辑

程序如下:
  1. clear all;
  2. y1=@(x) exp(-x.^2)+exp(-(x-4).^2);
  3.   y2=@(x) conj((exp(-(x+t).^2)+exp(-(x-4+t).^2)).*exp(-i.*2.*pi.*192.31.*t));
  4. z=@(x) y1(x).*y2(x);
  5. i=quadgk(z,-inf,inf);
复制代码
错误提示:
??? Error using ==> plus
Matrix dimensions must agree.
Error in ==> @(x)conj((exp(-(x+t).^2)+exp(-(x-4+t).^2)).*exp(-i.*2.*pi.*192.31.*t))

Error in ==> @(x)y1(x).*y2(x)

Error in ==> quadgk>evalFun at 357
            fx = FUN(x);
Error in ==> quadgk>f4 at 412
        [y,too_close] = evalFun(tt);
Error in ==> quadgk>vadapt at 269
            [fx,too_close] = f(x);
Error in ==> quadgk at 233
    [q,errbnd] = vadapt(@f4,interval);


本帖被以下淘专辑推荐:

回复
分享到:

使用道具 举报

发表于 2014-3-10 15:13 | 显示全部楼层
先给你一个参考程序:
  1. y1 = @(f) f.^2;
  2. y2 = @(f) f+1;
  3. y3 = @(f) y1(f).*y2(f);
  4. fg = quadgk(y3,0,1)
复制代码
对于你的这个程序,y2=@(x) conj((exp(-(x+t).^2)+exp(-(x-4+t).^2)).*exp(-i.*2.*pi.*192.31.*t));在求共轭的时候是有问题的,导致后续的程序有问题。这会,我没有时间,先大体上说一下自己的想法。具体的,有时间我再试一下……



点评

赞成: 5.0
赞成: 5
  发表于 2014-3-12 11:03

评分

1

查看全部评分

回复 支持 1 反对 0

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-15 18:03 , Processed in 0.072907 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表