 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2205 | 回復(fù): 6 | ||
| 【懸賞金幣】回答本帖問題，作者小蟲飛更高將贈(zèng)送您 5 個(gè)金幣 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
小蟲飛更高新蟲 (初入文壇)
|
[求助]
matlab光強(qiáng)編程，提示矩陣維度 必須一致，不知是哪里出現(xiàn)的問題，代碼在下面 已有1人參與
|
|
|
代碼如下： clear all; k=10; f=1004.8; x=-1.5:0.01:1.5; y=-1.5:0.01:1.5; [x1,y1]=meshgrid(x,y); r=sqrt(x1.^2+y1.^2); Phi=atan(y1./x1); I=(i*k*exp(i*(pi/2+k*f)).*sin(2.*Phi).*quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r.*sin(theta)),2.23,pi)).^2; subplot(2,1,1); mesh(x1,y1,I); axis([-1.5 1.5 -1.5 1.5 0 1]); subplot(2,1,2); subimage(I*255); axis off； 出現(xiàn)錯(cuò)誤如下： 矩陣維度必須一致。 出錯(cuò) formula>@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r.*sin(theta)) 出錯(cuò) quadl (line 62) y = feval(f,x,varargin{:}); y = y(  .';出錯(cuò) formula (line 9) I=(i*k*exp(i*(pi/2+k*f)).*sin(2.*Phi).*quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r.*sin(theta)),2.23,pi)).^2; |
回復(fù)此樓» 猜你喜歡
308求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料調(diào)劑，307分
已經(jīng)有3人回復(fù)
-
085600材料與化工調(diào)劑 280分
已經(jīng)有15人回復(fù)
-
求材料讀博院校
已經(jīng)有8人回復(fù)
-
申博
已經(jīng)有8人回復(fù)
-
復(fù)試調(diào)劑
已經(jīng)有3人回復(fù)
-
中科大材料299求調(diào)劑
已經(jīng)有16人回復(fù)
-
296求調(diào)劑
已經(jīng)有3人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有65人回復(fù)
-
一志愿中科大材料與化工。353分能調(diào)劑到哪？
已經(jīng)有3人回復(fù)
hzlhm
至尊木蟲 (著名寫手)
- 應(yīng)助: 387 (碩士)
- 金幣: 17808.5
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時(shí)
- 蟲號(hào): 1517335
- 注冊(cè): 2011-11-30
- 性別: GG
- 專業(yè): 常微分方程與動(dòng)力系統(tǒng)
5樓2020-09-11 22:54:25
hzlhm
至尊木蟲 (著名寫手)
- 應(yīng)助: 387 (碩士)
- 金幣: 17808.5
- 紅花: 53
- 帖子: 2879
- 在線: 606.3小時(shí)
- 蟲號(hào): 1517335
- 注冊(cè): 2011-11-30
- 性別: GG
- 專業(yè): 常微分方程與動(dòng)力系統(tǒng)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
問題出在，quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r.*sin(theta)),2.23,pi)，這個(gè)數(shù)值積分命令上，該函數(shù)要求點(diǎn)對(duì)點(diǎn)，即已知值為單一數(shù)值（不接受數(shù)組數(shù)據(jù)）。所以你的問題，應(yīng)用for循環(huán)語句來完成。 123.jpg |
» 本帖已獲得的紅花（最新10朵）
2樓2020-09-08 19:50:51
小蟲飛更高
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 20
- 帖子: 48
- 在線: 6.6小時(shí)
- 蟲號(hào): 14223907
- 注冊(cè): 2019-02-18
- 專業(yè): 光學(xué)和光電子學(xué)
送紅花一朵 |
你好，我按照你的編程輸入了下，還是有提示錯(cuò)誤，請(qǐng)你在幫我解答下， clear all; k=10; f=1004.8; x=-0.5:0.01:0.5; y=-0.5:0.01;0.5 [x1,y1]=meshgrid(x,y); r=sqrt(x1.^2+y1.^2); [K,L]=size(r); for m=1:k for n=1:L Phi=atan(y1(n)./x1(m)); r0=r(m,n); eq=quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r0.*sin(theta)),2.23,pi); I0=(i*k*exp(i*(pi/2+k*f))*sin(2*Phi)*eq).^2; I(m,n)=I0; end end Re=real(I);Im=imag(I); s.FaceColor='flat'; %subplot(2,1,1); C=rand(K,L); mesh(x1,y1,real(I).C);hold on 提示錯(cuò)誤如下： formula ans = 0.5000 未定義與 'double' 類型的輸入?yún)?shù)相對(duì)應(yīng)的函數(shù) 'Besselk'。 出錯(cuò) formula>@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r0.*sin(theta)) 出錯(cuò) quadl (line 62) y = feval(f,x,varargin{:}); y = y( .';出錯(cuò) formula (line 13) eq=quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r0.*sin(theta)),2.23,pi); >> |
3樓2020-09-11 11:15:13
小蟲飛更高
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 20
- 帖子: 48
- 在線: 6.6小時(shí)
- 蟲號(hào): 14223907
- 注冊(cè): 2019-02-18
- 專業(yè): 光學(xué)和光電子學(xué)
|
剛我又調(diào)試了下，出現(xiàn)如下錯(cuò)誤，您再幫我看下： k=10; f=1004.8; x=-0.5:0.01:0.5; y=-0.5:0.01:0.5; [x1,y1]=meshgrid(x,y); r=sqrt(x1.^2+y1.^2); [K,L]=size(r); for m=1:k for n=1:L Phi=atan(y1(n)./x1(m)); r0=r(m,n); eq=quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*besselk(2,k.*r0.*sin(theta)),2.23,pi); I0=(i*k*exp(i*(pi/2+k*f))*sin(2*Phi)*eq)^2; I(m,n)=I0; end end Re=real(I);Im=imag(I); s.FaceColor='flat'; %subplot(2,1,1); C=rand(K,L); mesh(x1,y1,real(I).C);hold on axis([-1.5 1.5 -1.5 1.5 0 1]); subplot(2,1,2); subimage(I*255); axis off 提示錯(cuò)誤： 定義變量 "real" 或類 "real"。 出錯(cuò) para (line 22) mesh(x1,y1,real(I).C);hold on |
4樓2020-09-11 15:30:28
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑 +7 | Seele_v5 2026-03-04 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿山東大學(xué)105500藥學(xué)專碩，總分302求調(diào)劑 +5 | 五維天空 2026-03-04 | 11/550 |
|
|
|
[考研] 調(diào)劑 +3 | 13853210211 2026-03-08 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +5 | 好好- 2026-03-07 | 5/250 |
|
|
|
[考研] 286求調(diào)劑 +12 | Faune 2026-03-06 | 14/700 |
|
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | G212 2026-03-03 | 6/300 |
|
|
|
[考研] 物理學(xué)求調(diào)劑 +3 | sx宋 2026-03-05 | 3/150 |
|
|
|
[考研] 材料調(diào)劑 +10 | ounce. 2026-03-04 | 12/600 |
|
|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 4/200 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
|
[考研] 332材料求調(diào)劑 +6 | zjy101327 2026-03-05 | 7/350 |
|
|
|
[考研] 考研282分求調(diào)劑，接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
|
[考研] 一志愿天津大學(xué)085600 319分 材料與化工 金屬方向 求調(diào)劑 +6 | 青科11 2026-03-02 | 6/300 |
|
|
|
[考研] 材料328求調(diào)劑 +10 | 一個(gè)蘿卜02 2026-03-03 | 10/500 |
|
|
|
[考研] 264求調(diào)劑 +8 | 26調(diào)劑 2026-03-03 | 8/400 |
|
|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿中科大080500總分324求調(diào)劑 +3 | jorna 2026-03-03 | 6/300 |
|
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
|
[考研] 0703化學(xué)306調(diào)劑 +4 | 26要上岸 2026-03-03 | 4/200 |
|