- 在线时间
- 15 小时
- 最后登录
- 2017-2-4
- 注册时间
- 2011-8-3
- 听众数
- 4
- 收听数
- 0
- 能力
- 0 分
- 体力
- 107 点
- 威望
- 0 点
- 阅读权限
- 20
- 积分
- 58
- 相册
- 0
- 日志
- 0
- 记录
- 3
- 帖子
- 54
- 主题
- 2
- 精华
- 0
- 分享
- 0
- 好友
- 1
升级   55.79% TA的每日心情 | 郁闷 2012-2-12 10:05 |
---|
签到天数: 23 天 [LV.4]偶尔看看III
 |
5体力
这是网上下载的一段程序,我放入m文件之后debug->Run,一直显示??? Input argument "k" is undefined.
Error in ==> xielou at 2
R=sqrt((-4)*k*t*log((C/((1-r)*Q))*(4*pi*k*t)^(1.5))),如果去掉k,显示的则是t,那说明或许是因为变量没有定义。求高手指点~
function []=xielou(k,Q,C,t,r,ws,wa,wd)
R=sqrt((-4)*k*t*log((C/((1-r)*Q))*(4*pi*k*t)^(1.5)))
for i=0:24:t/3600
if (wd*(i/24))<=R
R=R-wd;
R0=R;
t=0:0.1:0.5*pi;
p=0:0.1:2*pi;
wsx=ws*cos(wa)*i;
wsy=ws*sin(wa)*i;
[theta,phi]=meshgrid(t,p);
x=R*sin(theta).*cos(phi)+wsy;
y=R*sin(theta).*sin(phi)+wsy;
z=R*cos(theta);
hold on;
surfc(x,y,z);
end
end
R0=R0;
view(20,70)
title('核泄漏有风向影响');
xlabel('x');
ylabel('y');
zlabel('z');
axis([-10,15,-10,25,0,10]);
x0 = ws*cos(wa)*1000000/3600;
y0 = ws*sin(wa)*1000000/3600;
R1 = sqrt((x0)^2+(y0)^2);
end |
zan
|