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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab解方程问题求助
121/2返回列表12下一页
查看: 1628  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

上太空说

新虫 (初入文坛)

[求助] matlab解方程问题求助 已有2人参与

clc;
clear all;
close all;
cl=1500;cs=3000;bs=1800;
pl=1000;ps=2500;H=50;h=25;f=30;w=2*pi*f;
kal=w/cl;kas=w/cs;kbs=w/bs;
syms x;
sln=sqrt(kal^2-x^2);ssn=sqrt(x^2-kas^2);lsn=sqrt(x^2-kbs^2);
kn=solve('tan(sln*H)=(ps/pl)*((x/kbs)^4)*(sln/ssn)*(4*ssn*ssn/(x^2)-((x^2+lsn)/x^2)^2)','x');
这是源程序，在2014a中运行不报错，但是得不出一个准确的结果，只是一个4*1的矩阵，请教大神问题出在什么地方，如何解出具体的结果，不胜感激！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-05-03 20:11:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
上太空说: 金币+50, ★★★★★最佳答案 2016-05-04 11:27:43
CODE:
function jie_fangcheng
clc;
clear all;
close all;
format long
xx1=[linspace(0.1048,0.1256,1000)];
xx2=[linspace(-0.1256,-0.1048,1000)];
yy1=obj(xx1);
a1=0.1*(max(xx1)-min(xx1));
b1=0.1*(max(yy1)-min(yy1));

xx2=[linspace(-0.1256,-0.1048,1000)];
yy2=obj(xx2);
a2=0.1*(max(xx2)-min(xx2));
b2=0.1*(max(yy2)-min(yy2));

figure(1)
subplot(2,1,1)
plot([min(xx1)-a1 max(xx1)+a1],[0 0],'r-',xx1,yy1,'b-');axis([min(xx1)-a1 max(xx1)+a1 min(yy1)-b1 max(yy1)+b1]);


subplot(2,1,2)
plot([min(xx2)-a2 max(xx2)+a2],[0 0],'r-',xx2,yy2,'b-');axis([min(xx2)-a2 max(xx2)+a2 min(yy2)-b2 max(yy2)+b2]);


[x fval]=fsolve(@obj,[ -0.105 -0.12 0.105 0.12  ])   

function f=obj(x)

cl=1500;cs=3000;bs=1800;
pl=1000;ps=2500;H=50;h=25;f=30;w=2*pi*f;
kal=w/cl;
kas=w/cs;
kbs=w/bs;

sln=sqrt(kal^2-x.^2);
ssn=sqrt(x.^2-kas^2);
lsn=sqrt(x.^2-kbs^2);
f=(ps/pl)*((x/kbs).^4).*(sln./ssn).*(4*ssn.*ssn./(x.^2)-((x.^2+lsn)./x.^2).^2)-tan(sln*H);

解x =

  -0.104895162970198  -0.121516622857776   0.104895162970198   0.121516622857776


方程残差fval =

  1.0e-010 *

  -0.000201505478969  -0.051869619710487  -0.000201505478969   0.749196260585450

matlab解方程问题求助
赞一下(3人) 回复此楼

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

上太空说 上太空说
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
6楼2016-05-04 11:04:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

Mr__Right

专家顾问 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
1. x全都是以2次形式出现的，简化的话，不如先求x^2;
2. 符号解可能是求不出来的，求关于 t=x^2 的非线性方程 F(t)=0的根的时候，先画图出来可以看到在(0.01,0.016)之间 t有多根
赞一下 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
2楼2016-05-03 20:32:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

上太空说

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by Mr__Right at 2016-05-03 20:32:38
1. x全都是以2次形式出现的，简化的话，不如先求x^2;
2. 符号解可能是求不出来的，求关于 t=x^2 的非线性方程 F(t)=0的根的时候，先画图出来可以看到在(0.01,0.016)之间 t有多根

多谢您啊，对于第2点求根要用什么函数呢？

发自小木虫IOS客户端
赞一下 回复此楼
3楼2016-05-03 23:21:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Mr__Right

专家顾问 (著名写手)
引用回帖:
3楼: Originally posted by 上太空说 at 2016-05-03 23:21:54
多谢您啊，对于第2点求根要用什么函数呢？
...

数值解，用 fsolve, 自己看文档和说明 ， 用doc fsolve
赞一下 回复此楼
文章乃身外之物,要多考虑编辑、审稿人和读者的感受。
4楼2016-05-04 07:28:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

上太空说

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by Mr__Right at 2016-05-04 07:28:36
数值解，用 fsolve, 自己看文档和说明 ， 用doc fsolve...

按照您的说法试过了，还是做不出来，您能否把程序发给我看看，不胜感激。

发自小木虫IOS客户端
赞一下 回复此楼
5楼2016-05-04 10:48:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

上太空说

新虫 (初入文坛)
送红花一朵
引用回帖:
6楼: Originally posted by 月只蓝 at 2016-05-04 11:04:14
function jie_fangcheng
clc;
clear all;
close all;
format long
xx1=;
xx2=;
yy1=obj(xx1);
a1=0.1*(max(xx1)-min(xx1));
b1=0.1*(max(yy1)-min(yy1));

xx2=;
yy2=obj(xx2);
a2=0.1*(max(xx2)-mi ...

感激涕零，多谢大神！

发自小木虫IOS客户端
赞一下 回复此楼
7楼2016-05-04 11:26:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

上太空说

新虫 (初入文坛)
送红花一朵
引用回帖:
6楼: Originally posted by 月只蓝 at 2016-05-04 11:04:14
function jie_fangcheng
clc;
clear all;
close all;
format long
xx1=;
xx2=;
yy1=obj(xx1);
a1=0.1*(max(xx1)-min(xx1));
b1=0.1*(max(yy1)-min(yy1));

xx2=;
yy2=obj(xx2);
a2=0.1*(max(xx2)-mi ...

您好，我还有一个问题，当f由30变为其他数时候，比如我想做一个f由1到100分别求解，做的时候有部分数值是求解不出来的，请问应该做怎么样的修改呢？

发自小木虫IOS客户端
赞一下 回复此楼
8楼2016-05-05 11:19:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
8楼: Originally posted by 上太空说 at 2016-05-05 11:19:08
您好，我还有一个问题，当f由30变为其他数时候，比如我想做一个f由1到100分别求解，做的时候有部分数值是求解不出来的，请问应该做怎么样的修改呢？
...

初值得重设。
根据画出来的图，选取y=0附近的x数值，作为初值代入。
赞一下 回复此楼

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

上太空说
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
9楼2016-05-05 22:14:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

上太空说

新虫 (初入文坛)
送红花一朵
引用回帖:
9楼: Originally posted by 月只蓝 at 2016-05-05 22:14:54
初值得重设。
根据画出来的图，选取y=0附近的x数值，作为初值代入。...

您的意思是xx1和xx2是要先做图再去定义的，做图过程在您的程序中没有体现是么？那您下面的两个图是什么意思呢？新手不大懂，希望赐教……

发自小木虫IOS客户端
赞一下 回复此楼
10楼2016-05-06 10:32:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 上太空说 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 bjm2e4yior 2026-07-18 3/150 2026-07-19 16:43 by 0u3sh5hhfi
[考博] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 4/200 2026-07-19 16:10 by 0u3sh5hhfi
[教师之家] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 6/300 2026-07-19 16:03 by 0u3sh5hhfi
[考研] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 6/300 2026-07-19 16:03 by 0u3sh5hhfi
[硕博家园] 售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急 +3 ufgqpyvpi6 2026-07-18 4/200 2026-07-19 15:54 by 0u3sh5hhfi
[基金申请] 国社科通讯评审 +5 潜水冠军123 2026-07-18 6/300 2026-07-19 12:28 by lfy8008
[基金申请] 面上的水很深 +16 布布和一二 2026-07-18 21/1050 2026-07-19 10:07 by Vivilian
[基金申请] 小木虫没落了，除了祈祷帖子，几乎看不到有价值的帖子 +14 botar 2026-07-14 15/750 2026-07-19 00:15 by beefly
[基金申请] 5个%2F，非固定段22个字母 +8 jaasxulu 2026-07-15 9/450 2026-07-18 17:40 by dxcgbm
[基金申请] 生命口C13面上会评时间 +7 luckinging 2026-07-14 11/550 2026-07-17 20:43 by 虫卜1988
[基金申请] 生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？ +6 碧水00 2026-07-16 9/450 2026-07-17 15:21 by jaasxulu
[论文投稿] 投英文期刊如何查重 10+5 不知道取什么名� 2026-07-15 6/300 2026-07-17 09:07 by bobvan
[论文投稿] 跨出版社商投稿 20+5 倪好520 2026-07-13 5/250 2026-07-16 17:24 by 无视军团
[基金申请] 生生命学部会评 +4 wuke100666 2026-07-15 4/200 2026-07-16 00:12 by 啦啦小妮子
[基金申请] H口，面上，时间戳7月7日。 +3 plumage5 2026-07-15 3/150 2026-07-15 21:06 by 6543yes
[基金申请] 会评 +5 无名者登山 2026-07-15 5/250 2026-07-15 18:15 by 杠就是你对
[基金申请] E口会评 +11 布布和一二 2026-07-13 13/650 2026-07-15 16:34 by Kittylucky
[基金申请] E口青C会评结束了吗？ 40+3 小小书虫c 2026-07-14 9/450 2026-07-15 16:04 by 安静的知了
[论文投稿] 农业机械学报投稿周期 +3 G0DLIN 2026-07-14 7/350 2026-07-15 09:26 by xs74101122
[基金申请] 没收到消息，再次陪跑 +3 天空星辰fly 2026-07-14 3/150 2026-07-14 19:14 by 6543yes
信息提示
关闭
请填处理意见
关闭