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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】解方程组的错误
81/1返回列表
查看: 1396  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

luohubin

木虫 (正式写手)

[交流] 【求助】解方程组的错误 已有5人参与

大家好，我刚接触matlab，想求两条曲线的公切线方程，写了代码如下
function s=comtan(p)
syms x1 x2;
a=p(1);
b=p(2);
c=p(3);
d=p(4);
e=p(5);
g=p(6);
h=p(7);
i=p(8);
j=p(9);
k=p(10);
f1=a+b*exp(c*x1^0.3333334)+d*exp(e*x1^0.3333334);
f2=g+h*exp(i*x2^0.3333334)+j*exp(k*x2^0.3333334);
df1=diff(f1,x1);
df2=diff(f2,x2);
A1=simplify((f2-f1)/(x2-x1)-df1);
A2=simplify(df2-df1);
[x1,x2]=solve(A1,A2);
s=[x1,x2];

但是总提示我Error: Unexpected MATLAB expression，找不到错误，请大家帮忙，如果有更好的求公切线方程的方法希望能指点一下
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
消愁无计对金樽
1楼 2010-07-29 14:54:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ctgu_zheng

专家顾问 (著名写手)
★
nono2009(金币+1):鼓励应助。 2010-07-29 17:22:26
luohubin(金币+2):这些语句其实是要求出两条曲线f1,f2的两个公切点，结果是要得到两个点的横坐标x1,x2. 至于solve，应该没问题，用help solve可以找到范例，网上不少人也这样用。还是感谢！ 2010-07-29 18:09:08
程序我看了会
公切线方程没太明白到底是什么方程；
符号运算 simplify好像用的不太多；
然后 solve命令也不对 如果是想求A1=0,A2=0的解的话
应该是这样的 [x1,x2]=solve('A1=0','A2=0');
赞一下(1人) 回复此楼
2楼2010-07-29 16:15:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

icyspace

金虫 (小有名气)
luohubin(金币+2):啊？是吗？我用comtan(1:1:10)怎么没反应？不是这样输入吗？ 2010-07-30 21:48:27
我运行了下，感觉挺正常的啊~我舍 p=1:1:10;可以得到结果的~
赞一下 回复此楼
3楼2010-07-30 14:44:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gy9991
adu886886:专业版块请勿灌水，首次警告，下不为例！ 2010-07-30 17:02:30
4楼2010-07-30 15:17:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

icyspace

金虫 (小有名气)
★
robert2020(金币+1):鼓励应助！欢迎常来仿真模拟版交流！ 2010-08-03 00:53:48
luohubin(金币+1):谢谢 2010-08-04 09:12:51
试试这个

p=1:10;comtan(p);

[ Last edited by icyspace on 2010-7-31 at 13:23 ]
赞一下(1人) 回复此楼
5楼2010-07-31 13:21:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

robert2020

银虫 (著名写手)
luohubin(金币+2): 2010-08-04 09:10:19
使用comtan(1:1:10)可以运行（matlab7.0），把*.m文件最后一句“s=[x1,x2];”中的“;”去掉就会出结果
引用回帖:
Originally posted by luohubin at 2010-07-29 14:54:19:
大家好，我刚接触matlab，想求两条曲线的公切线方程，写了代码如下
function s=comtan(p)
syms x1 x2;
a=p(1);
b=p(2);
c=p(3);
d=p(4);
e=p(5);
g=p(6);
h=p(7);
i=p(8);
j=p(9);
k=p(10);
f1=a+b* ...

赞一下 回复此楼
Truthfulness,Compassion,Tolerance,Beauty
6楼2010-08-03 00:53:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luohubin

木虫 (正式写手)

还是有问题

引用回帖:
Originally posted by robert2020 at 2010-08-03 00:53:00:
使用comtan(1:1:10)可以运行（matlab7.0），把*.m文件最后一句“s=[x1,x2];”中的“;”去掉就会出结果


版主，还是有问题啊，我是 R2010版本，看你说的做了，下面是出来的信息
Warning: Explicit solution could not be found.
> In solve at 81
  In comtan at 19

s =

[ empty sym ]


ans =

[ empty sym ]
这是哪里出问题了呢？
赞一下 回复此楼
消愁无计对金樽
7楼2010-08-04 09:11:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

robert2020

银虫 (著名写手)
luohubin(金币+1):已经很感谢了，本来想多给金币，没想到不能增加,不好意思 2010-08-07 16:19:07
楼主，不好意思，偶没有装2010版的，不太清楚是何问题；
自己以前编过一个求半圆切线方程的程序，希望对楼主有所帮助。
引用回帖:
% 求两个半圆的公切线
clc;clear all;close all
syms x
y1=sqrt(1-x.^2);
y2=2-sqrt(1-x.^2);
% 求导
y11=diff(y1,x);
y22=diff(y2,x);
% 求切线斜率 k0 及切点坐标（x0，y0）
f1=y2-y1;
f2=y22-y11;
x0=solve(f1);
k0=solve(f2);
y0=subs(y1,x0);
y=k0*x+y0; % 切线方程
x1=-1:0.001:1; % 定义自变量
y=subs(y,x1,'r'); % 切线
y1=subs(y1,x1); % 半圆 1
y2=subs(y2,x1); % 半圆 2
% 绘图
plot(x1,y)
hold on
plot(x1,y1)
plot(x1,y2)
axis([-1 1 0 2])

引用回帖:
Originally posted by luohubin at 2010-08-04 09:11:58:

版主，还是有问题啊，我是 R2010版本，看你说的做了，下面是出来的信息
Warning: Explicit solution could not be found.
> In solve at 81
  In comtan at 19

s =

[ empty sym ]


ans =
...

赞一下 回复此楼
Truthfulness,Compassion,Tolerance,Beauty
8楼2010-08-04 13:28:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 luohubin 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 7/350 2026-08-21 11:45 by 6543yes
[基金申请] 时间戳又变了8-15 +15 archvillain 2026-08-15 27/1350 2026-08-21 11:32 by tim76
[基金申请] 今天基金会出结果吗？20260819 +15 kkkl_v 2026-08-19 16/800 2026-08-21 11:22 by 365372687
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
[基金申请] 只有每年这种时候来逛逛小木虫 +22 yaoyewhu2008 2026-08-20 24/1200 2026-08-21 07:44 by DanielDDDD
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +15 医学老男孩 2026-08-20 17/850 2026-08-20 20:35 by CaoLG
[基金申请] 时间戳变了，能看出什么问题？ +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息，中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[教师之家] 为什么余额宝的年化利率越来越低？主要原因有哪些？ +5 瞬息宇宙 2026-08-15 5/250 2026-08-19 20:23 by super2002521
[基金申请] 今天放榜吗？ +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖？ +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[基金申请] 明天放榜？ +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
关闭
请填处理意见
关闭