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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » 用matlab求解方程式
111/2返回列表12下一页
查看: 3473  |  回复: 10
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

钟声入翠微

铜虫 (小有名气)

[求助] 用matlab求解方程式 已有1人参与

假设a和b已知，a,b,x,y>0,求解x和y
a=[1+(0.0172*x)^2]^(-y/2)
b=y*arctan(0.0172*x)
希望可以提供程序代码，或者解！非常感谢！
回复此楼

» 猜你喜欢
1楼 2021-07-05 16:17:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
钟声入翠微: 金币+40, ★★★★★最佳答案 2021-07-08 21:19:15
引用回帖:
9楼: Originally posted by 钟声入翠微 at 2021-07-08 14:23:57
谢谢！（30000,0.05）这个初始值是自己根据实际情况设置，还是随意设置的呢？这次我运行出来了，但是感觉结果不是我想象中的那样。
...

这个需要尝试，可以用 fsolve 函数随机赋值求解。

结果不理想就要考虑方程是否正确，或者a,b取值是否合理。

function example_0708
clear;
a=0.709133065;
b=0.08555;
x0 = [1,1];
x = fsolve(@root2d,x0);
fprintf('\tx=%.8f\n',x(1));
fprintf('\ty=%.8f\n',x(2));

function F = root2d(x)
F(1) = (1+(0.0172*x(1))^2)^(-x(2)/2)-a;
F(2) = x(2)*atan(0.0172*x(1))-b;
end

end
赞一下(1人) 回复此楼

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

钟声入翠微
数值计算
10楼2021-07-08 16:30:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

zf015222
。

发自小木虫Android客户端
2楼2021-07-05 16:42:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

感谢参与，应助指数 +1
这个尽量不要做伸手党，求解符号方程组，看看solve函数用法

发自小木虫Android客户端
赞一下 回复此楼

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

钟声入翠微
数值计算
3楼2021-07-05 19:49:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
送红花一朵
引用回帖:
3楼: Originally posted by 独孤神宇 at 2021-07-05 19:49:57
这个尽量不要做伸手党，求解符号方程组，看看solve函数用法

我是去百度了solve用法，但是运行的时候一直显示错误，
我最先给a，b赋值进行计算，但是一直显示运算符错误
syms x y;
>>[x ,y]=solve(' 0.709133=[1+(0.0172*x)^2]^(-y/2)','0.08555=y*arctan(0.0172*x)');
然后使用了未赋值
clc,clear
syms a b y x
sola=solve(a==(1+(0.0172*x)^2)^(-y/2),x)   %待求解的变量是x
sol=solve(b==y*atan(0.0172*x),y)  %待求解的变量是y
还是运算不出想要的东西
因为只会一点点matlab，就一直解决不了问题。
赞一下 回复此楼
4楼2021-07-06 09:02:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
4楼: Originally posted by 钟声入翠微 at 2021-07-06 09:02:45
我是去百度了solve用法，但是运行的时候一直显示错误，
我最先给a，b赋值进行计算，但是一直显示运算符错误
syms x y;
>>=solve(' 0.709133=^(-y/2)','0.08555=y*arctan(0.0172*x)');
然后使用 ...

clear;
syms a b
syms x y positive
eqns=[(1+(0.0172*x)^2)^(-y/2)-a==0, y*atan(0.0172*x)-b==0];
[X,Y]=solve(eqns,[x y])

这个可能是方程有点问题吧。 a  b  的具体数值是多少
赞一下 回复此楼

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

钟声入翠微
数值计算
5楼2021-07-06 12:51:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
引用回帖:
5楼: Originally posted by 独孤神宇 at 2021-07-06 12:51:49
clear;
syms a b
syms x y positive
eqns=;
=solve(eqns,)

这个可能是方程有点问题吧。 a  b  的具体数值是多少...

0.709133065        0.08555
0.245947581        0.23527
0.679383641        0.63827
0.32593606        0.51032
0.208709677        0.53469

a和b的值如上，比较小，我也不清楚能不能解出来。
赞一下 回复此楼
6楼2021-07-06 14:42:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
送红花一朵
引用回帖:
5楼: Originally posted by 独孤神宇 at 2021-07-06 12:51:49
clear;
syms a b
syms x y positive
eqns=;
=solve(eqns,)

这个可能是方程有点问题吧。 a  b  的具体数值是多少...

我自己又代入了一下，得到这样的结果，是什么意思呢？

Warning: The solutions are parameterized by the symbols: z, z1. To include parameters and conditions in the solution, specify the
'ReturnConditions' option.
> In solve>warnIfParams (line 510)
  In solve (line 367)
Warning: The solutions are valid under the following conditions: 1/((1849*z^2)/6250000 + 1)^(z1/2) -
6387302814580195/9007199254740992 == 0 & z1*atan((43*z)/2500) - 1711/20000 == 0 & 0 < z & 0 < z1. To include parameters and
conditions in the solution, specify the 'ReturnConditions' option.
> In solve>warnIfParams (line 517)
  In solve (line 367)

X =

z


Y =

z1
赞一下 回复此楼
7楼2021-07-06 14:53:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 钟声入翠微 at 2021-07-06 14:53:48
我自己又代入了一下，得到这样的结果，是什么意思呢？

Warning: The solutions are parameterized by the symbols: z, z1. To include parameters and conditions in the solution, specify the
'ReturnCondit ...

初始值变化很大，用 vpasolve 函数赋初始值

clear;
a=0.709133065;b=0.08555;
syms x y positive
eqns=[(1+(0.0172*x)^2)^(-y/2)-a==0, y*atan(0.0172*x)-b==0];
[X,Y]=vpasolve(eqns,[x y],[30000;0.05])
赞一下 回复此楼

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

钟声入翠微
数值计算
8楼2021-07-06 17:42:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

钟声入翠微

铜虫 (小有名气)
送红花一朵
引用回帖:
8楼: Originally posted by 独孤神宇 at 2021-07-06 17:42:15
初始值变化很大，用 vpasolve 函数赋初始值

clear;
a=0.709133065;b=0.08555;
syms x y positive
eqns=;
=vpasolve(eqns,,)...

谢谢！（30000,0.05）这个初始值是自己根据实际情况设置，还是随意设置的呢？这次我运行出来了，但是感觉结果不是我想象中的那样。

发自小木虫Android客户端
赞一下 回复此楼
9楼2021-07-08 14:23:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 钟声入翠微 的主题更新
111/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 本子写完了，给DS兄弟看了，得了92分 +6 Doma 2026-03-01 10/500 2026-03-05 10:40 by ms629
[考研] 【求调剂】293分环境工程求调剂材料/化工，服从调剂，抗压能力强！ +6 xiiiia 2026-03-04 6/300 2026-03-05 09:44 by houyaoxu
[考研] 0856材料求调剂 +17 hyf hyf hyf 2026-02-28 18/900 2026-03-05 09:32 by oxidpl
[考研] 304分材料专硕求调剂 +11 qiuzhigril 2026-03-03 14/700 2026-03-05 09:30 by vfdbghn
[考研] 材料调剂 +3 ounce. 2026-03-04 4/200 2026-03-05 09:29 by oxidpl
[基金申请] 没有青基直接申请面上，感觉自己疯了 +6 kevin63t 2026-03-02 12/600 2026-03-04 22:23 by kevin63t
[考研] 环境调剂 +7 柒槿levana 2026-03-01 7/350 2026-03-04 22:00 by gwjxiaolang
[考研] 266求调剂 +7 哇塞王帅 2026-03-03 7/350 2026-03-04 20:30 by 一切OK
[考研] 求调剂 +4 熬夜的猫头鹰 2026-03-02 4/200 2026-03-04 17:22 by 梦天888
[考研] 本科太原理工采矿工程，求调剂 +3 onlx 2026-03-01 3/150 2026-03-04 15:57 by Stephen_ym
[考研] 能动297求调剂，本科川大 +4 邵11 2026-03-04 4/200 2026-03-04 15:29 by lylylyheart
[考研] 一志愿山东大学材料与化工325求调剂 +5 半截的诗0927 2026-03-02 5/250 2026-03-02 18:37 by 明亮9527
[考研] 284求调剂 +10 天下熯 2026-02-28 11/550 2026-03-02 11:03 by 无际的草原
[基金申请] 成果系统访问量大，请一小时后再尝试。---NSFC啥时候好哦，已经两天这样了 +4 NSFC2026我来了 2026-02-28 4/200 2026-03-01 22:37 by 铁门栓
[考博] 26申博 +4 想申博！ 2026-02-26 6/300 2026-03-01 17:32 by 想申博！
[考研] 290求调剂 +9 材料专硕调剂； 2026-02-28 11/550 2026-03-01 17:21 by sunny81
[考研] 311求调剂 +6 亭亭亭01 2026-03-01 6/300 2026-03-01 15:41 by 324616
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考研] 311求调剂 +9 南迦720 2026-02-28 10/500 2026-03-01 10:55 by sunny81
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
信息提示
关闭
请填处理意见
关闭