请教在循环中五次方程如何求啊。在线等哦

snowfish

我做的是 一个循环，比方说

for delta=[0.01:0.01:20];
  ep, b1,c1,a, b, c2 都是 delta的函数，所以每次循环都有一个确定的值。k是常数。
  现在的问题是，对于这样一个方程：
ep/sqrt(i2)+(b1*i2+c1)/(a*i2^2+b*i2+c2)-k=0
怎么在每次循环中给出Ii2的值呢
后面还有一些参数再用到这个i2的值。
最后要画的是i2的一个函数的图。f(i2）=a*i2+b*i2^5+c...,
其中，a,b,c也是delta的函数，每次循环有不同的值。

end

谢谢大家能给一个建议！:loveliness:

[ 本帖最后由 snowfish 于 2008-11-7 16:51 编辑 ]

snowfish

??? Error using ==> solve
Error, (in allvalues/rootseq) cannot evaluate with symbolic coefficients

rocwoods

help fzero
帮助文档里再搜索 anonymous function
或者nested function
把这两个东东相关文档仔细看下。你这个问题就解决了。
即fzero+匿名函数结构
或者fzero+nested function结构都能解决

[ 本帖最后由 rocwoods 于 2008-11-7 10:08 编辑 ]

snowfish

还是不会啊，呵呵，能否明示一下啊：:@L
k=0.01;
count_1=0;
for delta=[0.01:0.01:5];
x=delta;  
B1=3*x;
C1=2*x^2;
B2=3*x^4;
C2=5*x;
A=2*x;
B=x;
C=x;
ep=0.1*sqrt(x) ;

ep/sqrt(I2)+(B1*I2+C1)/(A*I2*I2+B2*I2+C2)-k=0 怎么求I2呢？？？

Y=I2
count_1=count_1+1;
Y0(count_1)=Y;
end

dingd

用1stOpt一试，不知对否？：
LoopConstant x=[0.01:0.1:5];
Constant k=0.01;
ConstStr b1=3*x,C1=2*x^2,B2=3*x^4,C2=5*x,A=2*x,B=x,C=x,ep=0.1*sqrt(x);
PlotLoopData I2;
Function ep/sqrt(I2)+(B1*I2+C1)/(A*I2*I2+B2*I2+C2)-k=0;

snowfish

k=0.01;
count_1=0;
for delta=[0.01:0.1:2];
x=delta;  
B1=3*x;
C1=2*x^2;
B2=3*x^4;
C2=5*x;
A=2*x;
B=x;
C=x;
ep=0.1*sqrt(x) ;

u=inline('ep/sqrt(x)+(B1*x+C1)/(A*x^2+B*x+C2)-k','x','ep','B1','C1','A','B','C2','k');
I2=fzero(u,150,[],ep,B1,C1,A,B,C2,k)

Y=I2
count_1=count_1+1;
Y0(count_1)=Y;
end
plot(delta,Y0)

不知道这样写可以吗？但是我运行之后每次都是显示i2是nan，不知道什么问题，请高手能详细的指点下。。:loveliness:

rocwoods

1. 
   
2. function snowfish
   
3. k=0.01;
   
4. count_1=1;
   
5. 
   
6.     function u = snowfish_nested(I2)
   
7.         B1=3*delta;
   
8.         C1=2*delta^2;
   
9.         B2=3*delta^4;
   
10.         C2=5*delta;
    
11.         A=2*delta;
    
12.         B=delta;
    
13.         C=delta;
    
14.         ep=0.1*sqrt(delta) ;
    
15.         u = ep/sqrt(I2)+(B1*I2+C1)/(A*I2*I2+B2*I2+C2)-k;
    
16.     end
    
17. 
    
18. Y = zeros(length(0.01:0.01:2),1);
    
19. for delta=0.01:0.01:2;
    
20. Y(count_1) = fzero(@snowfish_nested,300);
    
21. count_1=count_1+1;
    
22. end
    
23. plot(0.01:0.01:2,Y,'*-')
    
24. end
    

复制代码

保存，运行snowfish试试看

snowfish

:@L 虽然没能运行，但是还是非常感谢！！！:loveliness:

sigma665

ls的，你是什么版本，那个可以运行

snowfish

谢谢大家，我的问题解决了，用inline函数就可以了，可能我的变量取的不合适，才导致出现nan的情况。。:lol :lol :lol

snowfish

6.5的，是不是版本太低了啊？:@L

rocwoods

是比较低。虽然6.5很经典，毕竟是6年前的版本了。这几年，MATLAB推出很多新的结构和功能。比如说匿名函数和nested function就是其中代表。匿名函数的调用效率比内联函数高的多。
譬如以我的电脑为例:
CPU扣肉T8100，2G内存，MATLAB2008a：

1. 
   
2. clear
   
3. >> f = @(x) sin(x);
   
4. g = inline('sin(x)');
   
5. tic;for k = 1:100000
   
6. a = f(k);
   
7. end;toc
   
8. tic;for k = 1:100000
   
9. a = g(k);
   
10. end;toc
    
11. Elapsed time is 0.070580 seconds.
    
12. Elapsed time is 17.550051 seconds.

复制代码

分别构造匿名的sin(x)和内联的sin(x)函数，然后对其调用10万次。时间相差200多倍！当然，MATLAB的buit-in函数，像sin,cos,log,exp,+,-*-等等调用是最快的。

1. 
   
2. tic;for k = 1:100000
   
3. a = sin(k);
   
4. end;toc
   
5. Elapsed time is 0.007112 seconds.
   

复制代码

[ 本帖最后由 rocwoods 于 2008-11-8 11:06 编辑 ]

snowfish

:lol 看来我很落后啊，呵呵，换个7.0的试哈:lol