[問題] 求一個積分的問題已回收
==以下是程式碼==
l=zeros(119,119); x=zeros(1,120); p=zeros(1,120); q=zeros(1,120); c=zeros(1,
119); m = zeros(119); E=zeros(1,119); k=zeros(1,120);
>> for i=1:120 h=pi/120;x(i)=i*h;p(i)=1;k(i)=p(i)/h;
end
for i = 37:84
q(i) = 2.5;
end
for i = 1:118
l(i,i)=k(i)+k(i+1)+q(i)*h;l(i,i+1)=-k(i+1);l(i+1,i)=-k(i+1);
end
l(119,119)=k(119)+k(120)+q(119)*h;
>> for i=1:119
c(i)=1;
end
for i = 1:119
m(i,i)=c(i)*h;
end
>> [V,D]=eig(l,m);
for i = 1:119
E(i)=D(i,i);
end
E(1)
plot(V(:,1))
===以上===
因為要接下畢業學長的研究 這是他弄的MATLAB...
無奈小弟是MATLAB苦手 會用的指令又是基本中的基本...
想請版上的前輩們只點一下迷津 問題如下:
透過大概第六行的 i=37:84
第七行的 q(i)=2.5
那邊套入不同的值 可以PLOT出各種不同的曲線
目前是希望能夠積分出
隨著修改六.七行而出現的各個不同曲線下的曲面面積
估狗 估了半天..嘗試過用int(V(:,1)) 和 quad(V(:,1))
他最後還是會跑圖 可是 沒算積分....
有人有好的建議 或者 指點嗎...先在此跪謝了...
題外話..
有人有推薦 "極新手" 適用且容易上手的MATLAB工具書嗎...
對於MATLAB苦手 常常感到無奈....
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 123.192.14.241