| 查看: 5437 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研]
|
是辰啊 2026-04-02 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | Aniyaio 2026-04-02 | 7/350 |
|
|
[考研] 農(nóng)學(xué)考研求調(diào)劑 +3 | dkdkxm 2026-04-01 | 3/150 |
|
|
[考研] 學(xué)碩化學(xué)工程與技術(shù),一志愿中國海洋大學(xué)320+求調(diào)劑 +8 | 披星河 2026-04-02 | 8/400 |
|
|
[考研] 一志愿北京科技大學(xué)材料學(xué)碩328分求調(diào)劑 +6 | 1段時間 2026-03-31 | 7/350 |
|
|
[考研] 266分,一志愿電氣工程,本科材料,求材料專業(yè)調(diào)劑 +4 | 哇呼哼呼哼 2026-04-02 | 4/200 |
|
|
[考研] 找調(diào)劑 +6 | 楚喬喬 2026-04-01 | 6/300 |
|
|
[考研] 0856材料與化工調(diào)劑,339 +13 | 10213207 2026-03-31 | 13/650 |
|
|
[考研] 324求調(diào)劑 +5 | 想上學(xué)求調(diào) 2026-04-01 | 6/300 |
|
|
[考研] 303求調(diào)劑 +8 | DLkz1314. 2026-03-30 | 8/400 |
|
|
[考研] 材料調(diào)劑 +11 | 一樣YWY 2026-03-31 | 11/550 |
|
|
[考研] 材料調(diào)劑 +10 | Eujd1 2026-03-31 | 11/550 |
|
|
[考研] 086000生物與醫(yī)藥298調(diào)劑求助 +4 | 元元青青 2026-03-31 | 6/300 |
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
[考研] 343求調(diào)劑 +8 | 愛羈絆 2026-03-28 | 8/400 |
|
|
[考研] 274求調(diào)劑 +6 | xiao愛同學(xué) 2026-03-30 | 6/300 |
|
|
[考研] 0703化學(xué)321分求調(diào)劑 +10 | 三dd. 2026-03-30 | 11/550 |
|
|
[考研] 085600,專業(yè)課化工原理,321分求調(diào)劑 +5 | 大饞小子 2026-03-28 | 5/250 |
|
|
[考研] 復(fù)試調(diào)劑 +3 | raojunqi0129 2026-03-28 | 3/150 |
|
|
[考研] 295求調(diào)劑 +5 | 1428151015 2026-03-27 | 6/300 |
|