大家帮我看看这个程序有什么问题

daiweisc

请大家帮我看看这个程序有什么问题？见附件
谢谢

[ 本帖最后由 eight 于 2008-3-6 22:16 编辑 ]

eight

原帖由 daiweisc 于 2008-3-6 21:43 发表

                               
登录/注册后可看大图

请大家帮我看看这个程序有什么问题？见附件
谢谢

最好不要让大伙通过下载运行你的程序才知道问题所在

sigma665

1. function swirlnozzle
   
2. % xFsolve.m
   
3. % 喷嘴计算，采用《简式旋流喷嘴的设计》文章中的方法，作者“姚增权．马进”
   
4. 
   
5. clear all
   
6. clc
   
7. 
   
8. 
   
9. Q = 0.6689/3600; % 流量，单位NM3/H
   
10. alpha = 60; % alpha单位为度
    
11. m = 1; % 旋流腔半径与排放口直径之比
    
12. a = 2*pi*alpha/360; % a将alpha单位转为弧度
    
13. Pl = 1e5; % 压力为100kPa
    
14. density = 1000; % 水的密度
    
15. miu = 1.0034e-6; % 水的运动粘度
    
16. x0 = 0.5;
    
17. theta = fsolve(@NonlinEqs,x0,[],a);
    
18. Cd = theta*sqrt(theta/(2-theta)); % 流量系数
    
19. Dout = sqrt(4*Q/(pi*Cd*sqrt(2*Pl/density))); % 排放口直径
    
20. R = m*Dout; % 旋流腔半径
    
21. K = sqrt(2)*(1-theta)/(theta*sqrt(theta));
    
22. Dp = sqrt(2*R*Dout/K);
    
23. Re = 4*Q/(pi*miu*Dp); % 雷诺数
    
24. s = 10^(25.8/(log10(Re))^2.58-2); % 摩擦系数
    
25. Ky = R*(Dout/2)/((Dp/2)^2+s/2*R*(R-Dout/2)); % 几何常数
    
26. theta1 = fsolve(@NonlinEqs1,x0,[],Ky);
    
27. Cd1 = theta1*sqrt(theta1/(2-theta1));
    
28. Dout1 = sqrt(4*Q/(pi*Cd1*sqrt(2*Pl/density)));
    
29. R1 = m*Dout1;
    
30. K1 = sqrt(2)*(1-theta1)/(theta1*sqrt(theta1));
    
31. Dp1 = sqrt(2*R1*Dout1/K1);
    
32. Dp2 = Dp1/sqrt(0.9); % 修正后的入口直径
    
33. disp(' 单位：mm')
    
34. fprintf('排放口直径=%.1f,修正后的入口直径=%.1f',Dout1*1000,Dp2*1000)
    
35. 
    
36. 
    
37. % ------------------------------------------------------------------
    
38. function f = NonlinEqs(x,a)
    
39. f(1) = tan(a/2)-(2*sqrt(2)*(1-x)/(sqrt(x)*(1+sqrt(1-x))));
    
40. % -------------------------------------------------------------------
    
41. function f = NonlinEqs1(x,Ky)
    
42. f(2) = Ky-sqrt(2)*(1-x)/(x^1.5)

复制代码

程序见上

dexy

不是能运行出来吗

eight

这种提问方式论坛不太提倡，啥都不说就摆个程序上来，以为大家是神？除非有同行路过，否则，祝楼主好运