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間應加*號
求大神幫我看一下這個matlab程式怎麼改
yyy名 quadgk比quad好用一些 clear all clc fi 5.pi.180 f 25400 e 5 vx 10 w1 370.2.pi.60 w2 100.2.pi.60 r0 50.5 h 0.0005 a 0.00005 振幅alb 0.00001 abc 1 25.4 f1 ...
求大神幫我看看我電腦配置，求大神幫我看看我的電腦配置，求個建議，謝謝！
我說 你這跑分也太假了吧 a4的cpu最多跑1700分440的顯示卡也就500分 1333的記憶體不會超過300分 硬碟更扯不是sdd居然有700分 你這硬碟也就120分孩子你不會是開掛了吧 你這配置跑分不會超過3000分建議你換一個cpu 基於你這主機板只能上a8 5600k了520塊記憶體換成2...
求大神幫我看看這兩個玉鐲，求大神幫我看看這兩個玉鐲 150
效果看，是糯冰種的翡翠扁條手鐲，飄碎花。底子還行，棉多，水頭中檔。第一張看，還有一點黃翡色花，對整體影響不大。總的看，算中檔貨。各地 都不同，通常三五千元都是正常的 範圍。 lm艾瑪 我也不懂玉，但前段時間朋友給我帶了款玉。水色各方面和您 上這款挺像，價位三千多，說是進價。這件物品符合天然翡翠玻璃光...