matlab中符号常量用sym函数如何定义？ MATLAB中定义符号变量出错，sym函数怎么用啊？
matlab\u4e2d\u7b26\u53f7\u5e38\u91cf\u7528sym\u51fd\u6570\u5982\u4f55\u5b9a\u4e49\u5728MATLAB\u4e2d\uff0csym\u51fd\u6570\u7528\u4e8e\u521b\u5efa\u7b26\u53f7\u6570\u5b57\u3001\u7b26\u53f7\u53d8\u91cf\u3001\u7b26\u53f7\u5bf9\u8c61\u3002\u7b26\u53f7\u5bf9\u8c61\u7684\u7c7b\u578b\u662fsym\uff0c\u53ef\u4ee5\u901a\u8fc7class(S)\u6765\u9a8c\u8bc1\uff0c\u8fd9\u91ccS\u662f\u4e00\u4e2a\u7b26\u53f7\u5bf9\u8c61\u3002
\u7b26\u53f7\u53d8\u91cf\u7684\u4f18\u70b9\u662f\uff0c\u4f7f\u7528\u7b26\u53f7\u53d8\u91cf\u8fd0\u7b97\u5f97\u5230\u7684\u53ea\u662f\u4e00\u4e2a\u89e3\u6790\u89e3\uff0c\u4f8b\u5982\uff0c\u5728\u7b26\u53f7\u53d8\u91cf\u8fd0\u7b97\u8fc7\u7a0b\u4e2dpi\u5c31\u7528pi\u8868\u793a\uff0c\u800c\u4e0d\u662f\u5177\u4f53\u7684\u8fd1\u4f3c\u6570\u503c3.14\u6216
3.14159\u3002\u4f7f\u7528\u7b26\u53f7\u53d8\u91cf\u8fdb\u884c\u8fd0\u7b97\u80fd\u6700\u5927\u9650\u5ea6\u51cf\u5c11\u8fd0\u7b97\u8fc7\u7a0b\u4e2d\u56e0\u820d\u5165\u9020\u6210\u7684\u8bef\u5dee\u3002\u7b26\u53f7\u53d8\u91cf\u4e5f\u4fbf\u4e8e\u8fdb\u884c\u8fd0\u7b97\u8fc7\u7a0b\u7684\u6f14\u793a\u3002
\u8bed\u6cd5\u683c\u5f0f\uff1a
S = sym(A)\u5c06\u975e\u7b26\u53f7\u5bf9\u8c61\uff08\u5982\uff0c\u6570\u5b57\uff0c\u8868\u8fbe\u5f0f\uff0c\u53d8\u91cf\u7b49\uff09A\u8f6c\u6362\u4e3a\u7b26\u53f7\u5bf9\u8c61\uff0c\u5e76\u5b58\u50a8\u5728\u7b26\u53f7\u53d8\u91cfS\u4e2d\u3002
x = sym('x')
\u521b\u5efa\u7b26\u53f7\u53d8\u91cfx\uff0c\u5176\u540d\u5b57\u662f'x'\u3002\u793a\u4f8b\uff1aalpha = sym('alpha')
x = sym('x', 'real')
\u8fd9\u91cc\u5047\u8bbex\u662f\u5b9e\u6570\uff0c\u56e0\u6b64\u6709x\u7684\u5171\u8f6dconj(x)\u7b49\u4e8ex\u3002\u793a\u4f8b\uff1ar = sym('Rho','real')
k = sym('k', 'positive')
\u8fd9\u91cc\u521b\u5efa\u4e00\u4e2a\u6b63\u7684\uff08\u5b9e\u6570\uff09\u7b26\u53f7\u53d8\u91cf\u3002
x = sym('x', 'clear')
\u521b\u5efa\u4e00\u4e2a\u6ca1\u6709\u989d\u5916\u5c5e\u6027\u7684\u7eaf\u5f62\u5f0f\u4e0a\u7684\u7b26\u53f7\u53d8\u91cfx\uff08\u4f8b\u5982\uff0c\u521b\u5efa\u7b26\u53f7\u53d8\u91cfx\uff0c\u4f46\u662f\u5e76\u6ca1\u6307\u5b9a\u5b83\u662f\u6b63\u7684\u6216\u5b83\u662f\u4e00\u4e2a\u5b9e\u6570\uff09\u3002\u4e3a\u4e86\u517c\u5bb9\u65e7\u7684MATLAB\u7248\u672c\uff0cx = sym('x','unreal')\u7684\u529f\u80fd\u548cx = sym('x', 'clear')\u4e00\u6837\u3002
S = sym(A, flag)
\u628a\u4e00\u4e2a\u6570\u503c\u6807\u91cf\u6216\u77e9\u9635\u8f6c\u6362\u4e3a\u7b26\u53f7\u578b\u7684\u5bf9\u8c61\u3002\u8fd9\u91ccflag\u53c2\u6570\u7684\u503c\u53ef\u4ee5\u662f\uff1a'r', 'd', 'e', or 'f'\uff0c\u5b83\u6307\u5b9a\u4e86\u5bf9\u6d6e\u70b9\u6570\u8fdb\u884c\u8f6c\u6362\u65f6\u7684\u89c4\u5219\uff1a
'f':\u8868\u793a\u201cfloating-point\u201d\u3002\u8fd9\u6837\uff0c\u6240\u6709\u7684\u6570\u503c\u90fd\u7528N*2^e\u6216-N*2^e\u8fd9\u79cd\u5f62\u5f0f\u8868\u793a\u3002N\u548ce\u90fd\u4e3a\u6574\u6570\uff0c\u4e14N\u4e0d\u5c0f\u4e8e0\u3002\u4f8b\u5982\uff0csym(1/10, 'f')\u88ab\u8f6c\u6362\u4e3a3602879701896397/36028797018963968\u3002
'r':
\u8868\u793a\u201crational\u201d\uff08\u6709\u7406\u6570\u5f62\u5f0f\uff09\u3002\u8fd9\u91cc\uff0c\u6240\u6709\u7684\u6d6e\u70b9\u6570\u90fd\u88ab\u8868\u793a\u6210\u5f62\u5982p/q\uff08\u5206\u6570\u5f62\u5f0f\uff09\u3001p*pi/q\u30012^q\u300110^q\u3001sqrt(p)\u7b49\u5f62\u5f0f
\uff08p\u3001q\u90fd\u662f\u7b26\u53f7\u578b\u7684\u6709\u7406\u6570\uff09\u3002\u8fd9\u79cd\u8868\u793a\u65b9\u6cd5\u51cf\u5c11\u4e86\u8868\u8fbe\u5f0f\u4e2d\u56e0\u820d\u5165\u9020\u6210\u7684\u8bef\u5dee\u3002\u4f46\u6709\u65f6\u5019\u8fd9\u79cd\u65b9\u6cd5\u5e76\u4e0d\u80fd\u51c6\u786e\u8868\u793a\u4e00\u4e2a\u503c\u3002\u5982\u679c\u627e\u4e0d\u5230\u6700\u903c\u8fd1\u7684\u6709\u7406\u6570\u6765\u8868\u793a\u4e00
\u4e2a\u6d6e\u70b9\u6570\uff0c\u8fd9\u4e2a\u6570\u5c31\u4f1a\u88ab\u8868\u793a\u6210p*2^q\uff08\u5176\u4e2dp\u3001q\u90fd\u662f\u8f83\u5927\u7684\u6574\u6570\uff09\u8fd9\u79cd\u5f62\u5f0f\uff0c\u4ee5\u671f\u6765\u51c6\u786e\u7684\u8868\u793a\u5b83\u7684\u503c\u3002\u4f8b\u5982\uff0csym(4/3,'r')\u5c06\u4ea7\u751f\u7b26\u53f7\u91cf
'4/3'\uff0c\u4f46sym(1+sqrt(5),'r')\u5c06\u88ab\u8868\u793a\u62107286977268806824*2^(-51)\u3002
'e':\u8868\u793a
\u201cestimate
error\u201d\u3002\u8fd9\u79cd\u5f62\u5f0f\u4f1a\u5728'r'\u7684\u57fa\u7840\u4e0a\u6dfb\u52a0\u4e00\u4e2a\u7531\u5305\u542beps\u7684\u7b26\u53f7\u8868\u8fbe\u5f0f\u8868\u793a\u7684\u8bef\u5dee\u503c\u3002\u8fd9\u4e2a\u8bef\u5dee\u503c\u8868\u793a\u4e00\u4e2a\u6709\u7406\u6570\u7684\u7406\u8bba\u503c\u548c\u7528\u6d6e\u70b9\u6570\u5f62\u5f0f\u8868\u793a\u7684\u5b9e\u9645\u503c
\u4e4b\u95f4\u7684\u8bef\u5dee\u3002\u4f8b\u5982\uff0csym(3*pi/4,'e')\u5c06\u88ab\u8868\u793a\u62103*pi/4*(1+3143276*eps/65)\u3002
'd':\u8868\u793a
\u201cdecimal\u201d\uff08\u5341\u8fdb\u5236\uff09\u3002\u6211\u4eec\u77e5\u9053\uff0c\u5728\u5b9e\u9645\u751f\u6d3b\u4e2d\u6211\u4eec\u4e60\u60ef\u7528\u5341\u8fdb\u5236\uff0c\u4f46\u8ba1\u7b97\u673a\u4e2d\u5219\u7528\u4e8c\u8fdb\u5236\uff0c\u4e00\u4e2a\u7b80\u5355\u7684\u5341\u8fdb\u5236\u6d6e\u70b9\u65703.14\u7528\u4e8c\u8fdb\u5236\u8868\u793a\u4fbf\u4e0d\u80fd\u51c6\u786e
\u7684\u8868\u793a\u3002\u56e0\u6b64\uff0c\u66f4\u597d\u7684\u51cf\u5c11\u8bef\u5dee\u7684\u65b9\u6cd5\u4fbf\u662f\u5728\u8ba1\u7b97\u673a\u4e2d\u5f15\u5165\u5341\u8fdb\u5236\u6765\u8868\u793a\u4e00\u4e2a\u6570\u3002\u8fd9\u79cd\u60c5\u51b5\u4e0b\uff0c\u6570\u5b57\u90fd\u53d6\u81ea\u8c03\u7528vpa\u51fd\u6570\u540e\u4ea7\u751f\u7684\u6570\u5b57\u96c6\u5408\u3002\u867d\u7136\u752816\u4f4d\u6709\u65f6\u5019
\u5e76\u4e0d\u80fd\u51c6\u786e\u8868\u793a\u4e00\u4e2a\u6d6e\u70b9\u6570\uff0c\u4f46\u5927\u591a\u6570\u60c5\u51b5\u4e0b\uff0c\u6211\u4eec\u4e5f\u8bb8\u5e76\u4e0d\u5e0c\u671b\u7528\u8d85\u8fc716\u4f4d\u6570\u5b57\u6765\u8868\u793a\u4e00\u4e2a\u6d6e\u70b9\u6570\u3002\u4f8b\u5982\uff0c\u901a\u8fc7\u8c03\u7528digits(10)\uff0c\u6211\u4eec\u5f97\u5230\u7684\u6d6e\u70b9\u6570
\u90fd\u753110\u4f4d\u6570\u5b57\u6784\u6210\uff0csym(4/3,'d')\u5c06\u4ea7\u751f1.333333333\u8fd9\u4e2a\u7b26\u53f7\u6570\u5b57\u3002\u867d\u7136\u8c03\u7528digits(20)\u540e\u6211\u4eec\u53ef\u4ee5\u4f7f\u752820\u4f4d\u6570\u5b57\u6765\u8868\u793a
\u4e00\u4e2a\u6d6e\u70b9\u6570, \u8fd9\u6837sym(4/3,'d')\u5c31\u662f
1.3333333333333332593\uff0c\u4f46\u662f\u6211\u4eec\u770b\u5230\u53ea\u6709\u524d16\u4f4d\u6570\u503c\u662f\u51c6\u786e\u7684\uff0c\u540e\u9762\u7684\u51e0\u4f4d\u6570\u503c\u5df2\u7ecf\u4ea7\u751f\u4e86\u8bef\u5dee\uff0c\u56e0\u6b64\u5927\u591a\u6570\u60c5\u51b5\u4e0b\u7528\u8d85\u8fc716\u4f4d\u7684\u6570
\u5b57\u6765\u8868\u793a\u4e00\u4e2a\u6d6e\u70b9\u6570\u662f\u4e0d\u5fc5\u8981\u7684\u3002
\u56e0\u4e3a\u6ca1\u6709\u770b\u5230\u62a5\u9519\u4fe1\u606f\uff0c\u6240\u4ee5\u53ea\u80fd\u731c\u4e86\u3002
\u5728for\u5faa\u73af\u4e0a\u9762\u4e00\u884c\uff0c\u56e0\u4e3a\u4f60\u7528aw=zeros(1, n-m); \u5c06aw\u5b9a\u4e49\u4e3a\u4e86\u53cc\u7cbe\u5ea6\uff0c\u6240\u4ee5\u5b83\u5b58\u4e0d\u4e86\u7b26\u53f7\u3002
aw = sym(zeros(1,n-m)); % \u8fd9\u6837\u5c31\u53ef\u4ee5\u4e86
符号变量的优点是，使用符号变量运算得到的只是一个解析解，例如，在符号变量运算过程中pi就用pi表示，而不是具体的近似数值3.14或
3.14159。使用符号变量进行运算能最大限度减少运算过程中因舍入造成的误差。符号变量也便于进行运算过程的演示。
语法格式：
S = sym(A)将非符号对象（如，数字，表达式，变量等）A转换为符号对象，并存储在符号变量S中。
x = sym('x')
创建符号变量x，其名字是'x'。示例：alpha = sym('alpha')
x = sym('x', 'real')
这里假设x是实数，因此有x的共轭conj(x)等于x。示例：r = sym('Rho','real')
k = sym('k', 'positive')
这里创建一个正的（实数）符号变量。
x = sym('x', 'clear')
创建一个没有额外属性的纯形式上的符号变量x（例如，创建符号变量x，但是并没指定它是正的或它是一个实数）。为了兼容旧的MATLAB版本，x = sym('x','unreal')的功能和x = sym('x', 'clear')一样。
S = sym(A, flag)
把一个数值标量或矩阵转换为符号型的对象。这里flag参数的值可以是：'r', 'd', 'e', or 'f'，它指定了对浮点数进行转换时的规则：
'f':表示“floating-point”。这样，所有的数值都用N*2^e或-N*2^e这种形式表示。N和e都为整数，且N不小于0。例如，sym(1/10, 'f')被转换为3602879701896397/36028797018963968。
'r':
表示“rational”（有理数形式）。这里，所有的浮点数都被表示成形如p/q（分数形式）、p*pi/q、2^q、10^q、sqrt(p)等形式
（p、q都是符号型的有理数）。这种表示方法减少了表达式中因舍入造成的误差。但有时候这种方法并不能准确表示一个值。如果找不到最逼近的有理数来表示一
个浮点数，这个数就会被表示成p*2^q（其中p、q都是较大的整数）这种形式，以期来准确的表示它的值。例如，sym(4/3,'r')将产生符号量
'4/3'，但sym(1+sqrt(5),'r')将被表示成7286977268806824*2^(-51)。
'e':表示
“estimate
error”。这种形式会在'r'的基础上添加一个由包含eps的符号表达式表示的误差值。这个误差值表示一个有理数的理论值和用浮点数形式表示的实际值
之间的误差。例如，sym(3*pi/4,'e')将被表示成3*pi/4*(1+3143276*eps/65)。
'd':表示
“decimal”（十进制）。我们知道，在实际生活中我们习惯用十进制，但计算机中则用二进制，一个简单的十进制浮点数3.14用二进制表示便不能准确
的表示。因此，更好的减少误差的方法便是在计算机中引入十进制来表示一个数。这种情况下，数字都取自调用vpa函数后产生的数字集合。虽然用16位有时候
并不能准确表示一个浮点数，但大多数情况下，我们也许并不希望用超过16位数字来表示一个浮点数。例如，通过调用digits(10)，我们得到的浮点数
都由10位数字构成，sym(4/3,'d')将产生1.333333333这个符号数字。虽然调用digits(20)后我们可以使用20位数字来表示
一个浮点数, 这样sym(4/3,'d')就是
1.3333333333333332593，但是我们看到只有前16位数值是准确的，后面的几位数值已经产生了误差，因此大多数情况下用超过16位的数
字来表示一个浮点数是不必要的。
绛旓細鍦∕ATLAB涓紝sym鍑芥暟鐢ㄤ簬鍒涘缓绗﹀彿鏁板瓧銆佺鍙峰彉閲忋佺鍙峰璞銆傜鍙峰璞＄殑绫诲瀷鏄痵ym锛屽彲浠ラ氳繃class(S)鏉ラ獙璇侊紝杩欓噷S鏄竴涓鍙峰璞°俓x0d\x0a绗﹀彿鍙橀噺鐨勪紭鐐规槸锛屼娇鐢ㄧ鍙峰彉閲忚繍绠楀緱鍒扮殑鍙槸涓涓В鏋愯В锛屼緥濡傦紝鍦ㄧ鍙峰彉閲忚繍绠楄繃绋嬩腑pi灏辩敤pi琛ㄧず锛岃屼笉鏄叿浣撶殑杩戜技鏁板3.14鎴朶x0d\x0a3.14159銆備娇鐢...
绛旓細syms鏄畾涔夌鍙峰彉閲 sym鏄皢瀛楃鎴栬呮暟瀛楄浆鎹负瀛楃 姣斿 syms x y %灏辨槸瀹氫簡绗﹀彿鍙橀噺x y浠ュ悗x y灏卞彲浠ョ洿鎺ヤ娇鐢ㄤ簡锛屾湁浠栦滑杩愮畻鍑烘潵鐨勭粨鏋滀篃鏄鍙峰彉閲 褰撶劧涓婇潰鐨勪篃鍙互x=sym('x'),y=sym('y')sys('a+b')%灏辨槸灏哸+b杞寲涓虹鍙疯〃杈惧紡 涓よ呮湁涓嶅悓鐨勪娇鐢ㄥ満鍚 姣斿绗﹀彿鍙橀噺娑夊強鐨勪釜鏁版瘮杈冨皯锛...
绛旓細浣跨敤鏁板艰〃杈惧紡鐨勫彉閲忓繀椤婚鍏堣鍙橀噺璧嬪硷紝鍜鍦∕atlab涓瀹炵敤鏁板艰绠椾竴鏍凤紝鍦ㄨ繘琛岀鍙疯繍绠椾箣鍓嶏紝棣栧厛闇瑕佸畾涔夌鍙峰彉閲忋傚笇鏈涙垜鐨勫洖绛斿浣犳湁鎵甯姪銆傞棶棰樹簩锛氫粈涔堟槸绗﹀彿甯搁噺锛浣跨敤绗﹀彿甯搁噺鏈変粈涔堝ソ澶勶紵 绗﹀彿甯搁噺灏辨槸鐢╠efine瀹氫箟鐨勫父閲忥紝姣斿define a =1;浣犱互鍚庡氨鍙互鎶奱褰1鐢紝绋嬪簭閲屽彧瑕佹湁1鐨勫湴鏂逛綘閮藉啓...
绛旓細鍙互鐢 matlab 鐨勭鍙璁＄畻鍖(Symbolic Math Toolbox)锛屾瘮濡 syms x a b;f = a*x^2 + b;dfdx = diff(f, x)杈撳嚭 dfdx = 2*a*x鍏鐢 syms 澹版槑鎵鏈夌殑鍙橀噺銆俤iff 鍙互鎸囧畾姹傚鐨勫彉閲忥紝鍏跺畠灏卞綋鍙傛暟鐪嬩簡銆傚鏋滀綘纭疄鐭ラ亾鏌愪簺鍙橀噺鐨勫硷紝閭ｄ箞涔熷彲浠ヨ繖鏍凤細syms x a b;a = 2;f = a*x^2...
绛旓細绗﹀彿甯搁噺鐭╅樀A锛屽彲浠浣跨敤double鏉ヨ浆鎹负鏁板肩煩闃点>> A = sym([1/2 1/3;1/4 1/5])A = [ 1/2, 1/3][ 1/4, 1/5]>> A1=double(A)A1 = 0.5000 0.3333 0.2500 0.2000 鑻ユ槸绗﹀彿鍙橀噺鐭╅樀锛屾瘮濡傚惈鏈夊彉閲弜锛岄渶瑕佷娇鐢ㄥ彉閲忔浛鎹ubs锛屾墠鑳藉緱鍒版暟鍊肩煩闃 >> syms x >> A = [2*x ...
绛旓細绗﹀彿鍙橀噺瑕鐢⊿YMS鏉ュ０鏄庯紝MATLAB涓灏嗗叾璇嗗埆涓轰竴涓甯搁噺锛屽氨鍍忎竴涓粰浜嗗肩殑甯告暟涓鏍枫傚洜涓篗ATLAB鏄笉鑳借瘑鍒鍙风殑锛屽０鏄庝互鍚庣郴缁熷氨鎶婄鍙峰彉閲忓綋浣滀竴涓父閲忓寰呫備緥濡傚湪S鍙樻崲涓紝灏辫澹版槑S涓虹鍙峰彉閲忥紝涓嶇劧锛岀郴缁熶笉鑳借瘑鍒玈銆傜煡閬撲簡绗﹀彿鍙橀噺鐨勪綔鐢紝鐭╅樀鐨勫氨涓嶇敤鎴戝璇翠簡鍚с傚湪涓涓弬鏁板彇澶氫釜鍊兼椂锛屽彲瀹氫箟...
绛旓細1銆佹ゼ涓绘墍绀鐨鏂规硶鏄鍏堢粰A銆丅閮借祴鍊间簡锛岀劧鍚嶤鐨勮绠楃粨鏋滃氨浼氭槸涓涓叿浣撶殑鍊硷紝鑰屼笉鏄〃杈惧紡銆侰鏄彲浠ョ畻鍑烘潵鐨勶紝涓嶈繃缁撴灉鏄疌=5銆2銆佽嫢鏄兂鐩存帴寰楀埌绗﹀彿琛ㄨ揪寮忥紝鍙互鐩存帴璁剧疆绗﹀彿鍙橀噺锛屽涓嬫墍绀猴細>> syms A B C=[A,B]*[A;B]C = A^2+B^2 3銆佹槸鍙互鐢╯ubs鏇挎崲鍑芥暟瀵硅〃杈惧紡杩涜鏇挎崲璁＄畻鐨...
绛旓細鎯呭喌涓锛绗﹀彿鍙橀噺 蹇呴』瑕佸畾涔夛紝瀹氫箟涔嬪悗鎵嶅彲浠浣跨敤銆傚涓嬮潰鐨渚嬪瓙锛屾湭瀹氫箟灏变細鍑洪敊銆傜鍙峰彉閲忓畾涔夋柟娉曚緥濡: sym x 鎴栬 syms x 锛屼袱鑰呮湁鍖哄埆涔熸湁鍏卞悓鐐广傛儏鍐典簩锛氳绠楄繃绋嬪綋涓殑瀛樺偍鏁板肩殑鍙橀噺 鍙互涓嶅繀瀹氫箟锛岄殢鏃堕渶瑕侊紝闅忔椂瀹氫箟锛屼絾鏄湁鏃跺欏鏋滃彉閲忓緢澶氾紝閭ｄ箞鏈濂芥彁鍓嶅０鏄庯紝澹版槑鐨勬椂鍊欙紝鍙互鐩存帴璧嬩簣0鍊...
绛旓細A=magic(5)sym(A)
绛旓細1銆佷笉闇瑕浣跨敤syms瀹氫箟閭ｅ嚑涓绗﹀彿閲忋2銆佷娇鐢╯ubs鍑芥暟鎶婃暟鎹唬鍏ヨ〃杈惧紡銆3銆佷娇鐢╯olve鍑芥暟姹傝В锛屽苟浣跨敤double鍑芥暟杞崲涓烘诞鐐规暟銆4銆佺粰鍑虹粯鍥剧粨鏋滐紝鍙互鐩磋鍦伴獙璇佽В鐨姝ｇ‘鎬с傚弬鑰冧唬鐮 M=2e+5;G=2e+6;k=2;y=dsolve('Dy=M*vh/(G-k*vh)','y(0)=0','vh');% 寰楀埌y鍏充簬vh鐨勮〃杈惧紡y=subs...
