1樓:化學工程
這樣表達太不經濟了。改成矩陣形式,用迴圈。【假設a1=1000,具體是多少,你自己改過來】
a=[...
1000 256.3 30.91
979.4 198.9 8.9
926.8 298 18.45
1661 211.8 19.36
1923 180.5 8496
1393 174 20.84
2154 339.4 120.3
1442 121.9 76 ];
syms x
y=0;
for n=1:8
y= y+a(n,1)*exp(-((x-a(n,2))/a(n,3))^2);
endy x=1:1200;
y=subs(y);
plot(x,y)
結果:y =
1000*exp(-(100/3091*x-2330/281)^2)+4897/5*exp(-(10/89*x-1989/89)^2)+4634/5*exp(-(20/369*x-5960/369)^2)+1661*exp(-(25/484*x-5295/484)^2)+1923*exp(-(1/8496*x-361/16992)^2)+1393*exp(-(25/521*x-4350/521)^2)+2154*exp(-(10/1203*x-3394/1203)^2)+1442*exp(-(1/76*x-1219/760)^2)
2樓:匿名使用者
哇,學到了syms 和subs 兩個命令。
a=[...
1000 256.3 30.91
979.4 198.9 8.9
926.8 298 18.45
1661 211.8 19.36
1923 180.5 8496
1393 174 20.84
2154 339.4 120.3
1442 121.9 76 ];
x=1:1200
y=0for n=1:8
y= y+a(n,1)*exp(-((x-a(n,2))/a(n,3)).^2);
endplot(x,y)
3樓:山水遊鷹
你沒給出a1的值,這裡我隨便給了一個值:
a1 = 950 ;
b1 = 256.3 ;
c1 = 30.91 ;
a2 = 979.4 ;
b2 = 198.9 ;
c2 = 8.9 ;
a3 = 926.8 ;
b3 = 298 ;
c3 = 18.45 ;
a4 = 1661 ;
b4 = 211.8 ;
c4 = 19.36 ;
a5 = 1923 ;
b5 = 180.5 ;
c5 = 8496 ;
a6 = 1393 ;
b6 = 174 ;
c6 = 20.84 ;
a7 = 2154 ;
b7 = 339.4 ;
c7 = 120.3 ;
a8 = 1442 ;
b8 = 121.9 ;
c8 = 76 ;
x = 1:1200;
y = a1*exp(-((x-b1)/c1).^2) + a2*exp(-((x-b2)/c2).^2) + ...
a3*exp(-((x-b3)/c3).^2) + a4*exp(-((x-b4)/c4).^2) + ...
a5*exp(-((x-b5)/c5).^2) + a6*exp(-((x-b6)/c6).^2) + ...
a7*exp(-((x-b7)/c7).^2) + a8*exp(-((x-b8)/c8).^2) ;
plot(x,y)
matlab畫圖問題
x0 1000 1500 2000 2500 3000 3500 3800 4000 y 135.33 147.10 152.98 156.91 147.10 138.27 133.37 125.53 x0太大,以x的冪作為基函式會導致設計矩陣尺度太差,列變數幾乎線性相依。變換為 1.5 1.5 範...
matlab ezplot畫圖,matlab中關於ezplot畫圖問題
世界屋脊吧 f的取值範圍呢,你還是想畫關於f,b的影象。如果是f關於b的,如下 f1 0.5.0.5 0.05 b f 4 0.5.0.5 0.05 b f 4 0.25 f2 1.1 0.05 b f f3 diff f1,f 1 你上面這一行有誤吧 f4 1.1 0.05 b f 2 f5 1....
matlab如何描點畫圖,matlab 如何描點畫圖
智者總要千慮 x 5 5 40 y 0.02 0.04 0.058 0.081 0.103 0.121 0.135 0.154 plot x,y 導超 分析,首先畫出離散點,再進行擬合。用直線擬合。x 5 10 15 20 25 30 35 40 y 0.02 0.04 0.058 0.081 0....