| 查看: 5438 | 回復(fù): 12 | ||
nuomandi521金蟲 (小有名氣)
|
[求助]
Matlab求解a和b分別是不同值時 f(a,b,x)=0 的解 已有1人參與
|
|
|
希望求解 f(a,b,x)=0 的方程,求解a和b分別是不同值時候,x的值(大概率是復(fù)數(shù))。 Matlab 運行如下程序后, clear all syms a b x; a=0.9; b=0.3; x=solve('4*x^2*(x^2-1)*a^2+8*x^2*((sin(x*pi/2))^2-x^2)*a*b+4*((sin(x*pi/2))^2-x^2)^2*b^2+(sin(x*pi))^2','x'); x=subs(x); 出現(xiàn)問題, 警告: Support of character vectors that are not valid variable names or define a number will be removed in a future release. To create symbolic expressions, first create symbolic variables and then use operations on them. > In sym>convertExpression (line 1559) In sym>convertChar (line 1464) In sym>tomupad (line 1216) In sym (line 179) In solve>getEqns (line 405) In solve (line 225) 警告: Do not specify equations and variables as character vectors. Instead, create symbolic variables with syms. > In solve>getEqns (line 445) In solve (line 225) 警告: Cannot find explicit solution. > In solve (line 316) 多謝指教! |
版主 (知名作家)
|
真是無語,你直接將問題描述清楚不就行了。。 clear; a=-0.6; b=-0.9; x=-3:0.01:3; f=4*x.^2.*(x.^2-1)*a^2+8*x.^2.*((sin(x*pi/2)).^2-x.^2)*a*b+4*((sin(x*pi/2)).^2-x.^2).^2*b^2+(sin(x*pi)).^2; z=[]; j=0; for i=1:length(x) if abs(f(i))==0 j=j+1; z(j)=x(i); end end z %%%%************************* clear; a=-0.6; b=-0.9; x=-3:0.01:3; f=4*x.^2.*(x.^2-1)*a^2+8*x.^2.*((sin(x*pi/2)).^2-x.^2)*a*b+4*((sin(x*pi/2)).^2-x.^2).^2*b^2+(sin(x*pi)).^2; z=[]; j=0; for i=1:length(x) if abs(f(i))<0.0001 j=j+1; z(j)=x(i); end end z |

版主 (知名作家)
|
clear; f=@(a,b)@(x)4*x^2*(x^2-1)*a^2+8*x^2*((sin(x*pi/2))^2-x^2)*a*b+4*((sin(x*pi/2))^2-x^2)^2*b^2+(sin(x*pi))^2; x0=fzero(f(0.3,0.9),2) |

金蟲 (小有名氣)
版主 (知名作家)

金蟲 (小有名氣)
金蟲 (小有名氣)
版主 (知名作家)
|
matlab 基礎(chǔ)要好好學(xué)一下了。。 a=-0.6; b=-0.9; x = linspace(-3,3); y=4*x.^2.*(x.^2-1)*a^2+8*x.^2.*((sin(x*pi/2)).^2-x.^2)*a*b+4*((sin(x*pi/2)).^2-x.^2).^2*b^2+(sin(x*pi)).^2; plot(x,y); |

金蟲 (小有名氣)
金蟲 (小有名氣)
|
我希望解的精度再高一些,于是編寫程序如下, f=@(a,b)@(x)4*x.^2.*(x.^2-1)*a^2+8*x.^2.*((sin(x*pi/2)).^2-x.^2)*a*b+4*((sin(x*pi/2)).^2-x.^2).^2*b^2+(sin(x*pi)).^2; x=-3:0.01:3;z=[]; for i=1:length(x) if abs(f(x(i)))<0.01 z=[z x(i)]; end end z 運行后出現(xiàn)問題, 未定義與 'function_handle' 類型的輸入?yún)?shù)相對應(yīng)的函數(shù) 'abs'。 請問如何解決? |
版主 (知名作家)
|
clear; a=-0.6; b=-0.9; x=-3:0.01:3; f=4*x.^2.*(x.^2-1)*a^2+8*x.^2.*((sin(x*pi/2)).^2-x.^2)*a*b+4*((sin(x*pi/2)).^2-x.^2).^2*b^2+(sin(x*pi)).^2; z=[]; for i=1:length(x) if abs(f(i))<0.01 z(i)=f(i); end end z; %%%%********************** clear; a=-0.6; b=-0.9; x=-3:0.01:3; f=4*x.^2.*(x.^2-1)*a^2+8*x.^2.*((sin(x*pi/2)).^2-x.^2)*a*b+4*((sin(x*pi/2)).^2-x.^2).^2*b^2+(sin(x*pi)).^2; z=[]; j=0; for i=1:length(x) if abs(f(i))<0.01 j=j+1; z(j)=f(j); end end z; |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 266分,一志愿電氣工程,本科材料,求材料專業(yè)調(diào)劑 +7 | 哇呼哼呼哼 2026-04-02 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 085801 總分275 本科新能源 求調(diào)劑 +18 | bradoner 2026-04-01 | 22/1100 |
|
|
[考研] 316求調(diào)劑 +12 | 舟自梗 2026-04-01 | 16/800 |
|
|
[考研] 262求調(diào)劑 +5 | 勵志一定發(fā)文章 2026-04-02 | 6/300 |
|
|
[考研] 一志愿北京科技大學(xué)085601材料工程英一數(shù)二初試總分335求調(diào)劑 +9 | 雙馬尾痞老板2 2026-04-01 | 9/450 |
|
|
[考研] 求調(diào)劑!生物與醫(yī)藥專碩 +4 | 逆轉(zhuǎn)陸先生 2026-04-01 | 4/200 |
|
|
[考研] 一志愿9初試366 本雙非求調(diào)劑 +4 | 運氣來得若有似?/a> 2026-04-02 | 4/200 |
|
|
[考研] 材料求調(diào)劑 +10 | 呢呢妮妮 2026-04-01 | 13/650 |
|
|
[考研] 08生物與醫(yī)藥專碩初試346找調(diào)劑 +6 | dianeeee 2026-04-01 | 7/350 |
|
|
[考研] 266分,一志愿電氣工程,本科材料,求材料專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-04-01 | 11/550 |
|
|
[考研] 311求調(diào)劑一志愿合肥工業(yè)大學(xué) +14 | 秋二十二 2026-03-30 | 14/700 |
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +4 | Yam. 2026-03-30 | 5/250 |
|
|
[考研] 0703化學(xué) +20 | 妮妮ninicgb 2026-03-27 | 20/1000 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
[考研] 322求調(diào)劑 +10 | 宋明欣 2026-03-27 | 10/500 |
|
|
[考研] 296求調(diào)劑 +10 | 彼岸t 2026-03-29 | 10/500 |
|
|
[考研] 298求調(diào)劑 +4 | 種圣賜 2026-03-28 | 4/200 |
|
|
[考研] 312,生物學(xué)求調(diào)劑 +3 | 小譯同學(xué)abc 2026-03-28 | 3/150 |
|
|
[考研] 一志愿上海理工能源動力(085800)310分求調(diào)劑 +3 | zhangmingc 2026-03-27 | 4/200 |
|