| 5 | 1/1 | 返回列表 |
| 查看: 2064 | 回復: 10 | |||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||
bluethky木蟲 (著名寫手)
|
[求助]
matlab程序錯誤,提示如下【弄了半夜還是不行】 已有2人參與
|
||
|
function [x,y,Y]=levydistrintheorytest(t0,lambda,T,a,b); % 計算理論的二維levy分布的圖象 畫出來的3維曲面圖 % t0 :levy walk 趨于零的時間 % lambda :冪律值 % T:時間 % a: 離散的大步長的最大值 % b:小步長的最大值 error=0.125/180*pi; %假定平行于坐標軸的誤差角度 syms u c=-gamma(2-lambda)*t0^(lambda-1)*cos(pi*lambda/2)/lambda; h1=20; %大的離散步長 h2=0.1; %小的離散步長 x1=-a:h1:-b; x2=-(b-h2):h2 b-h2);x3=b:h1:a; x=[x1,x2,x3]; disp('%%%%%%%%%%%計算的次數(shù)為%%%%%%%%%%%%%%%') N=length(x) y=x; Z=zeros(N,N); for i=1:1 %x坐標 i for j=1:N % y坐標 theta=angle(complex(x(i),y(j))); if theta<error && theta>-error % x軸正半軸 Z(j,i)=double(int(cos(u.*x(i)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./pi+cos(1/u.*x(i)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1)); elseif theta<pi/2+error && theta>pi/2-error % y軸正半軸 Z(j,i)=double(int(cos(u.*y(j)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./pi+cos(1/u.*y(j)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1)); elseif theta>pi-error || theta<-pi+error %x軸負半軸 Z(j,i)=double(int(cos(u.*x(i)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./pi+cos(1/u.*x(i)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1)); elseif theta>-pi/2-error && theta<-pi/2+error %y軸負半軸 Z(j,i)=double(int(cos(u.*y(j)/(c*T)^(1/lambda))*exp(-(abs(u))^lambda)./pi+cos(1/u.*y(j)/(c*T)^(1/lambda))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1)); else Z(j,i)=double(int(cos(u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-(abs(u))^lambda)./pi+cos(1/u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1))*(double(int(cos(u.*y(j)/(c*T)^(1/lambda)/cos(theta))*exp(-(abs(u))^lambda)./pi+cos(1/u.*y(j)/(c*T)^(1/lambda)/cos(theta))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1))); %zy=(double(int(cos(u.*y(j)/(c*T)^(1/lambda)/cos(theta))*exp(-(abs(u))^lambda)./pi+cos(1/u.*y(j)/(c*T)^(1/lambda)/cos(theta))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1))); %Z(j,i)=zx*zy; end end end Y=Z; end 錯誤提示; Error using mupadmex Error in MuPAD command: DOUBLE cannot convert the input expression into a double array. If the input expression contains a symbolic variable, use the VPA function instead. Error in sym/double (line 710) Xstr = mupadmex('symobj::double', S.s, 0); Error in levydistrintheorytest (line 37) Z(j,i)=double(int(cos(u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-(abs(u))^lambda)./pi+cos(1/u.*x(i)/(c*T)^(1/lambda)/sin(theta))*exp(-1/(abs(u))^lambda)/u^2./pi,0,1))*(double(int(cos(u.*y( |
木蟲 (著名寫手)
木蟲 (著名寫手)
專家顧問 (職業(yè)作家)
![]() |
專家經(jīng)驗: +831 |

木蟲 (著名寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 總分322求生物學/生化與分子/生物信息學相關(guān)調(diào)劑 +6 | 星沉uu 2026-03-26 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 282求調(diào)劑 不挑專業(yè) 求收留 +3 | Yam. 2026-03-30 | 4/200 |
|
|
[考研] 0817化工學碩調(diào)劑 +7 | 努力上岸中! 2026-03-31 | 7/350 |
|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +3 | 鐵憨憨123425 2026-03-31 | 3/150 |
|
|
[考研] 一志愿 南京航空航天大學 ,080500材料科學與工程學碩 +5 | @taotao 2026-03-30 | 5/250 |
|
|
[考研] 看看我!301求調(diào)劑 081700 +3 | 一個小明同學 2026-03-24 | 3/150 |
|
|
[考研] 材料與化工調(diào)劑一志愿大連海事085600,349 +6 | 吃的不少 2026-03-30 | 6/300 |
|
|
[考研] 英一數(shù)一408,總分284,二戰(zhàn)真誠求調(diào)劑 +3 | 12.27 2026-03-30 | 5/250 |
|
|
[考研] 一志愿中海洋材料357 +3 | 麥恩莉. 2026-03-30 | 3/150 |
|
|
[考研]
|
鐘llll 2026-03-26 | 5/250 |
|
|
[考研] 085600 286分 材料求調(diào)劑 +11 | 麻辣魷魚 2026-03-27 | 12/600 |
|
|
[考研] 309求調(diào)劑 +15 | 誰不是少年 2026-03-29 | 15/750 |
|
|
[碩博家園] 求調(diào)劑 有機化學考研356分 +10 | Nadiums 2026-03-25 | 11/550 |
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
[考研] 2026年華南師范大學歡迎化學,化工,生物,生醫(yī)工等專業(yè)優(yōu)秀學子加入! +3 | llss0711 2026-03-28 | 6/300 |
|
|
[考研] 0856材料化工調(diào)劑 總分330 +14 | zhubinhao 2026-03-27 | 14/700 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 339求調(diào)劑,想調(diào)回江蘇 +6 | 烤麥芽 2026-03-27 | 8/400 |
|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 303求調(diào)劑 +6 | 藍山月 2026-03-25 | 6/300 |
|