jgchen1973 发表于 2008-5-27 19:28

这个小程序的错误在哪?

fid1=fopen('Hk1.dat','a+');
for z=0:0.1:k;
    f=z;
    g=inline('cos(f*sin(x))');
    =quad(g,-pi,pi);
    %J0=1/(2*pi)*S;
    fprintf(fid1,'% f % f\n',f,S);
end
fclose(fid1)
出现的错误为:
??? Error using ==> inline.subsref
Not enough inputs to inline function.

Error in ==> quad at 62
y = f(x, varargin{:});

Error in ==> fhkl at 5
    =quad(g,-pi,pi);
我不明白到底第五行到底出现了啥错误,请求指点!
谢谢了

ch_j1985 发表于 2008-5-27 20:44

原帖由 jgchen1973 于 2008-5-27 19:28 发表 http://www.chinavib.com/forum/images/common/back.gif
fid1=fopen('Hk1.dat','a+');
for z=0:0.1:k;
    f=z;
    g=inline('cos(f*sin(x))');
    =quad(g,-pi,pi);
    %J0=1/(2*pi)*S;
    fprintf(fid1,'% f % f\n',f,S);
end
fclose(fid1)
出现的错误 ...
doc quad
个人觉得好像不可以用inline,可以试试下面的代码:
fid1=fopen('Hk1.dat','a+');
z=0:0.1:2;
for i=1:length(z)
    f=z(i);
    g=@(x)cos(f*sin(x));
    =quad(g,-pi,pi);
    %J0=1/(2*pi)*S;
    fprintf(fid1,'% f % f\n',f,S);
end
fclose(fid1)

sigma665 发表于 2008-5-27 20:51

回复 楼主 的帖子

带参数积分,搜索论坛精华贴

jgchen1973 发表于 2008-5-27 23:05

谢谢2楼和3楼

我已经参考已有的帖子,做出来了!谢谢关心!
我看2楼的好象比我的还简单,再试试.
页: [1]
查看完整版本: 这个小程序的错误在哪?