应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助 matlab交点问题
52/1返回列表
查看: 1779  |  回复: 19
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
本帖产生 1 个 计算强帖 ，点击这里进行查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

ice17890

铁虫 (小有名气)

[求助] 求助 matlab交点问题

直线AB的方程为y=tan(θ)*(x-0.5)
直线CD的方程为y=-tan(61.01*pi/180)*(x+0.5)
求其中θ=15*pi/180：5*pi/180：85*pi/180变化时两直线的一系列交点，想问如何用程序实现，我能求出一个交点，就是当θ为一定值时，可是它变化了我就求不出来，请高手帮忙！
回复此楼

» 猜你喜欢
你不开心，没人替你开心的
1楼2012-04-27 13:45:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice17890

铁虫 (小有名气)
引用回帖:
14楼: Originally posted by baobiao007 at 2012-05-02 18:49:35:
你去论坛里“仿真模拟”版里找，有一些很不错的电子版书籍。

我也不知道是不是还可以问你别的问题，不过我还是问吧，嘿嘿，麻烦你了
Gf=1628.5
Gb=10184
Q=2700
q=500
L=60
a=500
xita=15*pi/180:5*pi/180:85*pi/180;
F=(Gf.*(xi+xj)+Gb.*(L/2.*cos(xita)+a+xj)+(Q+q+qs).*(L.*cos(xita)+a+xj))./bsyms a
l1=52.49
b=7.15
c=0.5
mf=1628.5
mh=80
lh=10.94
ml=1200
d=184
h=10
beta=61.01*pi/180
u=39.925
a=15*pi/180:5*pi/180:85*pi/180;
v=pi-57.48*pi/180-a
m=sqrt(2957.6932-2*14.23*52.49*cos(v))
y1=(m-u).*(l1.*sin(a)-(b-c)*tan(beta))./m
y=y1+(b-c)*tan(beta)
x=(y-(b-c)*tan(beta)).*(l1.*cos(a)+b)./(l1.*sin(a)-(b-c)*tan(beta))-b
xs=(-b+x)./2
xb=(x+l1.*cos(a))./2
mz=mh*lh
xi=(xs.*(m-u)*d+mz+xb.*ml)./mf
plot(a,xi)
得到，xj 就是上次你帮我弄的那个程序得到的，我该怎么样调用这两个值啊，谢谢了，非常感谢！
这里面的xi是通过
赞一下 回复此楼
你不开心，没人替你开心的
16楼2012-05-03 18:47:49
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

★
感谢参与，应助指数 +1
dbb627: 金币+1, 谢谢应助 2012-04-28 06:35:39
xita=15*pi/180:5*pi/180:85*pi/180;
k=61.01*pi/180.;
z=tan(xita)./tan(k);
x=1./(z-1) %x为存放所有根的数组
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
3楼2012-04-27 21:20:59
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色

【答案】应助回帖

★ ★ ★ ★ ★
ice17890: 金币+2, ★★★★★最佳答案, 很好 2012-04-28 09:37:14
cenwanglai: 金币+3, 感谢回复！ 2012-04-29 20:35:14
求根加画图：

xita=15*pi/180:5*pi/180:85*pi/180;
nxita=length(xita);
x=-30:0.01:30;
nx=length(x);
A=zeros(nxita,nx);
k=61.01*pi/180.;
z=tan(xita)./tan(k);
for ix=1:nx
    A(:,ix)=tan(xita).*(x(ix)-0.5);   
end
x1=1./(z-1)% 求到的根
plot(x,A);
hold on;
y=tan(k).*(x+0.5);
plot(x,y,'LineWidth',4,'Color','r');
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
4楼2012-04-27 21:51:35
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice17890

铁虫 (小有名气)
引用回帖:
4楼: Originally posted by baobiao007 at 2012-04-27 21:51:35:
求根加画图：

xita=15*pi/180:5*pi/180:85*pi/180;
nxita=length(xita);
x=-30:0.01:30;
nx=length(x);
A=zeros(nxita,nx);
k=61.01*pi/180.;
z=tan(xita)./tan(k);
for ix=1:nx
    A(:,ix)=tan(xita ...

谢谢你，我试试啊，不明白的再问你可以吗？非常感谢！
赞一下 回复此楼
你不开心，没人替你开心的
5楼2012-04-28 09:37:43
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭