 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2043 | 回復: 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( |
回復此樓» 猜你喜歡
招收博士1-2人
已經(jīng)有4人回復
-
08工科 320總分 求調(diào)劑
已經(jīng)有7人回復
-
生物學一志愿985，分數(shù)349求調(diào)劑
已經(jīng)有7人回復
-
0856材料專碩353求調(diào)劑
已經(jīng)有4人回復
-
山東省面上項目限額評審
已經(jīng)有4人回復
-
一志愿北京化工大學070300 學碩336求調(diào)劑
已經(jīng)有5人回復
-
生物學071000 329分求調(diào)劑
已經(jīng)有4人回復
-
一志愿華中科技大學071000，求調(diào)劑
已經(jīng)有4人回復
-
生物學調(diào)劑
已經(jīng)有4人回復
-
求調(diào)劑院校信息
已經(jīng)有4人回復
» 本主題相關價值貼推薦，對您同樣有幫助:
- matlab 運行ode45出錯，不知道什么原因
已經(jīng)有4人回復
- MATLAB解方程提示錯誤‘is not a valid expression equation’
已經(jīng)有10人回復
- matlab為啥出現(xiàn)以下這種錯誤
已經(jīng)有6人回復
- 以下matlab代碼出現(xiàn)問題
已經(jīng)有10人回復
- 請大俠指教matlab的這個錯誤提示到底是個神馬意思？
已經(jīng)有12人回復
- 用matlab求解非線性方程組說無解，一定是方程組本身無解，還是有可能程序有問題呢？
已經(jīng)有11人回復
- matlab的m文件關聯(lián)問題
已經(jīng)有17人回復
- 請matlab高手看一下這個問題怎樣編程才好，附程序代碼，謝謝！
已經(jīng)有9人回復
- matlab-出錯提示Undefined method for input arguments type 'double'.
已經(jīng)有5人回復
- matlab語句問題
已經(jīng)有5人回復
- 求助一個matlab關于分類討論的問題
已經(jīng)有9人回復
- 大家愛幫我看看這段MATLAB程序為啥出錯了
已經(jīng)有5人回復
- 問下這個問題怎么解決.....真心啊
已經(jīng)有4人回復
- Matlab兩個小問題
已經(jīng)有10人回復
- 同樣的MATLAB程序，為何在別人的電腦上運行不了
已經(jīng)有21人回復
- VS2005下matlab調(diào)用
已經(jīng)有8人回復
- 【求助】用MATLAB畫三維圖程序編寫的問題
已經(jīng)有6人回復
- 【求助】程序出錯
已經(jīng)有4人回復
wurongjun
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +831 - 數(shù)學EPI: 9
- 應助: 791 (博后)
- 貴賓: 0.308
- 金幣: 24609
- 散金: 310
- 紅花: 75
- 帖子: 3004
- 在線: 881.4小時
- 蟲號: 1368482
- 注冊: 2011-08-14
- 性別: GG
- 專業(yè): 計算數(shù)學與科學工程計算
- 管轄: 數(shù)學
3樓2014-05-10 06:53:34
bluethky
木蟲 (著名寫手)
- 應助: 1 (幼兒園)
- 金幣: 4698.8
- 散金: 849
- 紅花: 11
- 帖子: 1398
- 在線: 227.7小時
- 蟲號: 2609601
- 注冊: 2013-08-21
- 專業(yè): 植物學研究的新技術、新方
2樓2014-05-10 01:35:55
bluethky
木蟲 (著名寫手)
- 應助: 1 (幼兒園)
- 金幣: 4698.8
- 散金: 849
- 紅花: 11
- 帖子: 1398
- 在線: 227.7小時
- 蟲號: 2609601
- 注冊: 2013-08-21
- 專業(yè): 植物學研究的新技術、新方
4樓2014-05-10 10:03:23
bluethky
木蟲 (著名寫手)
- 應助: 1 (幼兒園)
- 金幣: 4698.8
- 散金: 849
- 紅花: 11
- 帖子: 1398
- 在線: 227.7小時
- 蟲號: 2609601
- 注冊: 2013-08-21
- 專業(yè): 植物學研究的新技術、新方
5樓2014-05-10 10:14:42
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
|
[考研] 化學工程321分求調(diào)劑 +18 | 大米飯！ 2026-03-15 | 22/1100 |
|
|
|
[考研] 326求調(diào)劑 +5 | 諾貝爾化學獎覬?/a> 2026-03-15 | 8/400 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 332求調(diào)劑 +3 | 鳳凰院丁真 2026-03-20 | 3/150 |
|
|
|
[考研] 機械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學大學 ，080500材料科學與工程學碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 求調(diào)劑一志愿南京航空航天大學289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 08工學調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 材料學碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|