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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 急急急！符号函数解代数方程遇到的问题
151/2返回列表12下一页
查看: 1781  |  回复: 14
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chaofan1231

银虫 (小有名气)

[求助] 急急急！符号函数解代数方程遇到的问题

各位高手帮忙看下这是怎么一回事呀！
rg=2;
>> T2=subs(fsolve('rg-exp((1/0.0000005461-1/0.000000700)*0.014388/x+log((x-993.74731)/2869.03078)-5*log(700/546.1))',1000))
??? Error using ==> inlineeval at 15
Error in inline expression ==> rg-exp((1/0.0000005461-1/0.000000700)*0.014388/x+log((x-993.74731)/2869.03078)-5*log(700/546.1))
Undefined function or variable 'rg'.
Error in ==> inline.feval at 36
        INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);
注：直接将rg=2代入方程式可以解的，但是我需要的是在循环中不断给rg赋值，然后再代入公式计算T2，不知道该如何做到。有其它能解上边代数方程的函数也可以用。谢谢啊！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-05-16 20:58:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
可以这样
CODE:
function f=do(rg)
      n=length(rg);
      x0=1000;
      f=zeros(n,1);
      for i=1:n
         f(i)=fsolve(@(x) fun(x,rg(i)),x0);
      end

function f=func(x,rg)
      f=rg-exp((1/0.0000005461-1/0.000000700)*0.014388/x+log((x-993.74731)/2869.03078)-5*log(700/546.1));
end

使用方法，将你要使用rg的值作为一个向量带入 ,运行do(rg),就ok
赞一下 回复此楼
2楼2012-05-17 12:01:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)

【答案】应助回帖

引用回帖:
2楼: Originally posted by lijie169 at 2012-05-17 12:01:53:
可以这样

function f=do(rg)
      n=length(rg);
      x0=1000;
      f=zeros(n,1);
      for i=1:n
         f(i)=fsolve(@(x) fun(x,rg(i)),x0);
      end

function f=func(x,rg)
      f=rg ...

CODE:
function f=do(rg)
      n=length(rg);
      x0=1000;
      f=zeros(n,1);
      for i=1:n
         f(i)=fsolve(@(x) func(x,rg(i)),x0);
      end

function f=func(x,rg)
      f=rg-exp((1/0.0000005461-1/0.000000700)*0.014388/x+log((x-993.74731)/2869.03078)-5*log(700/546.1));
end

有点小错误....fun应该是func
赞一下 回复此楼
3楼2012-05-17 12:03:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
chaofan1231: 金币+5, ★★★★★最佳答案, 谢谢 2012-05-17 13:43:29
xzhdty: 金币+2, 谢谢参与 2012-05-19 16:43:29
可以吧for 给去掉，只留一句  f=fsolve(@(x) func(x,rg),1000);后面的function再接上就可以了...不需要那么麻烦
CODE:
function f=do(rg)

         x0=1000;
         f(i)=fsolve(@(x) func(x,rg),x0);

function f=func(x,rg)
      f=rg-exp((1/0.0000005461-1/0.000000700)*0.014388/x+log((x-993.74731)/2869.03078)-5*log(700/546.1));
end

简洁就是美
赞一下 回复此楼

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

chaofan1231
4楼2012-05-17 12:08:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaofan1231

银虫 (小有名气)
送鲜花一朵
引用回帖:
4楼: Originally posted by lijie169 at 2012-05-17 12:08:00:
可以吧for 给去掉，只留一句  f=fsolve(@(x) func(x,rg),1000);后面的function再接上就可以了...不需要那么麻烦

function f=do(rg)

         x0=1000;
         f(i)=fsolve(@(x) func(x,rg),x0);

fun ...

再问一下，这是我循环中的一部分，用这个函数求解时每次循环都要迭代，运算量大，很 慢，有什么更好的办法吗？
赞一下 回复此楼
5楼2012-05-17 13:45:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)
引用回帖:
5楼: Originally posted by chaofan1231 at 2012-05-17 13:45:02:
再问一下，这是我循环中的一部分，用这个函数求解时每次循环都要迭代，运算量大，很 慢，有什么更好的办法吗？

你可以采用我最下面的代码，支持向量运算。一次性算完，采用矩阵的运算方法，听说matlab在矩阵运算里比for循环快N倍...你可以试试。
赞一下 回复此楼
6楼2012-05-17 16:26:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaofan1231

银虫 (小有名气)
引用回帖:
6楼: Originally posted by lijie169 at 2012-05-17 16:26:17:
你可以采用我最下面的代码，支持向量运算。一次性算完，采用矩阵的运算方法，听说matlab在矩阵运算里比for循环快N倍...你可以试试。

不能算啊！Tsolve函数是我创建的M文件。跟你说的Func是一样的。只能算出来一个值。
rg =

    2.2000    2.0000    1.0000    1.2000    1.1000    1.5000    0.8000
    1.3000    1.5000    2.5000    2.6000    0.9000    1.1000    2.2000
    1.1000    1.3000    2.2000    0.9000    2.2000    2.2000    2.2000

>> t=fsolve(@(x)Tsolve(x,rg),1000)
Warning: Trust-region-dogleg algorithm of FSOLVE cannot handle non-square systems; using Levenberg-Marquardt algorithm instead.
> In fsolve at 324

No solution found.

fsolve stopped because the last step was ineffective. However, the vector of function
values is not near zero, as measured by the default value of the function tolerance.


t =

  1.0628e+003
赞一下 回复此楼
7楼2012-05-18 08:44:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chaofan1231

银虫 (小有名气)
引用回帖:
6楼: Originally posted by lijie169 at 2012-05-17 16:26:17:
你可以采用我最下面的代码，支持向量运算。一次性算完，采用矩阵的运算方法，听说matlab在矩阵运算里比for循环快N倍...你可以试试。

我想用矩阵运算的，但是我的方程太复杂了，含有指数，不知道怎么化成矩阵运算的形式。
赞一下 回复此楼
8楼2012-05-18 08:47:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)
貌似是用3楼的代码可以，我理解错了，不能使用直接向量求解。使用规则是rg是向量，不是矩阵，若要是矩阵，则需要修改主函数代码
CODE:
function f=do(rg)
      [m,n]=size(rg);
      x0=1000;
      f=zeros(m,n);
      for i=1:m
      for j=1:n
         f(i,j)=fsolve(@(x) func(x,rg(i,j)),x0);
      end
     end
function f=func(x,rg)
      f=rg-exp((1/0.0000005461-1/0.000000700)*0.014388/x+log((x-993.74731)/2869.03078)-5*log(700/546.1));
end
end
赞一下 回复此楼
9楼2012-05-18 12:00:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijie169

铜虫 (著名写手)
引用回帖:
9楼: Originally posted by lijie169 at 2012-05-18 12:00:13:
貌似是用3楼的代码可以，我理解错了，不能使用直接向量求解。使用规则是rg是向量，不是矩阵，若要是矩阵，则需要修改主函数代码

function f=do(rg)
      =size(rg);
      x0=1000;
      f=zeros(m,n);
...

最后的代码rg可以是矩阵也可以是向量，3楼的rg不能是矩阵
赞一下 回复此楼
10楼2012-05-18 12:01:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chaofan1231 的主题更新
151/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国自然评审意见 +7 wangmingqi 2026-08-28 8/400 2026-08-28 18:42 by zhuzg0628
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 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
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭