第一題
>> p=[133 0 122 0 0 1] % p=133x^5+122x^3+1
p =
133 0 122 0 0 1
>> q=[2 0 100 0 1] % q=2x^4+100x^2+1
q =
2 0 100 0 1
>> M=conv(p,q)
M =
Columns 1 through 9
266 0 13544 0 12333 2 122 100 0
Column 10
1
>> % p, q 兩多項式的乘積為 266x^9+13544x^7+12333x^5+2x^4+122x^3+100x^2+1
執行結果
第二題
>> p=[100 23 1 45];x=magic(5)
x =
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9
以magic(5),x=每個元素帶入f(x)的答案為
>> f=polyval(p,x)
f =
498009 1395717 169 52725 342735
1228935 13125 35479 278967 415549
6817 22479 223645 809265 1075999
102355 176169 694267 936309 2955
135939 590715 1576945 939 74817
執行結果
第三題
>> p=[133 0 122 0 0 1] % p=133x^5+122x^3+1
p =
133 0 122 0 0 1
>> q=[2 0 100 0 1] % q=2x^4+100x^2+1
q =
2 0 100 0 1
>> [s,r]=deconv(p,q)
s =
66.5000 0
r =
1.0e+003 *
0 0 -6.5280 0 -0.0665 0.0010
>> % 商數s=66.5x 餘數r=e^3*(-6.528)x^3+e^3*(-0.0665)x+e^3*(0.001)
執行結果
第四題
>> p=[133 0 122 0 0 1] % p=133x^5+122x^3+1
p =
133 0 122 0 0 1
>> roots(p)
ans =
-0.0045 + 0.9578i
-0.0045 - 0.9578i
0.1039 + 0.1744i
0.1039 - 0.1744i
-0.1988
>> q=[2 0 100 0 1] % q=2x^4+100x^2+1
q =
2 0 100 0 1
>> roots(q)
ans =
0 + 7.0704i
0 - 7.0704i
0 + 0.1000i
0 - 0.1000i
執行結果
第五題
>> x=[1:9]; y=[1210, 1866, 2301, 2564, 2724, 2881, 2879, 2915, 3010];
>> p=polyfit(x,y,3)
p =
6.3047 -134.4603 994.3540 350.9127
>> f=@(a) (6.3047*a^3-134.4603*a^2+994.354*a+350.9127)
f =
@(a) (6.3047*a^3-134.4603*a^2+994.354*a+350.9127)
>> fplot(f,x)
>> hold on
>> plot(x,y)
>>
此題執行後，其polyfit的結果似乎只有x=1~2的時候比較符合，因為不管我限制多少範圍
做出來的圖ｘ座標都是從１～２
然後將其和原圖放在一起
執行結果
0 意見:
張貼留言