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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助 matlab交点问题
201/2返回列表12下一页
查看: 2214  |  回复: 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

铁虫 (小有名气)
自己顶起
回复此楼
你不开心，没人替你开心的
2楼2012-04-27 14:26:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

ice17890

铁虫 (小有名气)
引用回帖:
3楼: Originally posted by baobiao007 at 2012-04-27 21:20:59:
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为存放所有根的数组

syms x y
[x,y]=solve('tan(15*pi/180)*(x-0.5)=y','-tan(61.01*pi/180)*(x+0.5)=y','x','y')
ans=[x,y]
ans =

[ -.37072704483246323273931369246846, -.23331060849079649874922716666090]
为什么和你的程序求出来的第一个值不一样呢？可以告诉我一下吗？谢谢，麻烦你了 ，你的第一个x =

   -0.6744
赞一下 回复此楼
你不开心，没人替你开心的
6楼2012-04-28 10:47:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色
引用回帖:
6楼: Originally posted by ice17890 at 2012-04-28 10:47:36:
syms x y
=solve('tan(15*pi/180)*(x-0.5)=y','-tan(61.01*pi/180)*(x+0.5)=y','x','y')
ans=
ans =


为什么和你的程序求出来的第一个值不一样呢？可以告诉我一下吗？谢谢，麻烦你了 ，你的第一个x =

...

我没注意你那个CD直线方程前面有个负号，改成下面的代码就可以了：
clc;
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=-0.5.*(z+1)./(1-z)% 求到的根
plot(x,A);
hold on;
y=tan(k).*(x+0.5);
plot(x,y,'LineWidth',4,'Color','r');
赞一下 回复此楼
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
7楼2012-04-28 11:56:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baobiao007

木虫 (职业作家)

中国特色
请把上楼中“y=tan(k).*(x+0.5)”前面也加上负号，刚才又忘了。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

ice17890
我同意叔本华的观点，人们投身艺术和科学领域的强烈愿望之一就是逃离痛苦、残酷和枯燥无味的现实生活，逃离自己飘忽不定的七情六欲的桎梏。--爱因斯坦
8楼2012-04-28 12:01:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice17890

铁虫 (小有名气)
引用回帖:
8楼: Originally posted by baobiao007 at 2012-04-28 12:01:19:
请把上楼中“y=tan(k).*(x+0.5)”前面也加上负号，刚才又忘了。

太谢谢你了，我再试一下，非常感谢，我再给你点金币。谢谢你！不胜感激
赞一下 回复此楼
你不开心，没人替你开心的
9楼2012-04-28 13:53:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ice17890

铁虫 (小有名气)
送鲜花一朵
引用回帖:
8楼: Originally posted by baobiao007 at 2012-04-28 12:01:19:
请把上楼中“y=tan(k).*(x+0.5)”前面也加上负号，刚才又忘了。

不好意思，我不知道如何支付给你金币，我送你朵花吧，还有一个问题，我要是想让他显示x,y的坐标值，就是一个X对应一个Y的显示出来，需要加一个什么样的语句呢？谢谢了，太感谢你了!
赞一下 回复此楼
你不开心，没人替你开心的
10楼2012-04-28 14:34:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ice17890 的主题更新
201/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 18/900 2026-08-25 14:50 by Weed114114
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +6 yuleib84 2026-08-25 7/350 2026-08-25 14:40 by 六两废铜
[基金申请] 什么时候开奖？ +12 CrisMessi 2026-08-18 13/650 2026-08-25 14:40 by ZJTJZ
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 2026国自然函评费到账 +19 羊腰板 2026-08-21 22/1100 2026-08-25 10:38 by popular289
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗？20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭