关于inline函数的使用问题
想把dfdx0换成数组dfdx1,再积分,怎么不行呢?clear
clc
syms x z dfdx0
cfun=inline('5e13*cosh(dfdx0*z./1e6)','z','dfdx0')
sfun=inline('5e13*sinh(dfdx0*z./1e6)','z','dfdx0')
f=(1-x.^2).^(1/2);
dfdx=diff(f);
x1=-1:0.05:1;
dfdx1=subs(dfdx,x,x1);
Cfun=subs(cfun,dfdx0,dfdx1)
Sfun=subs(sfun,dfdx0,dfdx1)
A=quad(cfun,z,0,1e-7)
B=quad(sfun,z,0,1e-7)
??? Error using ==> inline.subsref
Not enough inputs to inline function.
Error in ==> quad at 62
y = f(x, varargin{:});
Error in ==> yy at 17
A=quad(cfun,z,0,1e-7)
[ 本帖最后由 eight 于 2007-8-30 13:28 编辑 ] 你的程序问题太多,例如连最基本的quad用法都有问题.
------这种问题参考我的一些示例贴很容易解决.
另:最好将原问题一并贴出来,这样别人就很容易给你好建议.
回复 #2 xjzuo 的帖子
今天看了您的示例帖受益匪浅。后来，又把自己的程序改了一下，好了一些。不过还是有点问题。刚学不久，很多问题都不懂。请版主多多指教。clear;clc
s0=0:0.05:0.3;
syms l1 x z dpdx
n2=100; l=16.52;L=l*0.001; R=0.015; E=2.2637e11;
v0=0.03527432;vi=2.19e-8; ss=3e6;
ff=zeros(6,length(s0));
for i=1:6
W=i*1000;
f=[];
for j=1:length(s0)
s1=s0(j);
u2=n2*pi*R/30; u1=((2-s1).*u2)./(2+s1); u=(u1+u2)./2; uc=u2-u1;
hmin=(2.65*vi^0.54*(v0*u).^0.7*R^0.43*L^0.13/E^0.03)./W.^0.13;
hc=hmin*4/3;
b=(8*W*R/pi/L/E).^(0.5);
ph=2*W/L/b/pi;
p=ph*(1-x.^2/b^2).^(1/2);
dpdx=vpa(diff(p));
ch=strcat(num2str(ss),'*cosh(dpdx*z/',num2str(ss),')/',num2str(v0));
sh=strcat(num2str(ss),'*sinh(dpdx*z/',num2str(ss),')/',num2str(v0));
cH=vpa(subs(ch,'dpdx',dpdx));
sH=vpa(subs(sh,'dpdx',dpdx));
CH=int(cH,z,0,hc);
SH=int(sH,z,0,hc);
c=vpa((uc*CH-SH*(CH^2-SH^2+uc^2)^0.5)/(CH^2-SH^2));
q=subs(strcat('dpdx*',num2str(hc),'+',num2str(ss),'*sinh(c)'),'dpdx',dpdx);
q=vpa(subs(q,'c',c));
fun=inline(q);
F1=quadl(fun,(-b),b)
F2=int(F1,l1,(-L/2),L/2);
f(j)=double(F2/W);
end
ff(i,:)=f;
end
ff
plot(s0,ff,'-*') 看来你并没有理解我的程序本意,而且有滥用的嫌疑。
你现在的程序看上去既不简洁,也不美观,更糟糕的是你用的还是符号计算。所以基本上我不认为这种程序有什么意义。
另: l1的出现基本上与y没有什么区别-----你能说出F1是l1的函数吗？
-----请将问题及公式用word上传一下吧,否则很难猜测你究竟想算什么。
有时间的话，我会再写一个示例程序，以便大家再次碰到这种问题时有个参考。
回复 #4 xjzuo 的帖子
其实我自己也这么认为。只是现在太菜了。先能算出结果来，增强一下信心，然后再追求简洁和美观。麻烦您了。:@)[ 本帖最后由 twomao 于 2007-8-31 17:01 编辑 ] 编程技术实在不行，所以只好加注释了。
另：l1或y确实没有用，所以，我也删掉了。请版主批评指正。:@)
% ============================================================================
clear;clc
s0=0:0.05:0.3;
syms x z dpdx
n2=100; l=16.52; L=l*0.001; R=0.015; E=2.2637e11;
v0=0.03527432;vi=2.19e-8; ss=3e6;
ff=zeros(6,length(s0));
for i=1:6
W=i*1000;
f=[];
for j=1:length(s0)
s1=s0(j);
u2=n2*pi*R/30; u1=((2-s1).*u2)./(2+s1); u=(u1+u2)./2; uc=u2-u1;
hmin=(2.65*vi^0.54*(v0*u).^0.7*R^0.43*L^0.13/E^0.03)./W.^0.13;
hc=hmin*4/3;
b=(8*W*R/pi/L/E).^(0.5);
ph=2*W/L/b/pi;
p=ph*(1-x.^2/b^2).^(1/2);
dpdx=vpa(diff(p));
%以上准备工作
%以下积分求CH、SH、c
ch1=strcat(num2str(ss),'*cosh(dpdx*z/',num2str(ss),')/',num2str(v0));
sh1=strcat(num2str(ss),'*sinh(dpdx*z/',num2str(ss),')/',num2str(v0));
ch=subs(ch1,'dpdx',dpdx);%因为最终要对x积分，所以需要将dpdx关于x的表达式代入
sh=subs(sh1,'dpdx',dpdx);%因为最终要对x积分，所以需要将dpdx关于x的表达式代入
cH=inline(ch)
sH=inline(sh)
CH=quadl(cH,z,0,hc) %这一句有问题。这一句完成后，CH应该是x的函数。
SH=quadl(sH,z,0,hc)
c=(uc*CH-SH*(CH^2-SH^2+uc^2)^0.5)/(CH^2-SH^2)
%以下将CH、SH、c代入求被积函数fun
q=subs(strcat('dpdx*',num2str(hc),'+',num2str(ss),'*sinh(c)'),'dpdx',dpdx);
q=vpa(subs(q,'c',c));
fun=inline(q);
%以下积分求F2，进而求f
F1=quadl(fun,(-b),b)
F2=quadl(F1,(-L/2),L/2);
f(j)=double(F2/W);
end
ff(i,:)=f;
end
ff
plot(s0,ff,'-*')
% ================================================================================
[ 本帖最后由 twomao 于 2007-8-31 22:27 编辑 ] 问题还是没有讲清楚.
从你上传的Word 文档看, F的被积函数＂只是x的函数＂，怎么可能求二重积分呢？难道你还要对常数积一次分？
回复 #7 xjzuo 的帖子
是的,就是再对常数积分一次.版主,我在6楼把程序分段了,看起来可能会好一些.[ 本帖最后由 twomao 于 2007-9-1 10:03 编辑 ] 你的程序我没有兴趣仔细看, 明显地, CH=quadl(cH,z,0,hc),SH=quadl(sH,z,0,hc)中"z"要去掉.
有时间我会写一个数值计算的示例程序,以便有人碰到此类问题时有个参考.
回复 #9 xjzuo 的帖子
可是cH和sH含有x与z两个变量。去掉后,就不知道对哪个变量积分了。出现提示:
??? Error using ==> inline.feval
Not enough inputs to inline function.
这几天,太麻烦您了.我自己再看看吧.再次感谢几天来您的关注.
[ 本帖最后由 twomao 于 2007-9-1 10:51 编辑 ] 这个问题可以这样描述：
求一个定积分，被积函数为x，其中含有两个中间变量。
这两个中间变量有点特殊，它们都是x的函数，而且是对z求定积分而得到。
页:
[1]