1樓：匿名使用者
我把步長改了改, 因為不然電腦執行太慢不出結果, 所以改長了, 但是沒有其他影響. 下面的**是沒改的.
t=0.1;
r=0.5;
a=2; %d的取值區間上限;
b=0.04; %d的步長, 我取的0.1;
c1=0; %o取值區間下限;
c2=10; %o取值區間上限;
d=0.1; %o的步長, 我取的0.5;
p21=1;r21=0.1;
syms y1;
g=zeros(((c1-c2)/d+1),(a/b*2+1));
n=0;
for o=c1:d:c2
if (o-fix(o)==0)
clear maplemex;
endn=n+1;
m=0;
for x=-a:b:a
m=m+1;
rs=(o^2*(r+2*r21))/(4*(x+0.002*y1)^2+(r+2*r21)^2);
d=4*rs^2+r^2+4*rs*p21+2*p21*r+p21^2;
s1=0.5*(-(2*rs+p21+r)+sqrt(d));
s2=0.5*(-(2*rs+p21+r)-sqrt(d));
r11=rs*(rs+p21)/sqrt(d)*(1/(rs+s2)*(exp(-rs*t)-exp(s2*t))-1/(rs+s1)*(exp(-rs*t)-exp(s1*t)))+exp(-rs*t); %ρ11
r22=rs/sqrt(d)*(exp(s1*t)-exp(s2*t));
r12=0.5*o*(r/2+r21)*(r11-r22)/((x+0.002*y1)^2+(r/2+r21)^2);
f=2.4*10^(-8)*y1.^2.*exp(-4.83*10^(-6).*y1.^2).*r12; %得到被積函式
f=@(y1)subs(f)+0.*y1; % 這裡加了一項0, 就能執行了, 結果肯定不受影響.
g(n,m)=quad(f,0,100);
endend
y=c1:d:c2;
x=-a:b:a;
z=g;
surf(x,y,z);
影象是精簡之後的, 不過你畫出來應該是類似的, 僅僅是更精細:
2樓：
原因是f是你的函式的名字，又是你函式裡變數的名字，導致衝突。你把你的函式改個名字就可以了。
求matlab 大神，幫我改改程式錯誤，我想讓l對t求導，但是總是有錯誤，求幫忙改改，謝謝 10
3樓：我行我素
首先，t沒定義就出現在m=24*t;句中
再者，aa=w1^2*z*sin(nt)-m2^2*z*sin(b2)+m2*x*e2*cos(b2);%曲柄2 角加速度 中，nt間應加*號
求大神看看我電腦配置有救嗎，求大神幫我看看我的電腦配置還有救嗎！！！
樓上的，要麼複製貼上，要麼就是狗不搭八。親你這配置打算玩什麼遊戲呢?cpu有點落伍，不過也不會太差，我的amd4000 七彩虹gt210 64bi位寬，都沒求救。你這配置，滿足吧親，至少有9600gso。cpu可以換成amd 245。那就完美搭配了。 黃啊黃的 你的主機板支援的是ddr2這個淘汰的記...
求大神看看會看的幫我看看手相,求大神幫我看看手相
1,緣主是男生,川子掌代表大膽,自由,放縱。在事業上勇於敢創,會有會成就回。答目前手相上財運並沒有顯示,建議做點固定工作,或者跟隨好的老闆會有好運,不建議自己創業。如果是女生的話,婚姻不好,家裡所有事情要自己擔。易離婚,太強勢 性格不受約束。2 感情線目前顯示還是比較好,記得多珍惜一人感情。一婚到老...
matlab作業急求大神,matlab作業急求大神
r 0.2 syms t phi 15 pi t 位移 x r cos phi r y r sin phi 速度來自 dx diff x dy diff y v sqrt dx dx dy dy v v 加速bai度 a diff v t1 0 4 1500 2 15 x1 subs x,t,t1 ...