版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進(jìn)行手機(jī)號驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

【懸賞金幣】回答本帖問題，作者小蟲飛更高將贈(zèng)送您 5 個(gè)金幣

小蟲飛更高

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 20
帖子: 48
在線: 6.6小時(shí)
蟲號: 14223907
注冊: 2019-02-18
專業(yè): 光學(xué)和光電子學(xué)

[求助] 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ù)此樓

» 猜你喜歡

A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向已經(jīng)有59人回復(fù)
0817找調(diào)劑已經(jīng)有5人回復(fù)
290 材料與化工求調(diào)劑已經(jīng)有3人回復(fù)
調(diào)劑的同學(xué)，走過路過，不要錯(cuò)過已經(jīng)有15人回復(fù)
化工270求調(diào)劑已經(jīng)有11人回復(fù)
考研求調(diào)劑，工科，總分304 已經(jīng)有4人回復(fù)
新疆大學(xué)地質(zhì)與礦業(yè)工程學(xué)院招生已經(jīng)有22人回復(fù)
交大前校長王樹國：現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室，而是在企業(yè)研究院已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有14人回復(fù)
0703求調(diào)劑已經(jīng)有8人回復(fù)

1樓 2020-09-08 11:19:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hzlhm

至尊木蟲 (著名寫手)

應(yīng)助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時(shí)
蟲號: 1517335
注冊: 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)對點(diǎn)，即已知值為單一數(shù)值（不接受數(shù)組數(shù)據(jù)）。所以你的問題，應(yīng)用for循環(huán)語句來完成。

matlab光強(qiáng)編程，提示矩陣維度必須一致，不知是哪里出現(xiàn)的問題，代碼在下面

123.jpg

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

小蟲飛更高

QQ：2120156492

2樓2020-09-08 19:50:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

小蟲飛更高

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 20
帖子: 48
在線: 6.6小時(shí)
蟲號: 14223907
注冊: 2019-02-18
專業(yè): 光學(xué)和光電子學(xué)

送紅花一朵

引用回帖:

2樓: Originally posted by hzlhm at 2020-09-08 19:50:51
問題出在，quadl(@(theta)(1+cos(theta))./(1-cos(theta)).*sin(theta).*Besselk(2,k.*r.*sin(theta)),2.23,pi)，這個(gè)數(shù)值積分命令上，該函數(shù)要求點(diǎn)對點(diǎn)，即已知值為單一數(shù)值（不接受數(shù)組數(shù)據(jù)）。所以你的問題，應(yīng)用 ...

你好，我按照你的編程輸入了下，還是有提示錯(cuò)誤，請你在幫我解答下，
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ù)相對應(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);
>>

贊一下

回復(fù)此樓

3樓2020-09-11 11:15:13

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

小蟲飛更高

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 20
帖子: 48
在線: 6.6小時(shí)
蟲號: 14223907
注冊: 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

贊一下

回復(fù)此樓

4樓2020-09-11 15:30:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

hzlhm

至尊木蟲 (著名寫手)

應(yīng)助: 387 (碩士)
金幣: 17808.5
紅花: 53
帖子: 2879
在線: 606.3小時(shí)
蟲號: 1517335
注冊: 2011-11-30
性別: GG
專業(yè): 常微分方程與動(dòng)力系統(tǒng)

【答案】應(yīng)助回帖

引用回帖:

4樓: Originally posted by 小蟲飛更高 at 2020-09-11 15:30:28
剛我又調(diào)試了下，出現(xiàn)如下錯(cuò)誤，您再幫我看下：
k=10;
f=1004.8;
x=-0.5:0.01:0.5;
y=-0.5:0.01:0.5;
=meshgrid(x,y);
r=sqrt(x1.^2+y1.^2);
=size(r);
for m=1:k
for n=1:L
Phi=atan(y1(n)./x1(m) ...

real是函數(shù)，應(yīng)該不存在定義變量的問題�？赡苁擒浖膯栴}，重啟試一試看。

贊一下

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

小蟲飛更高

QQ：2120156492

5樓2020-09-11 22:54:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

小蟲飛更高

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 20
帖子: 48
在線: 6.6小時(shí)
蟲號: 14223907
注冊: 2019-02-18
專業(yè): 光學(xué)和光電子學(xué)

送紅花一朵

引用回帖:

5樓: Originally posted by hzlhm at 2020-09-11 22:54:25
real是函數(shù)，應(yīng)該不存在定義變量的問題�？赡苁擒浖膯栴}，重啟試一試看。...

我試了試其他筆記本的電腦上的版本matlab，還是提示未定義變量‘real ’或類‘real’,請您運(yùn)行下我的程序，看看是什么問題呢？

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : lens.m

2020-09-13 12:38:12, 569 bytes

6樓2020-09-13 12:38:15

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

小蟲飛更高

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 20
帖子: 48
在線: 6.6小時(shí)
蟲號: 14223907
注冊: 2019-02-18
專業(yè): 光學(xué)和光電子學(xué)

引用回帖:

5樓: Originally posted by hzlhm at 2020-09-11 22:54:25
real是函數(shù)，應(yīng)該不存在定義變量的問題。可能是軟件的問題，重啟試一試看。...

定義變量的問題解決了，但是又遇到了一個(gè)索引函數(shù)編制的問題，問題在代碼后面：
clear all;
close all;
clc;
k=10;
f=1004.8;
x=-0.5:0.01:0.5;
y=-0.5:0.01:0.5;
[x1,y1]=meshgrid(x,y);%直角坐標(biāo)轉(zhuǎn)化為極坐標(biāo)
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ò)誤: 不能使用 {} 或 . 索引為函數(shù)編制索引。

出錯(cuò) para (line 25)
mesh(x1,y1,real(I).C);

贊一下

回復(fù)此樓

7樓2020-09-13 20:38:56

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主小蟲飛更高的主題更新

返回列表

不應(yīng)助 確定回帖應(yīng)助 (注意：應(yīng)助才可能被獎(jiǎng)勵(lì)，但不允許灌水，必須填寫15個(gè)字符以上)

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[教師之家] 交大前校長王樹國：現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室，而是在企業(yè)研究院 +3	zju2000 2026-03-08	4/200	2026-03-08 14:55 by Equinoxhua
[考研] 0856材料與化工290求調(diào)劑 +3	Nebulala 2026-03-08	4/200	2026-03-08 13:49 by 30660438
[考研] 085700資環(huán)求調(diào)劑，初始279，六級已過，英語能力強(qiáng) +4	085700資環(huán)調(diào)劑 2026-03-03	5/250	2026-03-08 13:35 by 30660438
[考研] 08工科 +4	li李樂成 2026-03-06	4/200	2026-03-07 21:46 by 韓應(yīng)許
[考研] 267化工調(diào)劑求助 +7	聰少OZ 2026-03-04	7/350	2026-03-07 17:33 by 無際的草原
[考博] 26申博-目前4篇SCI一作 +4	chen_2024 2026-03-02	4/200	2026-03-05 22:59 by sicilyl0001
[考研] 求調(diào)劑，學(xué)校研究所都可以，材料與化工267分 +6	wmx1 2026-03-05	6/300	2026-03-05 20:01 by houyaoxu
[考研] 求材料調(diào)劑 +4	berdmond 2026-03-05	4/200	2026-03-05 19:45 by 黑衣饅頭人
[考研] 304求調(diào)劑 +4	曼殊2266 2026-03-05	4/200	2026-03-05 17:10 by zhukairuo
[考研] 0856材料與化工，270求調(diào)劑 +17	YXCT 2026-03-01	20/1000	2026-03-05 09:22 by 一切OK
[考研] 0855機(jī)械工程324求調(diào)劑 +3	鄭SK 2026-03-02	7/350	2026-03-04 17:24 by congmingdou
[考研] 295求調(diào)劑 +6	等春來， 2026-03-04	6/300	2026-03-04 16:11 by zhukairuo
[考研] 能動(dòng)297求調(diào)劑，本科川大 +4	邵11 2026-03-04	4/200	2026-03-04 15:29 by lylylyheart
[考研] 331求調(diào)劑 +3	zzZ&zZ 2026-03-03	3/150	2026-03-04 10:37 by aaa顯卡批發(fā)王總
[考研] 292求調(diào)劑 +3	sgbl 2026-03-03	3/150	2026-03-03 17:30 by barlinike
[考研] 調(diào)劑 +5	13853210211 2026-03-02	7/350	2026-03-03 11:42 by 13853210211
[考研] 化工京區(qū)271求調(diào)劑 +7	11ing 2026-03-02	7/350	2026-03-03 07:30 by 利好利好.
[考研] 求調(diào)劑 +11	yunziaaaaa 2026-03-01	13/650	2026-03-02 21:59 by sunny81
[考研] 291 求調(diào)劑 +3	化工2026屆畢業(yè)?/a> 2026-03-02	3/150	2026-03-02 12:55 by houyaoxu
[基金申請] 此成果不能導(dǎo)入原因：元數(shù)據(jù)必填信息不完整，可進(jìn)行補(bǔ)充。 +4	Kittylucky 2026-03-02	5/250	2026-03-02 11:07 by jurkat.1640

24小時(shí)熱門版塊排行榜

[求助] matlab光強(qiáng)編程，提示矩陣維度 必須一致，不知是哪里出現(xiàn)的問題，代碼在下面 已有1人參與

» 猜你喜歡

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

【答案】應(yīng)助回帖

» 本帖已獲得的紅花（最新10朵）

» 本帖附件資源列表

[求助] matlab光強(qiáng)編程，提示矩陣維度必須一致，不知是哪里出現(xiàn)的問題，代碼在下面已有1人參與