1樓:我要那個妾
e的(-1/2(t的平方)的原函式
原函式啥意思?
你的這個表達啊。。。
syms t
f = exp(-1/2.*t.^2);
f = 1/2 +1/sqrt(2*pi)*int(f);
x=0.1:0.1:30;
result = subs(f,t,x);
這個就是答案。
2樓:匿名使用者
你的函式是f(x)還是f(t)啊
3樓:
1. latex 句法很簡單地在網頁中插入科學公式。word裡的equation雖然也能做,但在網頁裡只能貼圖了。
2. 已知函式f(x)是一個定義在某區間的函式,如果存在函式f(x),使得在該區間內的任一點都有
df(x)=f(x)dx,則在該區間內就稱函式f(x)為函式f(x)的原函式。∫f(x)dx = f(x)+c, 即不定積分。
4樓:匿名使用者
用matlab程式設計實現:
x = -2*pi:0.000001:2*pi;
y = (sin(x).*sin(x)+(1/2010)*1./sin(x)).*(cos(x).*cos(x)+(1/2010)*(1./cos(x)));
r = min(y)
r = -12675
r 就是所求的最小值
另外,虛機團上產品**,超級便宜
怎麼用matlab程式設計求(1/sqrt(2*pi))*exp(-t^2/2)在-inf到x上的積分吶?求大神幫忙。。。
5樓:匿名使用者
是這樣嗎:
syms x
int((1/sqrt(2*pi))*exp(-t^2/2),-inf,x)
如何用matlab求exp(x^2)在1到2 的積分?
6樓:計算天下
由於無法求得exp(x^2)的原函式,我們只能用數值演算法來求解,可以用復化梯形公式、romberg公式、gauss公式等,有好多種。我用matlab編了一個用gauss公式求解積分的函式。
function s=gaussintegrate()
%運用gauss求積公式計算數值積分
%f為被積函式,rho為權函式,二者均為符號函式
x=sym('x');
f=exp(x^2);
rho=1;
%a,b分別為求積區間的左界和右界
a=1;
b=2;
%n表示求積結點的個數,是一正整數
n=8;
%本程式利用線性變換將區間[a,b]變換到[-1.1],
%同時令g=f*rho為被積函式,然後利用
%古典的gauss求積公式進行計算,此時直交多項式即為legendre多項式
if n<=0||n~=floor(n)
error('錯誤,n必須是一個非負整數!');
end;
if a>b
error('錯誤,區間的左界a一定不大於右界b!');
end;
%計算n次legendre多項式
syms x;
p=1/(2^n*factorial(n))*diff((x^2-1)^n,n);
w=roots(sym2poly(p));
%計算數值積分
a=zeros(1,n);
s=0;
for k=1:n
a(k)=2/((1-w(k)^2)*(subs(diff(p),w(k))^2));
t=a+(b-a)/2*(w(k)+1);
g=(b-a)/2*subs(f*rho,t);
s=s+a(k)*g;
end;
--------------------------------
我取了8個結點,計算精度就已經達到了小數點後8位,效率還是很高的。
注意:由於matlab呼叫maple的符號計算工具箱,第一次執行時會載入一小會,耐心等待。
以後再執行速度就很快了。
7樓:我要那個妾
syms x
y = exp(x^2);
vpa(int(y,x,1,2)) %int的部分因為無法表達出來,所以matlab用erfi這個誤差函式來代替。vpa是將這個符號數,轉換成有限精度的浮點數。預設字長32。
8樓:我行我素
可按下面計算:
>> syms x
>> a=int(exp(x^2),1,2)a =
-(pi^(1/2)*(erfi(1) - erfi(2)))/2>> vpa(a)
ans =
14.989976019600048615932355867311
9樓:魯凌菁
syms x
isym=vpa(int(exp(x^2),x,0,1))
10樓:行者張強
x ='exp(x^2)';
int(x,1,2)
如果對你有用,請採納最佳!!
迭代法求解方程f(x)=e^x+10*x-2=0的跟。迭代方程為x(k+1)=(2-e^(xk))/10 用matlab程式設計求解
11樓:白楊龍
clc;clear
format long
x0=0;
tol=1e-6;
dx=1;
n=1;
while (abs(dx)>=1e-6) ||抄 (n>1000)x1=(2-exp(x0))/10;
dx=x1-x0;
x0=x1;
n=n+1;
endif n>1000 && (abs(dx)>=1e-6)disp('迭代
失敗')
else
disp('方程的解為')
disp(x1)
disp(['共進行迭代',num2str(n),'次'])end結果:
方程的解為
0.090525117687371
共進行迭代8次
matlab求積分:y=sin(x),t=y^2/(a+b*y),積分範圍(0,pi),
12樓:匿名使用者
首先,你的語
copy句沒問題。
但這個函式對於matlab來說過於複雜,matlab的數值計算能力很強,但是符號計算能力有限。所以他警告:warning:
explicit integral could not be found. (就是說積不出顯性表示式)。
如果可以的話,把a,b換成數字應該是可以積出數值結果的。
如果一定要積出符號解的話,可以用mathematica,或maple它們的符號計算能力比較強。
下面給出mathematica積出的結果(mathematica也積了好久,怪不得matlab積不出來)
13樓:匿名使用者
^可以解。
專>> syms x a b
>> y=sin(x)
y =sin(x)
>> s=int(y^屬2/(b*y+a),0,pi)s =-(-i*a^2*log(-i*a/(a^2-b^2)^(1/2))+i*a^2*log(i*a/(a^2-b^2)^(1/2))+a*(a^2-b^2)^(1/2)*pi+2*a^2*atan(1/(a^2-b^2)^(1/2)*b)-2*b*(a^2-b^2)^(1/2))/b^2/(a^2-b^2)^(1/2)
用matlab求函式在指定點的數值導數
14樓:深眠者
用matlab求函式在指定點的數值導數**如下:
clear;
clc;
syms x
f = sqrt(x^2+1);
diff_f = diff(f);
y1 = subs(diff_f,x,1)y2 = subs(diff_f,x,2)y3 = subs(diff_f,x,3)結果:y1 =0.7071
y2 =0.8944
y3 =0.9487
15樓:匿名使用者
clear;
clc;
syms x
f = sqrt(x^2+1);
diff_f = diff(f);
y1 = subs(diff_f,x,1)y2 = subs(diff_f,x,2)y3 = subs(diff_f,x,3)結果:y1 =
0.7071
y2 =
0.8944
y3 =
0.9487
16樓:傲super福蘭奇
如果得到的是sym格式,只需要新增double 命令就可將sym轉換為double型
這樣就可以看到數字
clear;
clc;
syms x
f = sqrt(x^2+1);
diff_f = diff(f);
y1 = double(subs(diff_f,x,1));
y2 = double(subs(diff_f,x,2));
y3 = double(subs(diff_f,x,3));
結果:y1 =
0.7071
y2 =
0.8944
y3 =
0.9487
17樓:我是幸運小當家
x=1:3
f=inline('sqrt(x.∧2+1)')dx=diff(f([x,4]))
dx=0.8219 0.9262 0.9608
用matlab程式設計計算1加到,用matlab程式設計計算1加到
s 0 for i 1 100 s s i ends 程式的功能是計算 1 2 100 的和,程式用到了一個 for 迴圈結構,基本語句比較簡單。擴充套件資料matlab常用的基本數學函式 abs x 純量的絕對值或向量的長度 angle z 複數z的相角 phase angle sqrt x 開平...
matlab程式設計怎麼求無窮限定積分
用如下程式 syms x f exp x 2 int f,x,0,inf 結果是pi 1 2 2,即2分之根號派 matlab程式設計求無窮限定積分 用如下程式 syms x f exp x 2 int f,x,0,inf 結果是pi 1 2 2,即2分之根號派 matlab中如何計算積分上限為無窮...
matlab中計算定積分時可否將函式作為積分變數
提供兩種解法供參考。1 解析解法作變數置換t 1 x,則積分上限為1,下限為inf syms x t f 1 1 x x x f t subs f,x,1 t int f t,inf,1 ans inf 在2007b和2013a上實測通過。2 數值解法直接對關於x的函式1 x求微分,但得到的函式表示...