luohubin

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 1562.8
散金: 50
紅花: 24
帖子: 604
在線: 132.9小時(shí)
蟲(chóng)號(hào): 429401
注冊(cè): 2007-08-04
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【求助】解方程組的錯(cuò)誤已有5人參與

大家好，我剛接觸matlab，想求兩條曲線的公切線方程，寫(xiě)了代碼如下
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，找不到錯(cuò)誤，請(qǐng)大家?guī)兔�，如果有更好的求公切線方程的方法希望能指點(diǎn)一下

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有16人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
材料工程專(zhuān)碩調(diào)劑已經(jīng)有6人回復(fù)
一志愿天大材料與化工（085600）總分338 已經(jīng)有4人回復(fù)
085700資源與環(huán)境308求調(diào)劑已經(jīng)有3人回復(fù)
求材料調(diào)劑已經(jīng)有8人回復(fù)
294求調(diào)劑材料與化工專(zhuān)碩已經(jīng)有5人回復(fù)
一志愿華中科技大學(xué)，080502，354分求調(diào)劑已經(jīng)有4人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有6人回復(fù)
085410人工智能專(zhuān)碩317求調(diào)劑（0854都可以）已經(jīng)有3人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

matlab求解非線性方程組已經(jīng)有16人回復(fù)
戈德斯坦課后習(xí)題，求解疑已經(jīng)有40人回復(fù)
matlab數(shù)值求解邊界條件微分方程組已經(jīng)有7人回復(fù)
求matlab求解一個(gè)方程組的問(wèn)題已經(jīng)有6人回復(fù)
修改稿已投，發(fā)現(xiàn)錯(cuò)誤，給主編發(fā)個(gè)郵件，是應(yīng)該重投還是等？謝謝，急！已經(jīng)有3人回復(fù)
最小二乘法解方程組系數(shù) 已經(jīng)有10人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題已經(jīng)有21人回復(fù)
【求助】常微分方程組求解中系數(shù)與某變量值關(guān)聯(lián)的問(wèn)題已經(jīng)有12人回復(fù)
【求助】用fortran求解大型線性方程組時(shí)出現(xiàn)的錯(cuò)誤【已解決】已經(jīng)有11人回復(fù)
【求助】求助：vb編程中用牛頓迭代解三次方程為什么只得到一個(gè)根？【已完成】已經(jīng)有15人回復(fù)
【求助】緊急求助有關(guān)MATHEMATICA的問(wèn)題已經(jīng)有9人回復(fù)

消愁無(wú)計(jì)對(duì)金樽

1樓 2010-07-29 14:54:19

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

ctgu_zheng

專(zhuān)家顧問(wèn) (著名寫(xiě)手)

專(zhuān)家經(jīng)驗(yàn): +133
應(yīng)助: 142 (高中生)
貴賓: 0.006
金幣: 3941.1
散金: 961
紅花: 115
沙發(fā): 1
帖子: 1337
在線: 831.7小時(shí)
蟲(chóng)號(hào): 751702
注冊(cè): 2009-04-18
性別: GG
專(zhuān)業(yè): 機(jī)械動(dòng)力學(xué)
管轄: 機(jī)械

★
nono2009(金幣+1):鼓勵(lì)應(yīng)助。 2010-07-29 17:22:26
luohubin(金幣+2):這些語(yǔ)句其實(shí)是要求出兩條曲線f1,f2的兩個(gè)公切點(diǎn)，結(jié)果是要得到兩個(gè)點(diǎn)的橫坐標(biāo)x1,x2. 至于solve，應(yīng)該沒(méi)問(wèn)題，用help solve可以找到范例，網(wǎng)上不少人也這樣用。還是感謝！ 2010-07-29 18:09:08

程序我看了會(huì)
公切線方程沒(méi)太明白到底是什么方程；
符號(hào)運(yùn)算 simplify好像用的不太多；
然后 solve命令也不對(duì) 如果是想求A1=0,A2=0的解的話
應(yīng)該是這樣的 [x1,x2]=solve('A1=0','A2=0');

贊一下(1人)

回復(fù)此樓

2樓2010-07-29 16:15:43

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

icyspace

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 794
帖子: 164
在線: 45.3小時(shí)
蟲(chóng)號(hào): 529026
注冊(cè): 2008-03-19
專(zhuān)業(yè): 金屬結(jié)構(gòu)材料

luohubin(金幣+2):啊？是嗎？我用comtan(1:1:10)怎么沒(méi)反應(yīng)？不是這樣輸入嗎？ 2010-07-30 21:48:27

我運(yùn)行了下，感覺(jué)挺正常的啊~我舍 p=1:1:10;可以得到結(jié)果的~

贊一下

回復(fù)此樓

3樓2010-07-30 14:44:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gy9991

adu886886:專(zhuān)業(yè)版塊請(qǐng)勿灌水，首次警告，下不為例！ 2010-07-30 17:02:30

4樓2010-07-30 15:17:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

icyspace

金蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 794
帖子: 164
在線: 45.3小時(shí)
蟲(chóng)號(hào): 529026
注冊(cè): 2008-03-19
專(zhuān)業(yè): 金屬結(jié)構(gòu)材料

★
robert2020(金幣+1):鼓勵(lì)應(yīng)助！歡迎常來(lái)仿真模擬版交流！ 2010-08-03 00:53:48
luohubin(金幣+1):謝謝 2010-08-04 09:12:51

試試這個(gè)

p=1:10;comtan(p);

[ Last edited by icyspace on 2010-7-31 at 13:23 ]

贊一下(1人)

回復(fù)此樓

5樓2010-07-31 13:21:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

robert2020

銀蟲(chóng) (著名寫(xiě)手)

仿真EPI: 2
應(yīng)助: 0 (幼兒園)
貴賓: 1.307
金幣: 292.3
散金: 9143
紅花: 19
沙發(fā): 6
帖子: 1950
在線: 620.1小時(shí)
蟲(chóng)號(hào): 861676
注冊(cè): 2009-10-02

luohubin(金幣+2): 2010-08-04 09:10:19

使用comtan(1:1:10)可以運(yùn)行（matlab7.0），把*.m文件最后一句“s=[x1,x2];”中的“;”去掉就會(huì)出結(jié)果

引用回帖:

Originally posted by luohubin at 2010-07-29 14:54:19:
大家好，我剛接觸matlab，想求兩條曲線的公切線方程，寫(xiě)了代碼如下
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* ...

贊一下

回復(fù)此樓

Truthfulness,Compassion,Tolerance,Beauty

6樓2010-08-03 00:53:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

luohubin

木蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 1 (幼兒園)
金幣: 1562.8
散金: 50
紅花: 24
帖子: 604
在線: 132.9小時(shí)
蟲(chóng)號(hào): 429401
注冊(cè): 2007-08-04
性別: GG
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

還是有問(wèn)題

引用回帖:

Originally posted by robert2020 at 2010-08-03 00:53:00:
使用comtan(1:1:10)可以運(yùn)行（matlab7.0），把*.m文件最后一句“s=[x1,x2];”中的“;”去掉就會(huì)出結(jié)果

版主，還是有問(wèn)題啊，我是 R2010版本，看你說(shuō)的做了，下面是出來(lái)的信息
Warning: Explicit solution could not be found.
> In solve at 81
In comtan at 19

s =

[ empty sym ]

ans =

[ empty sym ]
這是哪里出問(wèn)題了呢？

贊一下

回復(fù)此樓

消愁無(wú)計(jì)對(duì)金樽

7樓2010-08-04 09:11:58

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

robert2020

銀蟲(chóng) (著名寫(xiě)手)

仿真EPI: 2
應(yīng)助: 0 (幼兒園)
貴賓: 1.307
金幣: 292.3
散金: 9143
紅花: 19
沙發(fā): 6
帖子: 1950
在線: 620.1小時(shí)
蟲(chóng)號(hào): 861676
注冊(cè): 2009-10-02

luohubin(金幣+1):已經(jīng)很感謝了，本來(lái)想多給金幣，沒(méi)想到不能增加,不好意思 2010-08-07 16:19:07

樓主，不好意思，偶沒(méi)有裝2010版的，不太清楚是何問(wèn)題；
自己以前編過(guò)一個(gè)求半圓切線方程的程序，希望對(duì)樓主有所幫助。

引用回帖:

% 求兩個(gè)半圓的公切線
clc;clear all;close all
syms x
y1=sqrt(1-x.^2);
y2=2-sqrt(1-x.^2);
% 求導(dǎo)
y11=diff(y1,x);
y22=diff(y2,x);
% 求切線斜率 k0 及切點(diǎn)坐標(biāo)（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:

版主，還是有問(wèn)題啊，我是 R2010版本，看你說(shuō)的做了，下面是出來(lái)的信息
Warning: Explicit solution could not be found.
> In solve at 81
In comtan at 19

s =

[ empty sym ]

ans =
...

贊一下

回復(fù)此樓

Truthfulness,Compassion,Tolerance,Beauty

8樓2010-08-04 13:28:00

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 304求調(diào)劑 +6	司空. 2026-03-18	6/300	2026-03-18 23:03 by 星空星月
[考研] 085700資源與環(huán)境308求調(diào)劑 +3	墨墨漠 2026-03-18	3/150	2026-03-18 22:35 by bingxueer79
[考研] 一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑 +10	yangfz 2026-03-17	10/500	2026-03-18 20:14 by walc
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機(jī)用戶 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷 +10	曦熙兮 2026-03-15	10/500	2026-03-18 14:19 by 007_lilei
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3	Eember. 2026-03-17	9/450	2026-03-18 10:28 by Eember.
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛(ài)生物生物愛(ài)?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[基金申請(qǐng)] 今年的國(guó)基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤(rùn)色
[考研] 中科院材料273求調(diào)劑 +4	yzydy 2026-03-15	4/200	2026-03-16 15:59 by Gaodh_82
[考研] 機(jī)械專(zhuān)碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大本人本科雙一流 +4	自由的_飛翔 2026-03-13	5/250	2026-03-14 19:39 by hmn_wj
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 255求調(diào)劑 +3	李嘉慧， 2026-03-12	4/200	2026-03-14 16:58 by 有只貍奴
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 295求調(diào)劑 +3	小匕仔汁 2026-03-12	3/150	2026-03-13 15:17 by vgtyfty
[考研] 求調(diào)劑 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考博] 福州大學(xué)楊黃浩課題組招收2026年專(zhuān)業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655