2012-10-17 61 views
1

我是新來的MATLAB和不知道如何使這個多項式函數的值與W = 1 /(X^2 + 1)如何使用以下值在matlab中構建表格?

我嘗試現在

x= [-5,-3,-1, 1, 3,5] 

x = 

-5 -3 -1  1  3  5 

y= [0.0385, 0.10, 0.50, 0.50, 0.10, 0.0385] 

y = 

0.0385 0.1000 0.5000 0.5000 0.1000 0.0385 

yp=[0.0148,0.06,0.50,-0.50,-0.6,-0.0148] 

yp = 

0.0148 0.0600 0.5000 -0.5000 -0.6000 -0.0148 

hp = hermite (x, y, yp) 

hp = 

-0.0000 -0.0000 0.0001 0.0004 -0.0011 -0.0100 0.0072 0.0969 -0.0113 -0.4156 0.0051 0.8282 

比較表它只是在表中比較hp和w的值。

有人能幫我嗎?

感謝您的幫助

+0

請不要再使用[作業]標籤;它正在被刪除的過程中。 – geoffspear

+0

感謝您的信息 –

回答

0

我看不出有什麼w在你的問題,但總體上表達

hp == w 

會,如果向量具有相同的長度,返回0 S和1矢量s,表示兩個向量匹配(1)或不匹配(0)的情況。順帶一提,請注意浮點數的平等是比較充滿了「問題」,你可能會更好評價:

abs(hp-w) < 10^-6 

通過您的首選寬容更換10^-6

鑑於你定義爲w,你應該能夠編寫

hp == 1/(x.^2+1) 

注意使用的elementwise平方操作.^,該方法返回相同長度的向量作爲x與每個元素對應的平方元素在x。當然,表達式

hp - 1./(x.^2+1) 

將返回差異的向量,這可能是你想要的。

+0

w等於1 /(x^2 + 1) –

+0

您還需要一個'.'以及'1 ./(x。^ 2 + 1)'。另外,我認爲在OP的情況下'hp'與'x'的大小不同。至少從帖子中的數據來看似乎如此。 – angainor

+0

感謝@angainor,我更新了最後一個代碼片段。我將離開OP來整理不匹配的長度,他應該能夠弄清楚。 –