24小时热门版块排行榜    

查看: 1203  |  回复: 4

jiweibang

木虫 (正式写手)


[交流] 【求助】发贴求助,关于画一个隐函数的程序问题

我现在要解一个方程,像这样的,其中K1和K2是关于x,y的函数K1=(2x-y)^1/2。我考虑是个隐函数求根的问题,可是我Fortran学的不到家,不知道能用什么方法解决,请各位高手帮帮忙,提个思路就行,先谢谢了。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

scatcj

木虫 (小有名气)



jiweibang(金币+1): 2011-02-23 08:49:49
余泽成(金币+1): 谢谢参与应助! 2011-02-27 15:04:35
先把式子搞清楚了,解方程的话,Fortran书上基本都有很经典的方法的
2楼2011-02-23 00:02:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (知名作家)


★ ★
jiweibang(金币+5): 谢了,搜过,没有启发性的东西 2011-02-23 22:17:55
余泽成(金币+2): 谢谢参与应助! 2011-02-27 15:04:08
本帖仅楼主可见
3楼2011-02-23 17:05:41
已阅   申请程序强帖   回复此楼   编辑   查看我的主页

锐利的碎片

木虫 (正式写手)


★ ★
jiweibang(金币+10): 感谢关注,我就是光用二分法弄不出来 2011-02-23 22:17:14
余泽成(金币+2): 谢谢参与应助! 2011-02-27 15:04:16
可以试下多变量的优化算法,比如模拟退火来搜索局部或者全局最小点的位置,找到后可以在附近用二分法搜索。

[ Last edited by 锐利的碎片 on 2011-2-23 at 19:57 ]
4楼2011-02-23 19:53:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

leedobb

金虫 (正式写手)


★ ★ ★
jiweibang(金币+10): 先谢谢,貌似没这么简单,不过我目前也不是很清楚 2011-02-24 16:11:59
余泽成(金币+3): 谢谢参与应助! 2011-02-27 15:04:26
请问您想要求解哪个量,是x和y吗?那这样很简单啊,
以f(x,y)=2为例
可以有两种方案,方案一:
令 F(x,y) = (f(x,y)-2)^2
给一初值(x0,y0)
用下面的方程
dx/dt = -k * d F/d x
dy/dt = -k* dF/dy
演化x,y直至f(x,y)接近2即可,尝试不同的初值即可。

方案二:
直接编个二重循环
do x=-100, 100, 0.01
do y=-100, 100, 0.01
计算delt =f(x,y)-2
enddo
enddo
把最接近0的那个取出就行了。


我个人认为此题用方案2更恰当。因为变量少,没必要折腾。
5楼2011-02-24 12:20:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jiweibang 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 5/250 2026-09-01 05:45 by o3EMxHUAQDtS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 EZ67DG6CvcKL 2026-08-31 4/200 2026-09-01 01:30 by o3EMxHUAQDtS
[基金申请] 国社科又开始会评了,不知道这次命运如何 +7 雨打竹帘 2026-08-30 11/550 2026-08-31 23:16 by hittle2008
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 EZ67DG6CvcKL 2026-08-31 3/150 2026-08-31 22:23 by o3EMxHUAQDtS
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 基金不中,共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[基金申请] 为什么到现在没收到通知? +5 tannykie 2026-08-29 5/250 2026-08-30 21:05 by purplejack
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 2026年叶企孙基金 +4 bud_bud 2026-08-27 7/350 2026-08-29 07:23 by foolishmani
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 看板上这么多中的,有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
信息提示
请填处理意见