解方程

定义符号

syms x

等同于

sym('x')

解多元一次方程组


最后的A变量是一个struct,通过A.x查看x的值,X.y查看y的值

符号化

解方程

matlab自动认未知数,按26各英文字母排列,以距离x的远近定义,离x越近的,定为未知数。

例如解这个方程



可以用(注意在新版本matlab中,solve中的方程式不加引号)

syms x a b
solve(a*x^2-b)

结果(自动按照x为未知数)

ans =
 
  b^(1/2)/a^(1/2)
 -b^(1/2)/a^(1/2)

若要把b当作未知数,可以

syms x a b
solve(a*x^2-b,b)

求导

例如求

这个函数的导数

syms x
y = 4*x^5;
yprime=diff(y)

结果

yprime =
 
20*x^4

积分

例如求


syms x
y = x^2*exp(x);
z=int(y)

求解出

z =
 
exp(x)*(x^2 - 2*x + 2)

符号数值代入

我们如果希望把此处的x的值代入z中,可以用subs()函数

subs(z,x,0)

此处相当于求z(x=0)的值。
这样我们就能得到刚才积分的常数C了。
最终脚本

syms x
y = x^2*exp(x);
z=int(y);
z=z-subs(z,x,0)

注意,此处用int(y)求不定积分。如果要求定积分,可以用int(y,0,10)类似的语句定义上下限。

fsolve

Function Handle

matlab中,@表示Function Handle,即函数指针。
我们可以在函数调用的xy_plot(@sin,1:pi/10:pi),这里使用@相当于传入函数指针。

求解

例如求解


可以用fsolve(function_handle,init_guss_value)
求解脚本:

f2 = @(x) (1.2*x+0.3+x*sin(x));
fsolve(f2,0)

fzero

类似fsolve,都是求解为0时的解。但是fzero不能求解那些不穿过x轴,(即使和x轴相切)的方程。
例如

f =@(x) x.^2;
fzero(f,-0.1)
fsolve(f,0.1)

这里,fzero求不出解,而fsolve可以。

(单元非线性方程可用fzero,多元非线性方程可用fsolve)

optimset设置进阶选项

可以参考
optimset官方文档
optimset

例如关闭多余帮助显示

options=optimset('Display','None');
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 216,372评论 6 498
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 92,368评论 3 392
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 162,415评论 0 353
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 58,157评论 1 292
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 67,171评论 6 388
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,125评论 1 297
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,028评论 3 417
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,887评论 0 274
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,310评论 1 310
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,533评论 2 332
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,690评论 1 348
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,411评论 5 343
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,004评论 3 325
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,659评论 0 22
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,812评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,693评论 2 368
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,577评论 2 353

推荐阅读更多精彩内容

  • 目录 [TOC] 前言 sympy不仅在符号运算方面强大,在解方程方面也是很强大。本章节学习对应官网的:Solve...
    张一根阅读 7,326评论 0 3
  • 我们先从简单的来 例题1: 这是北师大版小学六年级上册课本95页的一道解方程练习题:   大家可以先口算一下,这道...
    天使猿ty阅读 1,880评论 0 1
  • 在大多数人的印象中,一谈到数学,脑海中蹦出来的词很能可能就是“无聊,枯燥,复杂”,在这种情况下,很多人对数...
    褚欢阅读 2,234评论 0 1
  • 程序的控制过程就是一个解方程的过程:用户输入有限的变量,程序根据变量求解出一些控制值,然后利用这些控制值完成指定的...
    ww4u阅读 342评论 0 0
  • 不知为什么,总有同学在解有关极坐标方程问题时,会把极坐标转化为直角坐标,用直角坐标的方法去解决问题。费时费力,得不...
    7300T阅读 3,781评论 0 3