Labfans是一个针对大学生、工程师和科研工作者的技术社区。 论坛首页 | 联系我们(Contact Us)
MATLAB爱好者论坛-LabFans.com
返回   MATLAB爱好者论坛-LabFans.com > 工程数学软件 > MATLAB论坛
MATLAB论坛 一切MATLAB相关问题在此讨论。
回复
 
主题工具 显示模式
旧 2012-02-25, 11:16   #1
sky mat
初级会员
 
注册日期: 2012-02-25
年龄: 36
帖子: 1
声望力: 0
sky mat 正向着好的方向发展
默认 求各位大侠帮忙

clear;clc;
hold on
f1=inline('z2','z1','z2','z3','z4');
f2=inline('-1200*z1/(z1^2+z3^2)^1.5','z1','z2','z3','z4');
f3=inline('z4','z1','z2','z3','z4');
f4=inline('-1200*z3/(z1^2+z3^2)^1.5','z1','z2','z3','z4');
z10=30;z20=0;z30=0;z40=7.5;Ek=0;Ep=0;W=0;L=0;
fprintf('x=%.5f y=%.5f Vx=%.5f Vy=%.5f\n',z10,z30,z20,z40);
h=0.05;
plot(z10,z30,'k.','markersize',6);
for t=0:0.1:720
k11=h*f1(z10,z20,z30,z40);
k12=h*f2(z10,z20,z30,z40);
k13=h*f3(z10,z20,z30,z40);
k14=h*f4(z10,z20,z30,z40);
k21=h*f1(z10+k11*0.5,z20+k12*0.5,z30+k13*0.5,z40+k14*0.5);
k22=h*f2(z10+k11*0.5,z20+k12*0.5,z30+k13*0.5,z40+k14*0.5);
k23=h*f3(z10+k11*0.5,z20+k12*0.5,z30+k13*0.5,z40+k14*0.5);
k24=h*f4(z10+k11*0.5,z20+k12*0.5,z30+k13*0.5,z40+k14*0.5);
k31=h*f1(z10+k21*0.5,z20+k22*0.5,z30+k23*0.5,z40+k24*0.5);
k32=h*f2(z10+k21*0.5,z20+k22*0.5,z30+k23*0.5,z40+k24*0.5);
k33=h*f3(z10+k21*0.5,z20+k22*0.5,z30+k23*0.5,z40+k24*0.5);
k34=h*f4(z10+k21*0.5,z20+k22*0.5,z30+k23*0.5,z40+k24*0.5);
k41=h*f1(z10+k31*0.5,z20+k32*0.5,z30+k33*0.5,z40+k34*0.5);
k42=h*f2(z10+k31*0.5,z20+k32*0.5,z30+k33*0.5,z40+k34*0.5);
k43=h*f3(z10+k31*0.5,z20+k32*0.5,z30+k33*0.5,z40+k34*0.5);
k44=h*f4(z10+k31*0.5,z20+k32*0.5,z30+k33*0.5,z40+k34*0.5);
z10=z10+(k11+2*k21+2*k31+k41)/6;
z20=z20+(k12+2*k22+2*k32+k42)/6;
z30=z30+(k13+2*k23+2*k33+k43)/6;
z40=z40+(k14+2*k24+2*k34+k44)/6;
Ek=0.5*(z20^2+z40^2);
Ep=- 1200/(z10^2+z30^2)^0.5;
W=Ek+Ep;
L=z10*z40-z20*z30;
fprintf('x=%.5f y=%.5f Vx=%.5f Vy=%.5f Ek=%.5f Ep=%.5f W=%.5f L=%.5f\n',z10,z30,z20,z40,Ek,Ep,W,L);
plot(z10,z30,'k.','markersize',6);
end

请大家忙忙忙,以t为横坐标,绘出Ek,Ep,W,L的图像,谢谢
sky mat 当前离线   回复时引用此帖
旧 2012-02-25, 11:52   #2
peterjim
游客
 
帖子: n/a
默认 回复: 求各位大侠帮忙

clear;clc;
hold on
f1=inline('z2','z1','z2','z3','z4');
f2=inline('-1200*z1/(z1^2+z3^2)^1.5','z1','z2','z3','z4');
f3=inline('z4','z1','z2','z3','z4');
f4=inline('-1200*z3/(z1^2+z3^2)^1.5','z1','z2','z3','z4');
z10=30;z20=0;z30=0;z40=7.5;
% Ek=0;Ep=0;W=0;L=0;
% fprintf('x=%.5f y=%.5f Vx=%.5f Vy=%.5f\n',z10,z30,z20,z40);
h=0.05;
i=1;
plot(z10,z30,'k.','markersize',6);
for t=0:0.1:720
k11=h*f1(z10,z20,z30,z40);
k12=h*f2(z10,z20,z30,z40);
k13=h*f3(z10,z20,z30,z40);
k14=h*f4(z10,z20,z30,z40);
k21=h*f1(z10+k11*0.5,z20+k12*0.5,z30+k13*0.5,z40+k14*0.5);
k22=h*f2(z10+k11*0.5,z20+k12*0.5,z30+k13*0.5,z40+k14*0.5);
k23=h*f3(z10+k11*0.5,z20+k12*0.5,z30+k13*0.5,z40+k14*0.5);
k24=h*f4(z10+k11*0.5,z20+k12*0.5,z30+k13*0.5,z40+k14*0.5);
k31=h*f1(z10+k21*0.5,z20+k22*0.5,z30+k23*0.5,z40+k24*0.5);
k32=h*f2(z10+k21*0.5,z20+k22*0.5,z30+k23*0.5,z40+k24*0.5);
k33=h*f3(z10+k21*0.5,z20+k22*0.5,z30+k23*0.5,z40+k24*0.5);
k34=h*f4(z10+k21*0.5,z20+k22*0.5,z30+k23*0.5,z40+k24*0.5);
k41=h*f1(z10+k31*0.5,z20+k32*0.5,z30+k33*0.5,z40+k34*0.5);
k42=h*f2(z10+k31*0.5,z20+k32*0.5,z30+k33*0.5,z40+k34*0.5);
k43=h*f3(z10+k31*0.5,z20+k32*0.5,z30+k33*0.5,z40+k34*0.5);
k44=h*f4(z10+k31*0.5,z20+k32*0.5,z30+k33*0.5,z40+k34*0.5);
z10=z10+(k11+2*k21+2*k31+k41)/6;
z20=z20+(k12+2*k22+2*k32+k42)/6;
z30=z30+(k13+2*k23+2*k33+k43)/6;
z40=z40+(k14+2*k24+2*k34+k44)/6;
Ek(i)=0.5*(z20^2+z40^2);
Ep(i)=- 1200/(z10^2+z30^2)^0.5;
% W=Ek+Ep;
L(i)=z10*z40-z20*z30;
% fprintf('x=%.5f y=%.5f Vx=%.5f Vy=%.5f Ek=%.5f Ep=%.5f W=%.5f L=%.5f\n',z10,z30,z20,z40,Ek,Ep,W,L);
% plot(z10,z30,'k.','markersize',6);
i=i+1;
end
t=0:0.1:720;
W=Ek+Ep;
subplot(2,2,1)
plot(t,Ek)
subplot(2,2,2)
plot(t,Ep)
subplot(2,2,3)
plot(t,W)
subplot(2,2,4)
plot(t,L)
  回复时引用此帖
旧 2012-03-01, 10:38   #3
aiyaolei
普通会员
 
aiyaolei 的头像
 
注册日期: 2011-09-16
年龄: 39
帖子: 30
声望力: 14
aiyaolei 正向着好的方向发展
默认 回复: 求各位大侠帮忙

目前我在深造matlab软件, 现在技术不错, 有需要帮忙的可以联系我, 数学博士出生,QQ:1443736742
aiyaolei 当前离线   回复时引用此帖
旧 2012-03-01, 17:15   #4
peterjim
游客
 
帖子: n/a
帖子 回复: 求各位大侠帮忙

你这在推销自己呀,呵呵,不错,顶一下,
  回复时引用此帖
旧 2012-03-17, 22:33   #5
民大matlab求知者
初级会员
 
注册日期: 2011-12-28
帖子: 4
声望力: 0
民大matlab求知者 正向着好的方向发展
默认 回复: 求各位大侠帮忙

你好,学生在此向您求教一个问题,将车牌转换成HSV颜色空间后,怎样提取出蓝底白字车牌的蓝色和白色的H、S、V三个分量的范围???谢谢。可以直接回复我,或者QQ联系:534136830.小弟在此感激涕零。
民大matlab求知者 当前离线   回复时引用此帖
回复


发帖规则
不可以发表新主题
不可以发表回复
不可以上传附件
不可以编辑自己的帖子

启用 BB 代码
论坛启用 表情符号
论坛启用 [IMG] 代码
论坛禁用 HTML 代码



所有时间均为北京时间。现在的时间是 15:59


Powered by vBulletin
版权所有 ©2000 - 2025,Jelsoft Enterprises Ltd.